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Abstract 

(N' 

This review describes in detail the essential techniques used in microscopic theories 

on spintronics. We have investigated the domain wall dynamics induced by electric 
00 . current based on the s-d exchange model. The domain wall is treated as rigid and 

planar and is described by two collective coordinates: the position and angle of wall 
C"-— ■ magnetization. The effect of conduction electrons on the domain wall dynamics 

is calculated in the case of slowly varying spin structure (close to the adiabatic 
^, limit) by use of a gauge transformation. The spin-transfer torque and force on the 

wall are expressed by Feynman diagrams and calculated systematically using non- 
equilibrium Green's functions, treating electrons fully quantum mechanically. The 
wall dynamics is discussed based on two coupled equations of motion derived for 
two collective coordinates. The force is related to electron transport properties, 
resistivity, and the Hall effect. Effect of conduction electron spin relaxation on the 
torque and wall dynamics is also studied. 
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Table 1 

List of important variables and their unit. Equation or equation number indicates 

def inition. 

S local spin vector 

n local spin direction S/S 

no local spin for domain wall configuration 

J exchange interaction between local spins 

K easy axis energy gain of local spin per spin 

K± hard axis energy loss of local spin 

J S( i s-d exchange interaction 

a Gilbert damping parameter of local spin 

(3 sr (3 arising from electron spin relaxation 

(3 W effective (3 or force acting on domain wall 

A domain wall thickness 

X center position of domain wall 

4>o collective angle out of easy plane of domain wall Eq. (82) 

iV w number of spins in domain wall 2XA/a 3 

M w mass of domain wall Eq. (96) kg 

R w resistance due to domain wall Q 

p s resistivity due to spin structure Eq. (234) 0, m 

Rs resistance due to spin structure p s L/A 

1 Introduction 

1.1 Magneto- electric effects and devices 



Present information technology is based on electron transport and magnetism. 
Magnetism has been most successful in high-density storages such as hard 
disks. For integration of magnetic storages into electronic circuits, mechanisms 
are necessary to convert electric current /volt age into magnetic information 
and vice versa. The most common and oldest electro-magnetic coupling is the 
one arising from Maxwell's equations. Ampere's law or Oersted's law, discov- 
ered in the early nineteenth century, describes the magnetic field created by an 
electric current (Fig. 1). This field can be applied to write information in mag- 
netic information storages. In fact, this mechanism is so far the only successful 



e electron charge (e < 0) 

j current density 

j s spin current density (without spin length |) 

j c threshold current density 

f c intrinsic threshold current density 
v electron density of states at Fermi energy per volume 

s electron spin density (without spin length |) 

s electron spin density vector 

s spin density vector in gauge transformed frame 

r electron lifetime (spin-dependent in §14) 

P spin polarization of current 

M spin polarization of conduction electron 

v s drift velocity of electron spin 

v c v s at intrinsic threshold 
kpo Fermi wavelength of conduction electron with spin a = : 

€f Fermi energy of conduction electron 

cfu Fermi energy with spin splitting included (a = ±) 

a lattice constant 

L system length along z direction 

A cross sectional area of system m 2 

mechanism used in commercial high-density magnetic devices. On the other 
hand, Faraday's law provides us means to convert magnetic information into 
electric current or voltage, for instance, detecting magnetic information by 
scanning a read head (a coil) on the stored magnetic bits. This mechanism is 
not, however, very useful in high density storages, and various magnetoresistive 
effects based on solid-state systems have been discovered and applied in the 
late twentieth century, such as anisotropic magnetoresistance (AMR), giant 
magnetoresistance (GMR), and tunneling magnetoresistance (TMR) effects 
(Fig. 2). AMR is a resistivity dependent on the angle between the magnetiza- 
tion and the electric current, discovered in 1857[1]. It arises from the coupling 
of magnetization and electrons' orbital motion due to spin-orbit interaction 
[2]. The resistivity change is of the order of only a few percent, but AMR is 
more efficient than using Faraday induction used in magnetic tape and hard 
disks in early days. Magnetic heads with higher sensitivity were developed by 
use of the GMR effect in thin magnetic multilayers discovered in 1988 [3,4]. In 
such multilayers, a strong magnetization dependence of the resistivity arises 
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from the spin-dependent scattering of electrons at the interface between a thin 
ferromagnetic layer and nonmagnetic metallic layers. A. Fert and P. Griinberg 
were awarded the Nobel Prize in 2007 for the discovery of the GMR effect. 
Quite recently GMR heads are being replaced by even more efficient TMR 
heads, where the nonmagnetic layer is replaced by an insulating barrier [5-7]. 
These rapid developments of read-out mechanisms by use of solid state sys- 
tems have made possible so far the rapid increase of recording density. These 
magnetoresistances are due to the exchange interaction between localized spin 
and conduction electrons, arising from the overlap of electron wave functions 
and their correlation. Present magnetic devices are therefore one of the most 
successful outcomes of material science. 
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Fig. 1. Ampere and Faraday's laws. 
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Fig. 2. GMR and TMR systems that realize large magnetoresistance using s-d in- 
teraction between localized spin and conduction electron. 



1.2 Magnetization switching by s-d exchange interaction 



Electron transport in magnetic metals and semiconductors is modeled by the 
so-called s-d model, where the conduction and magnetization degrees of free- 
dom are separated from each other. The conduction electrons we consider are 
non-interacting with each other, but are scattered by spin-independent im- 
purities and also by spin-dependent impurities (resulting in spin relaxation). 
The localized spin at position x at time t is described by a variable S(x,t). 
The localized spin is related to magnetization as 
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where jjl b = ^- is the Bohr magneton, g = 2 is the g-factor, and 7 = ^-{> 0) 
is the gyromagnetic ratio. The electron charge e is negative. In this paper S 
is treated as a classical variable, since quantum fluctuation of S is blocked 
by the strong exchange interaction, J, among localized spins, and besides, we 
are interested in a semi-macroscopic object made of many spins, the domain 
wall. The localized spin interacts with the conduction electron by an s-d type 
exchange interaction (Fig. 3), 

H sd = -J sd Jd 3 xS-(Jcrc). (2) 

Here the conduction electron is represented by creation and annihilation oper- 
ators c^ and c, and er = (a x , <J y ,a z ) where cr, are 2x2 Pauli matrices satisfying 
the commutation relation [ctj, cxj] = 2i€ijk<Jk- The description based on this s-d 
exchange picture is an effective one, treating localized spin S as a variable 
independent from the conduction electrons, i.e., neglecting the hopping of d 
electrons that form localized spin. Still, we will take this effective s-d model 
as the starting system for this investigation, and will not concern ourselves 
with the microscopic origin of the local moment. 




Fig. 3. The s-d interaction induces precession of localized spin and electron spin 
around each other. 

The s-d interaction is a coupling in spin space, which is decoupled from real 
space (as far as spin-orbit interaction is neglected). Nevertheless, this spin 
coupling can affect charge transport if the localized spin has inhomogeneity, 
and various magnetoresistive effects such as GMR arise. 

Since this exchange coupling describes the exchange of spin angular mo- 
mentum, the idea of spin reversal by spin-polarized current arises naturally. 
Namely, the injection of electron spin polarized in the opposite direction to a 
localized spin will cause flip of localized spin (Fig. 4). This simple idea was 
integrated into realistic magnetization switching of thin film magnets by Slon- 
czewski [8] and Berger [9]. The current- induced phenomena are expected to 
be applied to memory devices like magnetoresistive random access memory 
(MRAM) that operates without magnetic field, and intensive studies on pillar 
systems and domain walls have then started at the end of the last century. 
Domain-wall racetrack memory proposed by Parkin is one possibility of the 
high-density storage [10,11]. 

Compared with switching by use of the Ampere's field, the currend-induced 




Fig. 4. s-d exchange interaction describes the exchange of spin angular momentum 
between conduction electron and localized spin. The idea of spin flip using spin-po- 
larized current is thus quite natural. 

magnetization switching has a great advantage in downsizing. The field created 
by the Ampere's law is proportional to the current, which decreases when 
the system size is reduced with a constant current density. Therefore, higher 
current density is necessary for the Ampere's mechanism in smaller systems. 
In contrast, the current-induced switching rate is determined by the current 
density and material parameters, such as s-d coupling, spin relaxation, and 
anisotropy energies, and the efficiency remains constant when the system size 
is reduced. This is why current- induced switching becomes essential in high 
density devices. 

In this paper, we review recent developments in the theory of current-driven 
domain wall motion. In §2 and §3, the phenomeno logical argument and a brief 
history of current- induced domain wall dynamics are presented. The theoreti- 
cal study starts in §4 from the description of localized spin by the Lagrangian 
formalism. Collective coordinates to describe domain wall are introduced in 
§5. We consider the case of a rigid one-dimensional (planar) domain wall, 
sometimes called the transverse wall. The conduction electrons and s-d inter- 
action are introduced in §6. The equation of motion of a domain wall coupled 
to the conduction electrons is derived in §7. The equation is expressed using 
the conduction electron spin density, which acts as the effective field on the 
localized spins. Then the explicit equation of motion is obtained by calculating 
the conduction electron spin density in §8. The torque and force acting on the 
spin structure are obtained in §9. The adiabatic limit is briefly discussed in 
§10. The full equation of motion of a domain wall is finally obtained in §11 
and is solved in §12. The case of a wall having vorticity, called the vortex wall, 
is considered briefly in §13. The analysis in §4 to §12 is the main result of the 
paper, aiming at presenting our calculational method in a self-contained way. 

Another approach to current- induced domain wall dynamics is to use the 
Landau-Lifshitz-Gilbert (LLG) equation taking account of the effect of cur- 
rent (as done in §7.3). To do this, we need to calculate microscopically the 
torque induced by the electron. This is done in §14. In §15, electron trans- 
port properties in the presence of spin structures are discussed. The transport 
properties are shown to be the counter action of the current-induced forces. 
Another counter action of current-induced magnetization dynamics, the pump- 
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ing of current and spin current by magnetization dynamics, is briefly argued 
in §16. Details of the Lagrangian formalism of spin and Green's functions are 
explained in the appendices. 



2 Current-driven domain wall dynamics 

2.1 Switching using domain wall motion 



The idea of switching spin structure by electric current using the s-d inter- 
action was first discussed by Berger in 1978 [12], much earlier than works by 
Slonczewski and Berger in 1996. His idea was to push a domain wall by cur- 
rent. A domain wall is a twisted spin structure where spins gradually rotate, 
which appears between two magnetic domains [13-18] (Fig. 5). The thickness 
of the wall, A, is determined by the competition between the ferromagnetic 
exchange energy (between localized spins), J, which aligns the neighboring 
spins, and the magnetic anisotropy energy in the easy axis, K , which tends 
to reduce the wall thickness to minimize the deviation of spins from the easy 



axis, as A = Jj/K. (Here J has dimensions of J/m 2 .) Thus A depends on the 
material and also on the sample shape since K depends on the shape. In the 
case of 3d transition metals such as iron and nickel A ~ 500 ~ 1000 A[19,20], 
and A ~ 150 A in Co thin films [21]. This length scale is very large compared 
with the length scale of the electron, kp 1 ~ 0(1 A) (kp being the Fermi wave 
length of the electron) . 



x ttmvm^v 
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Fig. 5. Illustration of a Neel wall (left) and a Bloch wall (right), where the magnetic 
easy axis is along and perpendicular to the wall direction, respectively. A is the 
thickness of the wall. 

Let us consider how the motion of the wall is induced by electric current. 
When an electron is injected into the domain wall, there are basically two 
possibilities, reflection or transmission. In the transmission process, there are 
again two possibilities as indicated in Fig. 6, depending on the electron speed. 
If the electron is fast enough, it will pass through the wall without spin rota- 
tion, while the electron spin will be flipped by exchange coupling during the 
transmission if electron is slow. 

Corresponding to the above possibilities, there are two different mechanisms 
of domain wall motion induced by electric current and exchange interaction. 
The first one is due to reflection of the electron. The exchange interaction 
describes a spin-dependent potential created by a localized spin S, and so the 
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Fig. 6. Schematic figure describing a conduction electron going through a domain 
wall. The electron spin is denoted by small arrows and localized spins are denoted 
by large arrows. The upper figure is an adiabatic case, where the electron is slow 
and so can accommodate its spin as electron passes through the wall, (v is velocity 
of the electron.) The lower figure is the case of a fast electron, where the electron 
spin cannot follow the domain wall structure. 

electron is scattered if there is inhomogeneity, VS. Namely, the electron feels 
a force from domain wall, 

F ex Rj, (3) 

where R is the reflection probability for the electron. (For correct expression, 
see Eq. (136) and Eq. (235).) From the conservation of linear momentum, 
electron scattering by a domain wall indicates that the wall must move. This 
process is an exchange of linear momentum, and is sometimes called a momen- 
tum transfer process. This force is strong if the domain wall is thin, since then 
the electron scattering is significant. In reality, in most experiments, domain 
walls are thick, and the exchange interaction is strong, and so most of the 
electrons do not get scattered (except for systems with very thin walls [22] 
or in nano scale contacts [23]). This case is called the adiabatic case, and is 
suggested in experiments by small resistivity due to domain walls [17]. (For 
conditions of adiabaticity, see §6.5.) Thus this force is not a major driving 
mechanism in most cases as discussed by Berger [24]. 

The other mechanism arises from the adiabatic electron transmission. As seen 
in Fig. 6, the spins of slow electrons are nipped on transmission. The angular 
momentum of a conduction electron has changed by the amount h x (| — 
(— |)) = h when one electron goes through. From the conservation of angular 
momentum, the wall needs to shift by a distance of AX = ^ (a is the lattice 
constant) (Fig. 7). When a steady current density j is injected, the wall then 
moves at speed of 

. ^ r Pja 2 Pja 3 .„. 

e 2eb> 

where P = 3+ ~ 3 r is the spin polarization of current (j± represents the current 
carried by the electron with spin ± and j = j + + j-). This is so called spin- 
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transfer mechanism of domain wall motion. This argument applies to any 
adiabatic spin structure, and we can see that any spin structure tends to flow 
at the speed given by Eq. (4). The direction of wall motion is the same as that 
of the electron, and so is opposite to the current (since the electron charge e 
is negative). 
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Fig. 7. Domain wall motion by injection of spin-polarized electron in the adiabatic 
limit (spin transfer effect). 

In reality, two driving mechanisms exist and so the two equations, Eq. (3) 
and Eq. (4), need to be coupled. One may naively guess simply that M W (X + 
a'X) = F and X = ^|- (where M w is the wall mass and a' is the friction 
coefficient), but these are not correct, since the wall is not a simply a particle 
but has internal degrees of freedom. (One may notice also that these two 
equations for constant current have no solution.) Berger has proposed, based 
on phenomenological arguments, the correct equation in his series of papers 
[12,24]. Actually, a force on a domain wall induces not a simple acceleration 
(X) but an angle out-of the easy plane, (see Figs. 8 and 9). This has been 
known in the case where a magnetic field B is applied along the easy axis. 
Let us visualize the motion in this case. Under a magnetic field, each spin 
constituting the domain wall starts to precess around the field according to a 
torque equation of motion, 

S = ~fSxB, (5) 

where 7 (> 0) is the gyromagnetic ratio. (We will use magnetic flux density 
B instead of magnetic field H {B = (J,(H + M)), and B may be called the 
"magnetic field", as in Kittel's textbook [25].) The spin thus changes its direc- 
tion perpendicular to S and B. The translational motion of the domain wall 
is therefore coupled with the out-of plane dynamics, and this is an essential 
and complicated feature of the wall dynamics. The correct equation under a 
force is given by [15,26,27] (when friction is neglected) 

= VF, (6) 

where 77 is a numerical factor. (It turns out that 77 = ^ „ , where N w is number 
of spins inside the wall.) The effect of current in the adiabatic limit is to induce 
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a wall velocity as we saw, but the wall velocity is also related to the hard-axis 
anisotropy energy (if it exists), since the translational motion of the wall needs 
to be induced by the effective magnetic field perpendicular to the wall plane, 
again due to the precession equation (5). The other equation for the wall is 
therefore written as (without friction) 



X 



Pja 3 



+ v c sin 20, 



(7) 



where v c is a parameter that determines the hard axis anisotropy energy (Eq. 
(270)). (Here, assuming hard axis anisotropy of the standard sin 2 type, the 
effective field perpendicular to the wall plane is given as ex sin 20.) These two 
equations Eqs. (6) (7) are not yet correct, since they do not include the effect 
of damping (friction), which is known to be quite essential in spin dynamics 
[15]. Damping can be incorporated phenomenologically by the prescription by 
Gilbert or derived from spin relaxation processes [28,29] (see Eqs. (92) (93)). 




Fig. 8. Domain wall configuration where the easy axis is along the wall direction, 
called the Neel wall. The angle out of the easy plane is 0, which is a canonical 
momentum of the wall and which plays an essential role in dynamics. The dynamics 
of (called the chirality of the wall) is discussed in Refs. [30,31]. 




Fig. 9. Domain wall configuration where the easy axis is perpendicular to the wall 
direction, called the Bloch wall. The easy plane is perpendicular to the wall direction 
as denote by a small square. The angle out of the plane, 0, corresponds to the tilt 
of wall in the wall direction. 

Let us see how the wall motion changes when the sign of parameters changes. 
If J s d is negative, the spin-transfer torque gives a wall velocity opposite to the 
case of J s d > 0, resulting in wall motion in the current direction (if carrier has 
negative charge). Mathematically, this is because M, P, and j s change sign 
with J s d. The forces due to non-adiabaticity and spin relaxation (_F and F&) 
remain opposite to the current. In the case of a hole in semi-conductors, the 
charge of the carrier is positive, and so the reflection force is along the current 
direction. The spin-transfer torque pushes the wall in the same direction as in 
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the electron case if J s d has the same sign. For instance, in GaMnAs, exchange 
interaction between a p hole and the d localized spin is negative ( J p ^ < 0) [32] 
and so spin-transfer motion is opposite to the hole flow and the current. 

As is seen from the above arguments on spin transfer torque, what matters 
most is spin polarization of current, 

P= a+ - a - = n+r+ ~ n - r - , (8) 

a + + ct_ n + r + + n^T- 

where o~±, n±, and r± are the Boltzmann conductivity, density, and lifetime of 
electrons with spin ±, respectively. This polarization P should not be confused 
with other definitions of spin polarization, such as polarization of electron den- 
sity, n +~ n - or that of density of states, y +~ v - = , F+ 7,' F ~ (in three-dimensional 
free electron model). It would be interesting to control the sign of P in exper- 
iments and see how the wall motion changes. 

When friction is included, these two equations work fine for qualitative argu- 
ment close to the adiabatic limit. However, these phenomenological equations 
are too naive for quantitative study. For instance, the definition of 4> in Eq. 
(6) is not clear for a domain wall whose <fi can be position dependent, and 
the explicit form of F is not given. In addition, to include the effects of spin 
relaxation and non-adiabaticity correctly, one needs to use a many-body for- 
malism treating the electron quantum mechanically. This is what we are going 
to do in this paper. Based on a Lagrangian formalism, the equation of motion 
is derived without worrying about complicated dynamics of each spin. (The 
result is Eq. (268).) Mathematically, the coupling between the translational 
and out-of-plane dynamics is expressed by the fact that the canonical momen- 
tum of the wall is the average of (and not M W X), a fact arising from SU(2) 
commutation relation of spins. 



3 Brief history 

3.1 Berger's theories 



Berger considered a domain wall under an electric current, and saw that the 
s-d exchange coupling between the localized spin and conduction electron spin 
is the dominant interaction that drives the wall under a current in the case of 
a thin film (e.g., thickness less than ~ 0.1 /zm), where the effect of an induced 
magnetic field can be neglected [12]. In 1984 [33,34], he studied the effect of 
the force arising from the reflection of conduction electrons by the domain 



15 



wall caused by this exchange coupling. This force was obtained as 

F = M s a f {j-b f v w ), (9) 

where M s is the saturation magnetization, a/ and bf are coefficients introduced 
phenomenologically, j is the current, and v w is the wall velocity. The effect of 
the force was found to be small in most cases due to a very small reflection 
probability because the wall thickness is usually large compared with the Fermi 
wavelength. In 1978 [12,35], he argued that the exchange interaction produces 
a torque, 

h , , 

T = —j, (10) 

e 
(assuming full spin polarization of the conduction electron), which tends to 
cant the wall magnetization out of the easy plane (angle 0) and eventually 
induces a continuous rotation of a pinned wall under a large current [35]. 
This torque was found to push the wall by a different mechanism from the 
exchange force, which turned out to be the dominant driving mechanism [24]. 
The torque is nowadays called the spin transfer torque, after Slonczewski [8]. 
Based on the idea of torque-driven wall motion, an experimental study was 
carried out in 1993 [36] on a thin film of NigiFeig. There, a domain wall velocity 
of 70 m/s was reported at the current density of 1.35 x 10 10 A/m 2 applied as 
a pulse of duration 0.14 fjs. 

There has been a renewal of interest on the current-induced domain wall mo- 
tion for about a decade. Recent experimental studies have been carried out 
on submicron-size wires, and the domain wall motion induced by current has 
been confirmed [37-40]. The current density necessary for wall motion turned 
out to be rather high, of order of 10 12 A/m 2 . Measurement of the domain wall 
velocity was carried out by Yamaguchi et al. [41] by observing wall displace- 
ment by use of magnetic force microscopy (MFM) after each current pulse of 
strength 1.2 x 10 12 A/m 2 and duration of 5 /is. The average velocity was found 
to be 2 ~ 6 m/s. Other experiments also indicate a rather slow average wall 
velocity, of order of a few m/s under a steady current [42] . 



3.2 Recent theories 



Those experiments motivated theoretical studies to look into the problem in 
more detail. Microscopic derivation of the equation of motion of the domain 
wall under current was carried out by Tatara and Kohno [26,27,43]. They 
considered a planar (one-dimensional) wall and described the wall by the two 
collective coordinates, X and <po, i.e. within Slonczewski's description [44]. The 
variable X represents the position of the wall, and 0o describes the tilt of the 
wall plane. Considering a small hard-axis anisotropy case, other deformation 
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modes than 0o (such as change of wall width) were neglected (rigid wall ap- 
proximation). The equation of motion with respect to X and 0o was derived 
including the effect of conduction electrons via the s-d exchange interaction. 
The electron carrying a current was treated by the use of a non-equilibrium 
(Keldysh) Green's function. Qualitatively, the equation of motion derived was 
indeed the same as that obtained by Berger long ago [33,24], namely, Eqs. 
(6) (7) with the damping term included. Berger's theory was thus confirmed 
by microscopic calculation. The microscopic formalism made it possible for 
the first time to calculate the spin transfer torque and force systematically 
by representing these quantities by Green's functions and Feynman diagrams. 
Based on the obtained equation of motion, the wall motion under steady cur- 
rent was studied. It was found that in the adiabatic limit, where the reflection 
force can be neglected, and in the absence of spin relaxation of conduction 
electrons (/3 sr term below), there is an threshold current determined by the 
hard-axis magnetic anisotropy energy, K± as 

*■ = ■&&*>■*■ (11) 

This is the intrinsic pinning of the wall arising from the "pinning" of <$$ [26]. At 
larger currents, the wall gets depinned and its velocity becomes proportional 
to the spin current (spin polarization of the current flow), j s , as is required 
from the angular momentum conservation. 

Of practical importance would be the wall motion below the intrinsic thresh- 
old. Actually, the wall moves over quite a large distance (even in the absence 
of (3 term) if there is no pinning. According to the torque equation, when 
j <C jl the current induces a tilt of the wall, (fio ~ A- (see Eq. (272)). This 
tilt is associated with wall translation (by the second equation of Eq. (272) 
without pinning and (3 W ) over a distance of 

AX--H (j«j'i). (12) 

Since a is very small, e.g., a ~ 0.01, this distance can be quite large compared 
with A even for a current 10% of the intrinsic threshold. Such motion at 
very low current would be enough for applications. One should note, however, 
that the wall below threshold goes back to the original position as soon as 
the current is cut. One needs therefore a pinning site to maintain the wall 
displacement. 

Numerical simulation was performed based on an equation of motion of each 
localized spin by including the spin-transfer torque term in the adiabatic limit 
[45] . The equation of motion is given by 

S = 1 B s xS--SxS- —(j B • V)5. (13) 
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Here Bs is the effective field arising from the spin Hamiltonian, and a repre- 
sents damping. The equation S = 7-B5 x S — ^S x S has been well-known 
as the Landau-Lifshitz-Gilbert equation describing magnetization dynamics 
in a magnetic field Bs- Spin-transfer torque from current is represented by 
the last term of Eq. (13). The simulation result was similar to the analytical 
(collective-coordinate) study, indicating the existence of an intrinsic threshold 
current. The motion of the domain wall under magnetic field and spin-transfer 
torque was solved in Ref. [46]. 



3.2.1 The Landau-Lifshitz-Gilbert equation under current 

Later Zhang and Li [47] and Thiaville et al. [48] proposed to add a new torque 
term in the equation, which is perpendicular to the spin-transfer torque. After 
Thiaville et al. [48], we call this torque term the (3 term. Zhang and Li ar- 
gued that the (3 term arises from spin relaxation of conduction electrons [47]. 
Thus the phenomenological equation of motion of localized spin under current 
becomes 

S = 1 B s xS--SxS-—(j s -V)S- 1 ^lSx (j s • V)S] + r na . (14) 

The fourth term is the new (3 term. We explicitly wrote the coefficient (3 sr with 
a suffix sr to show that this term arises from spin relaxation. The last term, 
r na , denotes the non-adiabatic torque, which is spatially nonlocal [49,50]. (So 
far, this nonlocal torque has not been taken account in numerical simulations.) 
All these torques are derived in this paper based on Eq. (130). Parameters a 
and (3 sr are calculated in §14, and the nonlocal torque r na is given by the 
nonlocal part of Eq. (223). 




V)S 



Fig. 10. Torques acting on localized spin S: damping torque, S x S, spin transfer 
torque, (j ■ V)S, and (3 torque, S X (j ■ V)S. 

The j3 term turned out to modify the threshold current and the terminal veloc- 



ity of the wall significantly if it is not small compared with damping parameter 
a [47,48,51]. For a rigid domain wall, this (3 sr in the Landau-Lifshitz-Gilbert 
equation has exactly the same effect as the force due to non-adiabaticity (elec- 
tron reflection by wall) [48,51]. In other words, the nonlocal torque r na for a 
rigid wall is effectively represented by an additional /3 term, /3 na [51]. Therefore, 
the effective f3 or force that the rigid wall feels is given as 

£w=&r + /U (15) 

where /3 na denotes the contribution from non-adiabaticity (see Eq. (269) and 
Eq. (251)). (For other spin structures such as vortices, an expression similar to 
Eq. (15) holds.) The Gilbert damping parameter originates from electron spin 
relaxation and other sources, so we can write a = a sr + ao (see Eq. (352)). 

It has been shown that when j3 w /a > 1, the intrinsic threshold is smeared 
out and the true threshold current is determined by extrinsic pinning [48,51]. 
The terminal wall velocity is also determined by /3 w /a [47,48,51]. The param- 
eters j3 sr and /3 na are sensitive to the spin relaxation rate and wall thickness, 
respectively. It would be therefore very interesting to experimentally identify 
the origin of f3 w by changing material (spin-orbit interaction) and structure 
(wall thickness). 

Microscopic derivation of the /3-term has been carried out by several authors 
[28,29,52,53]. Tserkovnyak et al. [28] calculated /3 sr based on a one-band model 
considering spin-relaxation of conduction electrons semi-classically and as- 
suming spin dynamics of small amplitude. They considered the limit of weak 
ferromagnetism and found that (3 sr = a sr . They also showed that in general 
(3 ST 7^ a ST considering multiband effects or deviation from weak ferromag- 
netism. Their approach is, however, still phenomenological, treating the spin- 
flip process by a phenomenological spin-relaxation time in the equation of mo- 
tion of spin. The relation (3 ST = a ST was suggested also by another phenomeno- 
logical argument [54,55] (but see also Ref. [56]). Fully microscopic calculation 
of j3 and a due to spin relaxation was carried out on an s-d model by Kohno et 
al. [29,52] using standard diagrammatic perturbation theory, where the effect 
of spin relaxation are taken into account consistently and fully quantum me- 
chanically. The result indicated /3 sr ^ a sr . The same result was obtained later 
in the functional Keldysh formalism by Duine et al. [53] . Determination of (3 
and a values needs a careful microscopic calculation, since they are quantities 
smaller by a factor of l/(e^r) compared with conventional transport coeffi- 
cients. A phenomenological thermodynamic argument predicted j3 = a [54,55], 
but microscopic studies [29,53,52] indicate that that is invalid. The error may 
arise because the argument of Refs. [54,55] lacks consistent consideration of 
the work done by the electric current [56]. So far, the effect of spin relaxation 
from spin-orbit interaction has not been calculated. The effect is expected to 
be essentially the same as spin-flip impurities [52] . Such a term like (3 was also 
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theoretically found in ferromagnetic junctions [57]. 

Experimentally /3 W can be determined from domain wall dynamics [58,59]. A 
very large ratio of - ~ 8 was obtained from oscillating wall motion in a pinning 
potential [58], and the significant deformation of the domain wall observed by 
Heyne et al. [59] clearly indicated j3 w ^ a. Recent systematic experimental 
study [60] on the spin damping parameter in various ferromagnetic films re- 
vealed that a scales with (g — 2) 2 (g is the (/-factor), indicating that damping 
arises mainly from spin-orbit interaction [61]. 

Sometimes the f3 torque from spin relaxation is called the non-adiabatic torque 
[47] , extending the definition of non-adiabaticity to include spin deviation due 
to spin relaxation. The original meaning of adiabaticity in electron transport 
is the absence of scattering or the exchange of linear momentum. In this paper, 
we use the term non-adiabatic in this original, narrow sense. We therefore call 
non-adiabatic torques only those arising from momentum transfer (finite q 
contributions in Eqs. (215)(216)(217)(219) of torque given by Eq. (226)). Ac- 
cording to this definition, the (3 sr term from spin relaxation is still an adiabatic 
torque, since /3 sr is defined as the local torque in the Landau-Lifshitz-Gilbert 
equation (Eq. (14)). Correctly speaking, spin relaxation has additional truly 
non-adiabatic components that are nonlocal and oscillating (contributing to 
r na in Eq. (14)), which have not been discussed so far. 

The Landau-Lifshitz-Gilbert equation is useful in numerical simulation of spin 
dynamics. However, one should pay attention to the fact that it is a phe- 
nomenological equation utilized only within a local approximation of torques, 
and is not suitable for studies including non-adiabatic effects (i.e., non-local 
torques). The analysis in this paper is therefore not based on the Landau- 
Lifshitz-Gilbert equation. Instead, we will derive the equation of motion of 
the wall directly from the total Hamiltonian. 

Spin damping torque in the magnetic field can also be expressed using the 
Landau-Lifshitz prescription. The equation of motion then reads 

S = lLL B s x 5 + | 7 ll5 x (5 x B s ), (16) 

where 7ll is a constant. This equation is equivallent to the Landau-Lifshitz- 
Gilbert (LLG) equation 

S = jB s xS-^SxS. (17) 

In fact, Eq. (16) indicates 7ll(5 x Bs) = —5 + f 7ll5 x (5 x Bs), and thus 
the equation reads 

5 = 7 ll (l + a 2 )B s xS- -(5 x 5). (18) 
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Therefore the two equations are equivalent either by defining 7 in the LLG 
equation as 7 = 7ll (1 + « 2 ), or by neglecting 0(a 2 ). (Note that both equa- 
tions are approximations neglecting higher derivatives such as S.) From the 
microscopic theory point of view, the Landau-Lifshitz-Gilbert equation seems 
natural, since, as we will see in §14, the damping torque of the Gilbert type, 
|5 x S, is derived by a systematic gauge-field expansion. 



3.2.2 Non-adiabatic torque and spin-orbit interaction 

Waintal and Viret [62] and Xiao, Zangwill, and Stiles [49] studied the spatial 
distribution of the current-induced torque around a domain wall by solving 
the Schrodinger equation and found a nonlocal oscillatory torque (r na in Eq. 
(14)). This torque is due to the non-adiabaticity arising from the finite do- 
main wall width, or in other words, from the fast-varying component of spin 
structure. The oscillation period is ~ kp~ (kp is the Fermi wavelength) and is 
of quantum origin similar to the Ruderman-Kittel-Kasuya-Yoshida (RKKY) 
oscillation. This nonlocal torque looks complicated, but represents in fact a 
force acting on the wall, studied in Ref. [26]. In this paper, we demonstrate 
that this oscillating torque on each spin (Eq. (227)) is indeed summed up to 
a force (Eq. (233)) when looked collectively [50]. 

Ohe and Kramer [63] studied wall motion solving the torque due to the ex- 
change interaction numerically, including non-adiabaticity. Non-adiabaticity 
was studied further in Refs. [50,64,65]. Thorwart and Eggar [65] applied a 
one-dimensional method to evaluate the conduction electrons, and carried out 
a gradient expansion with respect to slowly varying localized spins. The torque 
in the second order of the spatial derivative was derived and was shown to de- 
form the wall significantly if the wall is very narrow, A ~ 5.5a. 

Nonlocal oscillating torques were numerically studied by taking account of 
the strong spin-orbit interaction based on the Luttinger Hamiltonian (i.e., in 
magnetic semiconductors) [66]. It was shown there that due to the spin-orbit 
interaction, the oscillating torque becomes asymmetric around domain wall 
and that this feature results in high wall velocity. Current- induced domain wall 
motion in the presence of Rashba spin-orbit interaction was recently studied 
[67]. It was found that for one type of Bloch wall, a very large effective (3 term 
is induced by the Rashba interaction, which acts as an effective magnetic field, 
and wall mobility is strongly enhanced. A strong Rashba interaction is known 
to arise not only in semiconductors, but on the surface of metals doped with 
heavy ions [68,69], and such systems would be very interesting in the context 
of domain wall dynamics. 
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3.3 Recent experiments 



3.3.1 Metallic systems 

So far experimental results on metallic samples all show threshold currents of 
the order of 10 12 A/m 2 . If we use K±/ks ~ 0(1 K) estimated experimentally, 
[41,70,71] the observed threshold is orders of magnitude (10~ 2 — 10 _1 times) 
smaller than the intrinsic threshold, j l c . For instance, a sample of Yamaguchi 
[41,70,71] showed j c — 1 x 10 12 A/m 2 . The anisotropy energy is estimated 
to be K±/k B = 2.4 K, and using S ~ ~, a ~ 2.2 A and P ~ 0(1), we 
obtain j\ = 5.8 x 10 13 A/m 2 , i.e., jc/Jl ~ 0.02. The observed low threshold 
currents in metals thus should be regarded as due to an extrinsic pinning. 
Actually, direct evidence excluding intrinsic pinning in permalloy wires so 
far was given by Yamaguchi et al. [72]. They prepared permalloy wires with 
different geometries, and realized different perpendicular anisotropy energies 
^K ± ~ (0.1 - 7.6) x 10 5 J/m 3 , which corresponds to K±/k B ^ 0.03 - 2.4 K 
(per 1 spin). The intrinsic pinning, Eq. (11), predicts then a threshold current 
of 5 x 10 11 — 4 x 10 13 A/m 2 . In contrast to such a difference in the predicted 
values of the intrinsic threshold current, experimental values of the threshold 
for these samples do not vary so much, (3 — 8) x 10 11 A/m 2 , and are smaller 
than the predicted intrinsic threshold by factors of 2 to 100. Besides, data by 
Yamaguchi et al. indicate that these experimental values do not scale with K±, 
although there is a weak dependence on K±. Therefore the observed threshold 
in Ref. [41,71] would be of some extrinsic origin. 

The wall speed is another important quantity to determine the driving mech- 
anism and efficiency Under long (> fis) current pulses in metallic samples, 
the wall speed so far is very small, less than 1 m/s [42] or 2 ~ 6 m/s at a 
current density of 10 12 A/m [71]. This is far below the perfect spin transfer 
limit, v w = ^gPj — 100 m/s at j = 10 12 A/m . Processes involving strong 
dissipation of angular momentum or linear momentum during wall motion, 
such as deformation by extrinsic pinning centers or spin wave emission, might 
explain the discrepancy. Measurements on clean samples are necessary. (In 
semiconductors, in contrast, perfect spin transfer seems to be realized. See 
below.) 

Direct observation of the spin structure indicates that the wall is considerably 
deformed upon motion [42,73-75]. It was shown [42] that the initial state is 
not a planar wall but more like a vortex (called a vortex wall), which is the 
case in film or wide wires, the vortex wall moves by applying a current pulse 
of 2.2 x 10 12 A/m 2 , and the wall is deformed to become a transverse wall 
after some pulses. It is interesting that although the vortex wall moves more 
easily, the transverse wall does not move at the same current density. This 
would be explained by the absence of intrinsic pinning for the vortex wall 
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(Eq. (320)). As far as the transverse wall is concerned, rigidity or non-rigidity 
does not seem essential since the analytical results for a rigid wall [26,50] 
and numerical simulation including deformation effect [48,76] predict similar 
behaviors. 

Heating effect in metallic samples has indeed been found to be crucially impor- 
tant [70,77]. For applications, heating may assist wall motion. Sub ns pulses 
were reported to be quite efficient in driving the wall at a low current den- 
sity of ~ 10 10 A/m 2 [78]. This motion could be the motion due to j3 term or 
below the intrinsic threshold. Short pulses could be useful for applications. 
Quite recently, Laufenberg et al. measured the temperature dependence of 
the threshold current and found that it decreases at low temperatures, for 
instance, from 2.4 x 10 12 A/m 2 at T ~ 170 K to 1.9 x 10 12 A/m 2 at T ~ 100 
K [79]. Dissipation of spin-transfer torque by spin waves was suggested as a 
possible explanation, but a theoretical study is yet to be done. 

For recent experimental results, see Ref. [18]. 



3.3.2 Thin wall 

Quite an interesting result was obtained recently by Feigenson et al. [22] 
in SrRuOs, an itinerant ferromagnet with perovskite structure. The current 
density needed to drive the wall was 5.3 x 10 9 A/m at T = 140 K and 
5.8 x 10 10 A/m 2 at T = 40 K. A small threshold current at 140 K would be 
due to reduction of magnetization close to T c = 150 K. The threshold current 
is about 2 orders of magnitude smaller than in other metals. This high effi- 
ciency would be due to a very narrow domain wall, A ~ 3 nm, as a result of 
very strong uniaxial anisotropy energy (K) corresponding to a field of 10 T. 
They defined a parameter determining the efficiency as a ratio of the depinning 
field and depinning current density, A = B c /j c . Their results were A = 10 -12 
T m 2 /A. They compared this value with the threshold current of extrinsic pin- 
ning [26] (given by Eqs. (308) and (310)). Using -f^ ~ 0.5 x 10~ n T m 2 /A 
and S ~ 3/2, we see that A = 10~ 12 T m 2 /A is realized if (3 W ~ 0.5. This value 
would be too large if interpreted as due to spin relaxation. Using the measured 
resistivity of the domain wall, the non-adiabatic force contribution to /3 W was 
estimated and the result of j c was of similar order as the observed ones but 
with a discrepancy of a factor of around 6 at low temperature (Fig. 4(a) of Ref. 
[22]). This discrepancy seems not very crucial considering the crude rigid and 
planar approximation of the wall. There is another extrinsic pinning threshold, 
j* a ) (Eq. (307)). If we use this expression, A = 10 -12 T m 2 /A corresponds to 
K± ~ ^bB c = 4 x (10~ 3 — 10 -2 ) K (per site). This number seems also possible 
although K± was not measured. 

Controlling magnetic anisotropy in metallic thin films is an interesting pos- 
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sibility. An easy axis perpendicular to the film is realized using a Pt layer, 
and in such case, the domain wall is a Bloch type and the thickness is known 
to be small, A ~ 12 nm [80]. Ravelosona et al. [80] found in such perpen- 
dicular systems that the threshold current density is reduced by a factor of 
~ Yq compared with other metallic systems, j c ~ 10 11 A/m . The reduction 
could be explained by the intrinsic threshold, Eq. (11), as due to reduction 
of A by roughly a factor of -^ if K± is the same order as in-plane anisotropy 
systems. It could be possible also that K± is very small in perpendicular sys- 
tems, since the system is quite symmetric within the plane perpendicular to 
the easy axis, resulting in reduced threshold current. Although non-adiabatic 
effects are mentioned in Ref. [80], a 12 nm wall would be still in the adiabatic 
regime, since kp\ ~ 0(100). Reduction of threshold current in perpendicular 
anisotropy films was also confirmed in numerical simulation [81,82]. 



3.3.3 Magnetic Semiconductors 

Beautiful experiments were carried out at low current in ferromagnetic semi- 
conductors by Yamanouchi et al [83,84]. They fabricated a wall structure of 
20 /zm width made of GaMnAs with different thicknesses, which determines 
the ferromagnetic coupling and transition temperature, and trapped a domain 
wall. The wall position was measured optically after applying a current pulse, 
and the average velocity was estimated. The current necessary was ~4x 10 9 
A/m 2 , which is 2-3 orders of magnitude smaller than in metallic systems. 
This is due to the small average magnetization, S ~ 0.01, carried by dilute 
Mn ions, and small hard-axis anisotropy K± [84]. The obtained velocity was 
rather high, ~ 22 m/s at j = 1.2 x 10 10 A/m 2 . This velocity is consistent with 
the adiabatic spin-transfer mechanism, Eq. (4), and the threshold appears to 
be consistent with the intrinsic pinning mechanism [26] with anisotropy energy 
obtained from band calculations. 

However, there are some puzzles. First, the theory of intrinsic pinning [26] and 
adiabatic spin transfer does not take account of strong spin-orbit interaction in 
semiconductors. So the agreement with these theories might be a coincidence. 
The second puzzle is the validity of using purely adiabatic theory. In fact, 
quite a large momentum transfer (force) is expected from the wall resistance, 
R w = 1 Q [85], corresponding to j3 w ^> 1 in terms of (3 W [85]. 

Another puzzle, which was solved just recently, is the temperature dependence 
of the wall velocity. The observed velocity scaled as In?; ~ — (T c — T) 2 ^" 1 / 2 , 
similar to the creep behavior under a magnetic field [86], but this fractional 
power of j has not been explained in the current-driven case. A simple theory 
of thermal activation assuming a rigid wall under the spin-transfer torque 
predicts a different behavior, lnt> ~ j/T [87], and thus creep motion would 
be essential in the experiment by Yamanouchi et al. Successful explanation of 
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creep behavior was just recently offered by Yamanouchi et al. [88], by taking 
account of growth of at the pinning center. The thermal effect and creep by 
current was studied in more detail by Duine et al. [89]. 

Nguyen et al. studied theoretically the domain wall speed in magnetic semicon- 
ductors based on the 4-band Kohn-Luttinger Hamiltonian [66]. It was shown 
there that the wall speed can be enhanced by the spin-orbit interaction by a 
factor of 10 3 — 10 4 due to the increase of mistracking, hence reflection, of con- 
duction electrons. Rashba spin-orbit interaction is also expected to be useful 
for efficient magnetization switching [67]. 



3.3.4 Excitation of wall 

Time-resolved study of excitation of wall provides rich information on the wall 
character and driving mechanism. 

Under alternating current (AC), the domain wall shows another aspect not 
seen in the direct current (DC) case. AC can drive domain walls quite ef- 
fectively at low current if the frequency is tuned close to the resonance with 
the pinning frequency. This resonance was realized in a recent experiment by 
Saitoh et al. [20]. They applied a small AC (of amplitude of 10 10 A/m 2 ) in 
a wire with a domain wall in a weak pinning potential controlled by a mag- 
netic field. Although the current is well below the threshold, the wall can 
shift slightly as we see below (distance is calculated to be around /im, which 
might be an overestimate). Under a small current, 4>q remains small, and the 
equation of motion reduces to that of a "particle" ; 

M w X + ^X + M w n 2 pin X = F(t), (19) 

where M w is the wall mass, r w oc a -1 is a damping time, f2 pin is the (extrinsic) 
pinning frequency, and F(t) is a force due to current. For AC, I(t) = Ioe tuJt , 



i(t) 



where to is the frequency, the force is given F(t) = -^ ^Pw — iPf\ ' jf-j 
where j3 w , given by Eq. (269), is the total force from momentum transfer and 
spin relaxation, and the last term, proportional to o>, is from the spin-transfer 
torque. The wall under a weak current thus shows a forced oscillation of a 
particle. By measuring the energy dissipation (from complex resistance), a 
resonance peak would then appear when oj is tuned closely around fi pin . From 
the resonance spectra, the mass and the friction constant were obtained as 
M w = 6.6 x 1CT 23 kg, r w = 1.4 x 1CT 8 sec. The experimental result seems to 
be well described by the rigid-wall picture, and this would be due to a low 
current density (by factor of 10~ 2 compared to DC experiments on metals), 
resulting in a small deformation. What is more, from the resonance line shape, 
the driving mechanism of the domain wall was identified to be the force (/3 W ) 
rather than the spin-transfer torque. This finding was surprising at that time, 
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when the adiabatic spin-transfer torque was considered as the main driving 
mechanism. The observed force corresponds to the value of /3 W ~ 1.5, which 
is too large if j3 w arises from spin relaxation j3 sr (j3 sr is considered to be of 
the same order as a, both arising from spin relaxation). If it comes purely 
from the momentum transfer, the wall resistance is estimated to be R w = 
3 x 10 -4 Q, a quite reasonable value. A striking point in this experiment is 
a significant enhancement of the effect of the force due to resonance, which 
made possible the low-current operation. On the other hand, the spin-transfer 
torque is suppressed in the MHz range (as seen from the factor of ia in the 
spin-transfer torque term of F(t)). 

Quite recently, Thomas et al. [58,90] succeeded in detecting periodic oscillation 
of a wall in a confining potential by using a ns current pulse at j = 6.9 x 
10 11 A/m . The motion was consistent with the rigid wall description in terms 
of X and 4>q. Periodic variation of chirality, </>o, of a wall was observed in the 
presence of magnetic field and current pulse of 10ns at 1 x 10 12 A/m 2 [91]. 
The results indicated that the chirality, <po, plays an important role in the wall 
propagation, as predicted theoretically [44,26,48]. 



4 Localized spins 

4-1 Spin Lagrangian 



Throughout this paper, we use the Lagrangian formalism, as it is useful in de- 
scribing collective objects such as domain walls. The basics of the Lagrangian 
spin formalism are summarized in §A. The localized spin S(x, t) is a function 
of position x and time t. We define polar coordinates as (Fig. 11) 

S = S(sin 9 cos 0, sin 9 sin 0, cos 9) . (20) 



We will consider the effect of damping later. The Lagrangian of the spin system 
with Hamiltonian H% is given by 

^-hSJ>(cos6-l)-H s , (21) 

a 6 

The first term, 

L B = f^hS Mcos 9-1), (22) 

J a d 

is a term that describes the time development of the spin. This term has the 

geometrical significance of a solid angle in spin space, or a spin Berry phase 
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Fig. 11. Polar coordinates (0, (f>) of localized spin. 

[92]. As is easily seen, this Lagrangian (21) results in the correct Landau- 
Lifshitz equation of motion: 



where 



HS = k-iBs x S, 

tff. 



(23) 



(24) 



h-f ss ' 

is the effective magnetic field acting on the spin (from the localized spin Hamil- 
tonian ifg). (The energy density due to field is given as —MBs = -^S-Bs) 

Let us demonstrate this. The equations of motion for the spin derived by 
taking derivatives of Lg with respect to 9 and <f> are 



d 5Ls SLc 



dt 59 
dSL s 



89 



dt 



(25) 
(26) 



By use of the explicit form of Lb , they become 



HS sin i 



-HS sin 99 



5_Hs 
' 59 



(27) 
(28) 



Eqs. (27) (28) become more concise in terms of vector equations for S. The 
time-derivative of S is given by 



S = S(9e e + sm9<pe ft 



(29) 



where 
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cos 9 cos (t 



eg 



cos 9 sin < 
— sin# 



(30) 



( ■ A 

- sin ly- 



cos < 




(31) 



/ 



are unit vectors in the 9- and ^-directions, respectively (Fig. 12). For any 




Fig. 12. Unit vectors n, eg,e^. They satisfy n x eg = e^, n x e ? 
eg x e^ = n. 

function F(9,(f>) of (9,(f>), the following identity holds: 



-eg, and 



5F / 5F\ 
l6 =SBm9 [ e *'SS. 



Thus Eqs. (27) (28) (29) become 



(32) 
(33) 



HS = S \e e \e 4 



SH* 



5S 



e-s e e 



5S 



s ¥ xn 



which is Eq. (23). Here we used 



eg x es = n, 



(34) 



(35) 
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Fig. 13. Motion of localized spin under an effective field, B, without and with 
damping a. Without damping, the spin precesses forever, while with damping, the 
spin direction eventually relaxes to that of the field. 



where 



sin 9 cos ( 



S 



\ 



sin 9 sin < 
cos 9 



(36) 



is a unit vector along S. 



4-2 Spin algebra 



The meaning of the " spin Berry phase" term can be understood if one notes 
that the canonical structure is contained in this kinematical term in the La- 
grangian. Let us demonstrate this within classical mechanics. The canonical 
momentum conjugate to is defined as 



5U 



hS 7 — hs. 



(37) 



Defining the Poisson bracket (times h) by {A, B} PB = (8A/d(j))(dB/dS i 
(dB/d(j))(dA/dS z ), we have {0, SJpb = 1. By using S x ±iS y = 



IS 2 - SI e 
we can derive the correct SU(2) algebra of the spin angular momentum as 



±i<[ 



{Si, Sjjpb — etjkSk- 



(3* 



4-3 Dissipation 



In reality, the spins are subject to energy dissipation arising from various 
sources such as phonons. This damping is rather essential in spin dynamics, 
for instance, to relax to the equilibrium configuration (Fig. 13). The equation 
of motion with damping taken account is called the Landau-Lifshitz-Gilbert 
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(LLG) equation [93], which was discussed on a phenomenological basis. It 

reads 

dS „ „ a _ dS 1 . . 

Here a is a dimensionless damping constant, called the Gilbert damping pa- 
rameter. This a here includes all the effects of conduction electrons and non- 
electron origins. The electron contribution to a is calculated in §14 (Eq. (352)). 
The parameter a is not therefore a phenomenological parameter, but can be 
derived [28,29]. 

For the study of domain wall dynamics, the effect of damping is incorporated 
in the Lagrangian formalism by use of Rayleigh's method as in classical me- 
chanics [94] by introducing a dissipation function 

Ws = J ^^ 2 = J ~hS (9* + sin* 6^), (40) 

where a is a dimensionless damping parameter. The equation of motion with 
damping included is modified to be 

±5_Ls _ 5Ls = J_Ws (41) 

dt 5q 5q 5q 

where q represents 9 and <fi and the last term describes the energy dissipated. 
Explicitly, they read 



^ 59 59 

-HSsm99 = -^-—^, 43 

5(j) 5<p 

In terms of vector equations for S Eqs. (42) (43) read 

Noting -j^ = %£, we see that the Landau-Lifshitz-Gilbert equation (39) is 
obtained. 

Let us see that the damping term really corresponds to energy dissipation. By 
multiplying Sx and then S- in Eq. (39), we obtain immediately the energy 
dissipation rate proportional to aS 2 ; 

The above procedure of introducing dissipation seems artificial. We will later 
show that Gilbert damping can be derived microscopically by considering a 
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whole system including spin relaxation. Actually, we will calculate the damp- 
ing arising from the exchange coupling to conduction electron and spin relax- 
ation [29]. 



4-4 Spin Hamiltonian 



In this paper, we consider a spin system with an easy axis and a hard axis, 
chosen as the z and y directions, respectively. The energies of one spin pointing 
in the easy and hard axes are denoted by — \KS 2 and ^K±S 2 , respectively, 
where K and K± are positive. With this two-axis anisotropy, the analysis 
here can be applied to both wires and films with perpendicular anisotropy 
(Fig. 14). We also include a pinning potential, V P i n . The explicit form of V p i n 
will be discussed later. The Hamiltonian is thus 




muM^Utf , 




> 



Fig. 14. Systems we consider: wires with easy axis along the wire and film with 
perpendicular easy axis. The domain wall become the Neel type and the Bloch 
type, respectively (see Sec. 4.6). Coordinates in spin space are shown (which do not 
necessarily coincides with spatial coordinates). The current is applied along the wall 
direction. 



Hs = I $ (i (V5)2 " j {Sz)2 + ir ( ^ )2 ) + v ^ 

^ (^((W) 2 + sin 2 £(V0) 2 ) + ^f- sin 2 0(1 + k sin 2 




where k = K±/K. The equation of motion Eq. (39) is then given by 



- sin 9<p-a6 = ^- (\ 2 (-V 2 9 + - sin 20(V0) 2 ) + - sin 29(1 + k sin 2 <f>)\ 

+J_^w (47) 

HS 69 y } 

9-asin9<j) = % l — f-A 2 V(sin 2 9V<p) + - sin 2 9 sin 20) 

n sm9 \ 2 

+ tJ—J^- (48) 



HS sin 9 



31 



Here we have introduced a length scale determined by 

(49) 

This length governs the spatial scale of the magnetic structure and turns out 
to be the thickness of the domain wall. 

4-5 Static domain wall solution 



Throughout this paper, we consider a planar (i.e. one-dimensional) domain 
wall as realized in a narrow wire or film, where the spin configuration depends 
only on the coordinate z in the wire direction. 

We consider in this subsection the case without magnetic field and pinning. 
Equations (47) (48) for a static configuration of a constant 6 are then given by 

^| = J-siii0cos0, (50) 

dz A 




which yields after integration 

2 1 
= —sm 2 9 + C, (51) 

where C is a constant. For the configuration we consider, sin 9 —* at z — ► 
±oo, we see that C = and thus 

dz A 

This equation is easily integrated to obtain tan | = e T( - z ~ x ^ x , where a constant 
X represent the wall position. The solution for 9 is thus obtained as cos 9 = 
±tanh^^. The sign here corresponds to a "topological charge" of domain 
wall. We consider in this paper domain wall with positive charge, i.e., spin 
texture changing from — z at z = — oo to +z at z = oo. The wall solution we 
consider is thus given by 



cos# = tanh — - — (53) 

A 



1 
cosh 



^ 9 = ..^z-X ( 54 ) 



A 

with 6 = 0. 
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4-6 N eel wall and Block wall 



We considered above the case of a domain wall where magnetization is chang- 
ing in the spatial coordinate z, which coincides with the spin easy axis. The 
wall structure looks in this case as in the left of Fig. 14, and is called the Neel 
wall. Another type of domain wall, called the Bloch wall as depicted in the 
right of Fig. 14, is also possible. This corresponds to a case where the spatial 
coordinate along the wire is y, and the easy plane (zx-plane) is perpendicular 
to the wire direction. These differences of structure do not affect the electron 
transport or the spin torque in the absence of an interaction that correlates 
the spin space and real space, such as spin-orbit interaction. 



4-7 Pinning potential 

Let us specify here the pinning potential. Pinning arises from various origins. 
Here we consider a simple case of a point defect that modifies the easy axis 
magnetic anisotropy, K. Other cases can be treated in a similar manner, and 
the result of pinning potential for X is essentially the same. A defect at x = 
is assumed to cause local enhancement of easy axis anisotropy of 5K(> 0), 
and then the pinning potential is given by 

V pin = -[ ^5K^5\x) sin 2 0(x), (55) 

J a 6 2 

where 5 3 (x) = 5(x)5(y)5(z) represents a 5-function in three dimensions. In 
terms of domain wall variables, this potential reads 

V + {X) = -—«£*?■ (56) 

For generality, we model this potential by a harmonic one given by 

v Pm = ^M w ^ m (x 2 - emz - i*d = ^V 2 - em - i*d, (57) 

where 6{x) is a step function, V = 5KS 2 £ 2 /(2N W ) is the pinning strength 
per spin, M w = £^ 3 = j^g is the wall mass, and fi pin = y/2K±V -^ 
corresponds to the oscillation frequency at the potential minimum. The range 
of the pinning £ is equal to A if the defect is point-like. 



33 



5 Description of rigid planar wall 



5.1 Collective coordinates 



To derive the equation of motion of a rigid domain wall, we here consider 
the collective coordinate description [95,96]. This treatment and the results 
are essentially the same as the one considered by Slonczewski [44,97] in the 
context of dynamics under a magnetic field. 

The idea is to consider the constant X in Eq.(54) as a dynamical variable, 
X{t). During wall motion, the angle 4>(z,t) can be excited, as suggested by 
the equation of motion, Eq.(39) and Eqs. (47) (48). Another collective variable, 
</>o, defined as the spatial average of </>(x) [31], is therefore required to describe 
the wall dynamics. In fact, the features of spin dynamics are taken account 
of by the fact that X and 0o are canonically conjugate to each other, as 
indicated by the fact that the first term of Eq.(21) is written as oc X<p . 
Naively considering only the variable X as dynamical results in the wrong 
answer in general. These two variables are called collective coordinates since 
they describe the collective dynamics of many localized spins. 

In the absence of sample inhomogeneity and a driving force, X describes a 
gapless zero mode owing to the translational symmetry of the system. If the 
transverse anisotropy energy K± is zero in addition, <fi is also a gapless mode. 
We will in this section consider this case (Vo = 0, Kj_ = 0) and show that 
the two variables X and 0o indeed become dynamical variables when we take 
account of the fluctuation around the classical solution. 



5.2 Coherent representation of spin 



To describe fluctuations (spin waves) around the domain wall, coherent repre- 
sentation of spin is convenient. We could instead define a fluctuation in rather 
a naive way as 9 = 9 + 59 and = </> + 5(f), as done in Ref. [31], but this 
definition results in a rather complicated spin-wave Hamiltonian. So let us 
here use a coherent representation and define a complex variable £(x,t) as 

£ = e**tan|. (58) 

This variable is related to 9, d> as 
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sin v cos < 



l + l£P 



sin 6 1 sin = — i 



COS0: 



1 + KI 2 

i-iei 2 



i + iep 



-2 



^=^^u-^e-o- (59) 

Thus the spin Lagrangian is written as (without pinning and K±) 

Ls = 1 l^TTW h U ~ U) ~ (ITW ( A2|V?|2 + 1?|2 )) ' (60) 

We now study fluctuation around a static domain wall sotultion, 



, Z ~ X 
cos V = tanh — ; — 



(61) 



where X and 0o are arbitrary constants, since we here assume vanishing ex- 
trinsic pinning and hard-axis anisotropy. We define a complex variable r\ rep- 
resenting fluctuation as 

t = e -u(z,t)+i<f> +ri(z-X,t) /Q2) 

where u(z,t) = ^-^-- The fluctuation r](z — X,t) here is thus defined with 
respect to the wall position. The derivatives of £ are 



d t Z = (d tV )C (63) 

Writing rj = rfe + ir]i, we obtain 

^-^ = 2i7Zr|e| 2 

A 2 |V^| 2 = [5 hz (1 - 2XV zVR ) + A 2 |V^| 2 ]|e| 2 . (64) 

We expand the Lagrangian (60) up to the second order in r\ by use of 
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lei 2 i i / ^ 2+ , 

( \ + 2t7 R tanh u + f 2 ^-) t/r 2 ) , (65) 

V cosh uJ J 



(1 + lei 2 ) 2 4 cosh 2 u 
and obtain 



^s = ~hS j ^— ^ M*) - KS" J *£ *_^| V„| 2 - H w , (66) 
J or cosh u J a 6 cosh u 2 

where H w = N W KS 2 is the domain wall energy. Since H w is a constant, we 
neglect it below. In Eq. (66), linear terms in 77 vanish since 77 is defined around 
a classical solution, Eq. (61). The fluctuation 77 is not a proper mode since the 
cubic terms have a weight factor — K— in the integration. This is removed by 
redefining the fluctuation as 

77 = 2fj cosh u. (67) 

Then we have, for example, 



dtf) = 2 cosh udtt] 

<9j77 = 2 cosh u ( <9j7/ + tanh u—5 i:Z fj J . (68) 

The Lagrangian is then written as (after symmetrization of the time- derivative 
term) 



d 3 c 



L s = f ^ \ihS(m - m) - 2KS 2 ( A 2 | Vt?| 2 + ( 1 L-) |t?| 2 

J a J L V V cosh w/ 

The dispersion of fluctuation is determined by the quadratic term. 



.(69) 



5.3 Spin-wave dispersion 



The spin-wave part is written conveniently by use of eigenf unctions satisfying 

where u is the eigenvalue. The eigenfunctions are well known [98]. There is a 
single bound state with uj = 0, 

<*« = ds*' (71) 
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and continuum states labeled by k: 

uj = 1 + k 2 X 2 = u k , 
with 



V^fc(^) 



-ikX + tanh — | 



■tkz 



(72) 



(73) 



The wave function <£>o is called a zero mode after its zero eigenvalue. The 
most important feature of the fluctuation is that the bound state function is 
a derivative of domain wall solution; (po(z) = —XV z 9q. Other wave functions 
are orthogonal to zero-mode wave function: 



dzip k (z)ip (z) = 0, (74) 

and form an orthogonal base (dropping an oscillating term at \z\ — > oo), 



dzip k '(z)(p k (z) 



jL= r dz (e« k - k '>(l + X 2 kk') - \4- fe^-^tanh^l 



2n^Ju k 
■ 8{k-k'). 



(75) 



5.4 Zero mode 



By use of these eigenf unctions, we can expand 77 as 



V(x, t) = I -Vo(t)ip (z) + J2 Vk(t)<p k (z) J 



(76) 



(Here we neglect dependences of rj (t) and r] k (t) on the wave vector perpen- 
dicular to the wire direction, considering a narrow wire.) The Lagrangian is 
then written as 



iV a 



ihS 



(Wo - Wo) + Yl (iftS(r] k r] k - 77*77*.) - 2KS 2 uj k \r) k \ 



(77) 



The zero mode, 770, has an important role. This can be seen by rewriting Eq. 
(62) by use of fluctuation modes r] and r) k s; 



z=* 



-i(j>o 



cxp 



z 1 
2 cosh t-(-7/Wo + Yl Vkfk) 
x l k 



z-X- 

~e * 



¥^ e i(^+i™vo(t)) exp 



2 cosh - Y Vk<fik 
A k 



(78 
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Thus the real part of 770 corresponds to a translation of domain wall, X and 
the imaginary part to 0o- Therefore, taking account of zero-mode fluctuations 
is equivalent to treating X and 0o as dynamical variables, X(t) and 4>o{t), i.e., 

X(t)=X + XRe7] (t) 

0o(t)=0o + lmr/ o (t). (79) 

The spin Lagrangian is now written as the sum of the domain wall part and 
the spin wave part as 

L s = L® + L sw , (80) 

where 



4 O) = ^Wo-X0 o ) 



L sw = N W ^ (iftS(r} k f)k ~ VkVk) ~ 2KS 2 u k \r) k \ 



5. 5 Condition of rigid wall 



When K± is finite, the definition of 0o is determined by the condition of 
vanishing linear coupling to fluctuations [31]: 

dz 



/CLZ 
— sin 2 # o 0M) (82) 



where sin#n 



cosh A ^ 



If a pinning potential Vq is present, the energy scale of X motion is Vq. Sim- 
ilarly, the energy scale of the 0o- m ode is given by K±. Since the energy gap 
of the spin-wave mode is ~ ^KKj_, the modes described by X and <pQ are at 
low energy compared to others if the following condition is satisfied: 



V < y/KK ± , K ± < K. (83) 

In this case, the low-energy wall dynamics is described by the two variables, 
X and 4>q. Otherwise, the pinning and/or Kj_ leads to deformation of the wall, 
whose description requires other variables than X and 0o- The condition (83) 
gives a criterion that such deformations can be neglected. 

Precisely speaking, we need one more condition for justifying collective de- 
scription, namely, vanishing linear coupling of spin-wave modes to X or O - 
In reality, when Vq an d K± are finite, such linear couplings arise, and the wall 
dynamics is not closed in X and 0o in a strict sense. This is quite natural 
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since the pinning and K± result in a deformation of the wall whose descrip- 
tion requires other variables than X and <j>q. However, the condition (83) also 
assures that such linear couplings are small. We assume the condition (83) in 
this paper. 

The wall solution we start from is thus given by 9 = 9o(x,t), where 

i z-X(t) 

cos 9 = tanh 

A 

sing o^ _ m , (84) 

cosh — y^- 

and 

<P = Mt)- (85) 

5. 6 Rigid domain wall Lagrangian 



From these considerations, the Lagrangian for the low-energy dynamics of a 
rigid wall is given by using n = (6*o,0o)- As a result, L s describing many 
spins (with Hamiltonian given by Eq. (46)) reduces to Lffl of two dynamical 
variables (^ denotes without electrons): 



4°> = hN w S (j<p ~ ^sin 2 o ) - Vpi>o], 



(86) 



where A^ w = 2XA/a 3 is the number of spins in the wall. (A is the cross-sectional 
area of the system.) 

The equation of motion of the wall is obtained simply by taking variations 
with respect to X and <po, including the dissipation function, 

aN w hS (X 2 • 2 \ 
^ = ^—(^ + 00 J- (87) 

In the absence of electrons and a magnetic field, the equation becomes (using 
Eq. (41) and Eq. (57) for pinning) 



^o + a-- -^xe^-\x\) 



X 2AVb 
A ~~hStf 2 

X -a\<i) = ^—Ssm2<f) . (89) 

2n 

The solution without any driving force is of course trivial: <$$ = and X = 0. 
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5. 1 Domain wall particle ? 



The equation of motion (88) (89) is one example of the Hamilton equation 
of motion for a classical object. For a classical particle with mass M, the 
Hamilton equation for the position X and momentum P is given by 



P = F (90) 

X = — (91) 

M V ; 

where F is the force. For a domain wall, 0o plays the role of momentum P, 
and so the relation between X and momentum is not simply Eq. (91) but is 
non-linear (Eq. (135)). In this sense, the domain wall is not a classical particle. 
This is natural, since the domain wall (even in the rigid case) has an internal 
degree of freedom, O - 

Let us consider the equation of motion with harmonic pinning, 



*°+4^( F - 2Ar -| A ')' (92) 

X- a \4, = ^Ssm2^ + jr^gT. (93) 

Here we included force F and torque r (from magnetic field or current), which 
are constant in time here. 

There is a limiting case, where the domain wall behaves like a particle de- 
scribed by one linear equation of motion. This occurs if 4>q remains small, i.e, 
if the hard-axis anisotropy is strong. 



5. 1. 1 Weak extrinsic pinning 

Let us first consider the case of large Kj_. In this case 0o ~ and so we can 
linearize the shape anisotropy term as sin20 o ~ 20 o . Then we obtain from 
(93) (neglecting a 2 < 1) 

^-i6rs^- a i^s{ F ~ 2N 4 x )->irs^ (94) 

Substituting O i n Eq. (92), we obtain the equation for X as 

M W X + M w —X + M w Q pin 2 X = F, (95) 
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where 



*.= ££ (96) 

°+-JW (97) 



= -(^ + 25^72 , (9? 



are mass, pinning frequency, and reciplocal decay time of the wall, respectively. 
In this limit, the domain wall becomes a simple particle. The existence of mass 
of the wall was first discussed by Doring [99]. 



5.7.2 Strong pinning or large <po 

When the pinning potential is harmonic, we can always delete X from the 
equations of motion (92) (93) and obtain one equation for O as 

M+io + ahSN w 4> (l + -^i cos 20 o ) + N W K ± S 2 ^ sin 20 o = -r, (99) 



where 

n 2 s 2 e 

2V X 2 

is the mass of the 0-particle. The <fi variable thus is a particle moving in a sine 
or cosine potential [30,31]. 



M <$> = ^wi^tt^ (100) 



Comparing two masses M w and M^, we see that M w is lighter than M^ when 
pinning is weak, Vq <C K±. This indicates that X is a better variable to 
describe the weak pinning case. In contrast, M w is heavier if pinning is strong, 
Vq ^> K±, and in this case, the system is better described in terms of the 
dynamics of <fto [30,31]. This applies at the quantum level, too. In a path 
integral representation, the vacuum transition amplitude for the Lagrangian 
Eq. (86) is given as (approximating (f) <C 1) 



^ 



-*2l f d t(ns4<t>o-^- sin 2 ^o-^-X 2 , 

VXV(/) e h J V xV0 2 «* J. (101) 

When the hard-axis anisotropy is strong, 0o is a high-energy mode, and should 
be integrated out first, resulting in 

Z w = |Me4/4f^-f^ 2 ). (102) 

Thus the wall behaves like a standard particle. When pinning is strong, we 
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can integrate out X and obtain 



Z w = jV<p e »-»™V» w 2 ^7. (103) 

In this case, 0o (chirality of the wall) becomes a good variable, and tunneling 
of chirality can occur as discussed in Refs. [30,31]. 



6 Conduction electrons 



In this section, we consider conduction electrons. We use the second quantized 
representation, and write annihilation and creation operators for the electron 
with spin a — ± at site x and at time t by c a (x, t) and c^ a (x, t), respectively. 
The electron density is c f (x)c(x). We sometimes suppress spin indices when 
obvious, e.g., c^c denotes J^a^aCa- The Hamiltonian for the electron we con- 
sider is given by a free part H , an impurity scattering H imp , a spin flip part 
H ST , and an interaction with the electric field H em which drives the current. 
(The interaction with localized spin is explained later in §6.5). The electron 
Hamiltonian is thus 




H c = J d 6 x I 2^|Vc|' - e F dc ) + H imp + H sv + H em , (104) 

with €f being Fermi energy. 

6.1 Free electron 



The free part, H , is given in Fourier space as H = Yjk e kCk c k where e^ = 
■^ — €f and the Fourier transform is defined as c{x) = -4^J2k elkXc k- Let us 
here briefly examine the behavior of the free electron Green's function. Con- 
ventional Green's functions such as retarded and advanced Green's functions 
are not physical quantities, but are useful in calculation. Let us consider the 



behavior of a function p(x,t,x',t') = (c^(x',t')c(x,t)) , which is the elec- 
tron density if at equal times and equal positions. The average (• • ■) H is for 
free states determined by H . This function satisfies an "equation of motion" 



d tP (x, tx', t') = -- (cV, [#o, c(x, t)] /H 



i(^7 2 x + 6 F /h\p(x,t,x',t'). (105) 
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This equation is in fact a continuity equation, but is not convenient for cal- 
culating non-free cases. A more useful function can be denned if one includes 
anti-commutator and step function in time as 



<f{x, t, x', t') = --9(t - ({cV, 0, c (*> t)}) Ho ■ (106) 



The equation of motion is then 



(ihd t + ( £-Vl + e F /h j j g T (x, t, x', t') = 5{t - t')5 3 {x - x'), (107) 

which is a Green's function with a source term. This type of equation is math- 
ematically useful allowing systematic calculation. This function is called a 
retarded Green's function due to a factor 6{t — t'). The equation (107) is eas- 
ily solved using the Fourier transform, g T (x, t, x', t) = y J 7^Y,k e ^ k x ~ dki^) 
as 

£(«) = l — ^ (108) 

to — e k + tO 

where the small imaginary part denoted by iO is to reproduce the retardation 
factor, 6{t-t'). 



6.2 Impurity scattering (spin-independent) 



We include here the effect of electron scattering by impurities. The scattering 
is treated as elastic and spin- independent, and is then described by a term 



•Wimp 



H imp = Y^ / d 3 xv(x - Ri)c\x)c(x) 

i=l 
N- 

1 " imp 

= EE^)e i(fc - fc ' ) - R *ct fc , Cfc , (109) 



t=l kk' 



where v represents the potential due to an impurity, Ri represents the po- 
sition of random impurities, A^ mp is the number of impurities, and v(q) = 
y J d 3 xe tq ' x v(x). We approximate the potential as an on-site type, v(x — Ri) = 
Vi mp a 3 S 3 (x — R^ (fi m p is a constant), i.e., v(q) = ^p (N = V/a 3 is number 
of sites) . To estimate physical quantities, we need to take the random average 
over impurity positions as 
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E (*-** > = E / ^e~^ = N imp 5 q , = (for g ^ 0) 

i i 

y^ / e -iq 1 'Ri e -iq 2 'R j \ _ V^ f " Rj f d Rj c -i qi .R lc ~iq 2 -R, 
ij ij 

— N imp qi+q2 fi + -Wi mp 0qi,o0q2,0 = ^imp0qi+q2,Cb (HO) 

where the term linear in iVi mp is a contribution at % — j. Taking account of 
successive impurity scattering by ladders, the electron Green's function in the 

presence of random impurities is given by e.g., g\{oj) = ^— -, where the 

inverse lifetime is given as 



ij 91 92 



v f 1 ) I-E E (e^^e— ^)^ + q» - f wL p ,a 3 . 



(Ill) 
Here z/ is the density of states per volume (v = N(0)/Na 3 , N(0) is the density 
of states) and rii mp = — j^p is the impurity concentration. In the derivation 
we used J^kQki 00 ) — ^(0) J-^oo ^ e _ rig — ~~ iftN(0). We consider a conducting 

e +2r 

system, and so 

^>1. (112) 

When spin polarization due to (uniform component of) localized spin is taken 
into account, the lifetime becomes spin- dependent. For the most part of this 
paper, we neglect this spin-dependence, to avoid unnecessary complication. 
The spin-dependence is correctly included for evaluation of delicate quantities 
such as a and (3 below in §14. 



6.3 Spin relaxation 



We also include spin flip interaction due to random magnetic impurities to 
take account of the spin relaxation effect. The interaction is represented by 

H ST = u s J d 3 x J2 S imPl 5(x - R! t )^*c) x , (113) 

i 

where S[ mpi represents the impurity spin at site R^. (The spin-orbit interaction 
leads to essentially the same results as spin flip case here.) A quenched average 
for the impurity spin was taken as 



^imp, (a, j3 = x,y) 



SnxpfStap^Ma/* X -^ ■ ( 114 ) 

Simp^ (a, P = z) 

The effect of spin relaxation is considered in §14. 
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6.4 Electric field 



In this study, the applied electric current has the most important role of 
inducing spin dynamics. The effect of current is calculated as a response to 
the applied electric field. The interaction is expressed by use of uniform charge 
current density j and electromagnetic gauge field A em as 

H em = - d 3 xA cm • j, (115) 



The gauge field is given by use of E = — A as 

E 



A - n p iQ 0t 



(116) 

U ->0 



where E is the applied electric field assumed to be spatially homogeneous. For 
calculation purpose, it is treated as having finite frequency, Qq, which is chosen 
as Qq — > at the last stage of the calculation (this is a standard technique in 
linear response calculation [100]). The current density is given in the presence 
of the gauge field by (e < is electron charge) 



m k 



1 r 7 o eh i , + £♦ . e 2 h sr ^ + 

- / d x—-{c l V, c) A em: „ > d k c k 

V J m 2 m '-r 

Y^ (hk v - eA emiU ) 4c fe . (117) 



Within the linear response to E, the last term is neglected. We used the 
coupling (A ■ j) between gauge field and current, and not the one between 
electric charge and scalar potential ($p), since the A ■ j-description is known 
to be convenient in describing the system as spatially uniform as in the Kubo 
formula. (In contrast, $p description is useful in a Laudauer type description 
treating the spatial difference of chemical potential explicitly.) 

The electron part of the Lagrangian is given by use of the Hamiltonian H e as 

L e = ih f d 3 xc<c-H e , (118) 

where the first term is a dynamical term that correctly reproduces the Schrodinger 
equation. 

6.5 s-d exchange interaction and adiabatic condition 



The most important interaction for us is the exchange interaction between 
electron spin and localized spin, which is a source of all current-driven spin 
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dynamics and magnetic transport properties. This term is given by 

H sd = -J sd f d 3 xS ■ (cW) = -M f d 3 xn ■ (J ere), (119) 

where J sd is the strength of s-d exchange interaction and M = J sd S is half the 
exchange splitting. An important point is that J sd is rather strong in 3-d fer- 
romagnets: J s d/^F ^ 0(0.1 — 1). These values are indicated from experimental 
observations of large magnetoresistances such as the GMR. 

The most non-trivial part of the theory is the treatment of this strong ex- 
change interaction when the localized spin has a spatial structure and/or is 
dynamical. Fortunately, spin structures in 3-d ferromagnets are slowly vary- 
ing compared to the scale of conduction electrons. This is a consequence of 
the strong exchange interaction, J, between localized spins, which is of order 
of 1000 K as indicated by the high critical temperature of 3-d ferromagnets 
(For Fe, T c ~ 1043 K). (Correctly, the typical length scale, A, is determined 
by the ratio of exchange energy and magnetic anisotropy (Eq. (49)).) Since 
many localized spins within the scale of A are coupled, the spin structure is 
(semi-) macroscopic and its time scale is slow compared to that of electrons. 
From these considerations, the electron can go through the spin structures 
adiabatically. 

In our study, we carry out the expansion with respect to the gauge field rep- 
resenting the non-adiabaticity. The gauge field contains the space and time 
derivatives of localized spins. The expansion parameters are given by 



e F 1 £ 1 e F 2 r 

MkpJx ~ {k F X) 2 KM ^ 

-j=-<l, (120) 

where u w is the frequency scale of the domain wall motion. This is understood 
by noting that the electron spin density lowest order in the gauge field is given 
by Figs. 19 and 20, and the correction to these processes contains a factor of 
-\ (fe • A + )(k ■ A~)g k+ g k _ or Aq Aq g k+ g k ^ in the adiabatic limit, where g 
is either the retarded or advanced Green's function. Noting A^ ~ \ViS\ ~ 
A -1 , A ~ u w , and g kj+ g k - ~ «^j when Mr/h ^> 1, we have identified the 
expansion condition as given by Eq. (120). The gauge field expansion is thus 
justified if either the spin splitting is large or the spin structure is slowly 
varying, and it results in a different series expansion from the simple gradient 
expansion assuming (kpX)' 1 ^C 1 and u> w r ^C 1. (In Ref. [65], the gradient 
expansion condition was argued to be A st /A <^ 1 and \ st u w /v F <^ 1, where 
they introduced a phenomenological parameter of spin transport length scale 
A st [47]. However, the result (Eq. (39) of Ref. [65]) appears to be that of a 
simple gradient expansion assuming (/ci?A) _1 ,a; w r ^C 1. ) 
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In actual 3-d ferromagnets, M/e F < 1 and therefore the gauge field expansion 
would become essentially the same as a gradient expansion. Nevertheless, it 
would be formally useful to carry out the gauge field expansion first and then 
consider a slowly varying limit as we will do in §9.3 in estimating the reflection 
force. 

The adiabatic condition obtained above has an extra factor of £/X due to 
disorder scattering if we compare with the condition proposed by Waintal and 
Viret [62]. They obtained in the ballistic case the adiabaticity condition of 

1 €F l «1, (121) 



k F XM "" (k F+ - k F ^)\ 

where the left-hand side is a ratio of the precession time of conduction electron 
due to the exchange interaction, h/M, to the time needed for the electron to 
pass through the spin structure, X/v F . 

In the context of quantum electron transport, Stern [101] introduced a different 
condition in the disordered case, 

h/(Mr) < 1. (122) 

This would be satisfied in 3-d ferromagnets, but is not a necessary condition in 
our calculation. For quantum transport, other conditions have been proposed 
[102], which appear to depend on the system considered. 



7 Equation of motion of domain wall under current 

7.1 Effective Lagrangian 



As we have discussed, the total system we consider is described by the La- 
grangian L = L s + L c - H sd (eqs.(21)(46)(104)(118)(119)). We are interested 
in localized spin dynamics, and for this purpose, we derive the effective La- 
grangian for localized spin by integrating out the electron. "Integrate out" 
here means taking the trace over the quantum mechanical states of the elec- 
tron. In the path-integral formalism [103], this process corresponds indeed to 
an integration in the following way. The partition function of the system is 
represented as 

Z= f V9V(f)sm9VcVce i I dt(Ls[S]+Lc[5 ' c] ~ Hsd[s '~ c ' c]) . (123) 

Here T>9 denotes integration over the field variable [T>9 = U Xit d9(x,t)), and 
c and c are Grassmann numbers corresponding to creation and annihilation 
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operators for the electron. The time integration is on the real axis, but the 
argument here applies also to the case of the Keldysh contour C. By integration 
over the electron, Z reduces to 

Z = t V9V4> sin 9e l I dtL f [s] , (124) 

where L c <f [S] = L$ + ALs is the effective Lagrangian for localized spin and 

fdtAL s [S] = -ihx f VcDce*! dt{L «fc c] - H * d[SAc]) = -i\nZ e , (125) 

is the contribution from electrons, which includes formally everything from 
the electrons exactly. The equation of motion of spin with all the effects from 
electrons included is then written as (neglecting dissipation) 

if = °- (126) 

Let us look into this equation in more detail. The electron contribution is 
written as 



SAL, 



— jVcVc(-J sd {c(Tc))e^ dt{L ^ c] ~ H ^ s ^ c]) 



ss 

= -Js d s (127) 

where we noted that the right-hand side of the first line is the definition of 
electron spin density, 

,t 



s(x, t) = (c\x, t)ac(x, t)) . (128) 

(We define spin density without the factor of \ representing electron spin 
magnitude.) The average here is taken using the full electron Hamiltonian, 
H c + H s d, namely taking account of background localized spin structure, elec- 
tric field, impurity scattering, and spin relaxation. Let us define the effective 
field from the electron, B e , as 

B e = --^J sd s. (129) 

nrf 

We now understand that the effect of the electron is taken into account sim- 
ply by adding the effective field from the electron, B c , to the equation of 
motion (44). The full equation of motion (126) thus reduces to (now including 
dissipation) 

f = 7BsxS - 7S x S .-| S xf. (130) 

The effective Lagrangian for localized spin is therefore given by 

L s ea = L s -J sd [d 3 xS-s. (131) 
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Note carefully, however, that if we solve for the electron spin density s (Eq. 
(170)) and put the result in Eq. (131), we obtain the trivial answer of no 
effect from the current. To discuss spin dynamics, we have to first derive the 
equation of motion regarding s as independent variable as S, and then apply 
the result of s. This is what we do below. 



7. 2 Equation of motion 

From the considerations above, the effective Lagrangian of the domain wall in 
the presence of electrons is given by Ls + (H s d). Replacing the localized spin 
direction n by the domain wall configuration no (whose polar coordinates are 
(0 O , 0o) (Eq- (84))) yields 

L w = hNS (jfa - ^S sin 2 O ) - V pin [n ] + M Jo?xn Q - s. (132) 

Noting 

r 

— -^- = sm9 e (t> - s = (n x s) z , (133) 

b<po 

the equation of motion of the domain wall under current is given by 



4>o + a j = j^g(F + F pin ), (134) 

X -a \(f)o = —-Ssm2(f) + — r z . (135) 

Here F pin = P 7x > an< ^ ^ ne f° rce an d torque due to electrons are defined 

as 



'-<£ 



-M I d 3 xV z n - s, (136) 



T = -Jd 3 x/^-\xS = -Mjd 3 x{n xs). (137) 

We stress here again that these equations contain all the effects of the electron 
without any approximation so far. We note also that this set of equations, 
(134) and (135), is essentially the same as those obtained by Berger [33,24]. 
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What is new and essential in the present theory is that we have formal but 
exact expressions of force and torque, which we can evaluate by a systematic 
diagrammatic method. 



7.3 Equation of motion from the Landau-Lifshitz- Gilbert equation 



The equation of motion of the domain wall can be derived from the LLG 
equation, Eq. (14), by using the domain wall solution including collective co- 
ordinates, Eqs. (84) (85). Here we neglect nonlocal terms for simplicity. (These 
terms are calculated in §9.3.) Using d^S = S((d^9o)eg + sin^c^o)^), and 
#o = ysin# , V z 0q = — isin# , the LLG equation reduces to (see also Eqs. 
(47)(48)) 



8 — a sin 9q4>q = sin 9 I — — a0 o J 

-j s sin H — ;rz— sin sin 20 o 



2eS\ J * u 2h 
sin 6*000 + &9o = sin 9q 0o + a— 



2eSX J u 2h 

Integrating over position x, we obtain 



X ■ a 3 K ± S 



a 3 KS 

/3 sr j s sin^ — sin20 o (l + /tsin 2 o ). (138) 



X a 3 

h + a T = 2eJ\ (3s ^ (139) 

which is Eqs. (134) (135) with force from the /3 sr term included (see §9.3). 

This derivation of the domain wall equation from the LLG equation is useful 
as far as only local torque is concerned, but is wrong when non-adiabaticity 
is to be considered. In this paper, we proceed based on equations 134 to 137, 
which include all the torques without local approximation. 



7.4 Spin conservation law 



Let us briefly look into the conservation law of spin. The equation of motion 
of localized spin is given by Eq. (130). The spin part of the effective field, Bs, 
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is described by Eq. (46) as (neglecting pinning) 



h^B s = - JV 2 5 - Ke z S z + K ± e y S, 



.'/• 



(140) 



and so the LLG equation Eq. (130) is written as 



S + V • J< 



J, 



S(.I 



■ js 



ha 3 



S X S + T S , 



where the spin current associated with localized spin is given by 

J4 = ^((V„S) x S) a , 



(141) 



(142) 



and the spin source or sink is given by anisotropy and Gilbert damping as 



a 



t s = -SxS-K 



Q Q 

— OyO Z 

Jx&z 
\ I 



-K, 



OyD Z 



^—S X Sy j 



(143) 



The equation for the spin density of the conduction electron, defined by Eq. 
(128), can be derived by considering its time derivative, 



ds i 
dt h 



Jct[H,c]- [H,c ] ]rrc 



(144) 



and evaluating the commutation relation with the total Hamiltonian H. Using, 
e.g., [c^ x 'C x >, c x ] = —c x 5(x — x'), we obtain 



1 ds 1 _ . J s d „ 

--T = V • 7 S + T—^S X S + T s . 

Idt 2 J ha 3 



(145) 



Here jf s is the spin current density (divergence here is with respect to spatial 
coordinate), defined as 



« = ^<^ v ^- 



(146) 



where A V ^B = A(V^B) — (V fi A)B, and t s represents relaxation of electron 
spin. 

Combining Eq. (141) and Eq. (145), we see that s-d exchange torques cancel 

3 

each other in the equation of motion for the total spin, S to t = S + \s, and 
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total spin current, J^tot = Js+ \j s - This is natural since exchange interaction 
is the internal exchange of angular momentum, which does not change total 
spin dynamics. The continuity equation thus becomes 



Stat + V ■ J 



5, tot 



-T S + T S . 



(147) 



This continuity equation is another representation of current-induced torques, 
where the torque due to current is included in the V • j s term [104]. 



8 Calculation of electron spin density 



In this section, the electron spin density is calculated. The calculation in this 
section is done for general spin structures, not restricted to domain walls. 



8.1 Gauge transformation 



Our task now is to calculate the electron spin density s(x,t). This is non- 
trivial, since the electrons are interacting with the background spin, which 
is spatially and temporally non-uniform. For estimating s, we consider first 
the free part of the electron with the exchange coupling. The corresponding 
Lagrangian (we call L° e ) is 



1 
h 



d 3 x[ihc j c- — |Vc| 2 



1m 



epc^c 



M I d 3 xn- {c ] rrc). (14 



As we discussed in §6.5, we are interested in the adiabatic regime, and the 
treatment using a gauge transform becomes useful in this case. 

The idea of a local gauge transformation is simply to diagonalize locally the 
s-d exchange interaction. This is always possible by choosing an appropriate 
2x2 unitary matrix U(x,t) such as 



W(x,t)(n(x,t) ■ cr)U(x,t) = a z 



(149) 



(Since the localized spin direction depends on position and time, the matrix 
U also is, i.e., U(x,t).) This transformation is implemented by choosing 



U(x,t) = m(x,t) ■ cr, 



(150) 
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where m is a real three-component unit vector given as 



e 







m = sin - cos 0, sin - sin <p. cos - 
V 2 r ' 2 Y ' 2, 



The matrix [/ satisfies ?7 2 = 1, or 



f/(a;,t)- 1 = f/(a;,t). 



(151) 



(152) 



This unitary transformation corresponds to defining a new electron operator 
a = (a+ja^)* (t represents transpose) as 



c(cc,£) = U(x,t)a(x,t). 



(153) 



Now the exchange interaction is diagonalized, but this redefinition affects the 
kinetic term of the electron. In fact, we immediately see from the identity 

<V0, *) = U(x, t){dp + U(x, ty l d^U(x, t))a = U(x, t)(d„ + lA^a, (154) 

that a gauge field appears, defined as 

Ap = -iU(x,t)-%U(x,t). (155) 



By use of (151), A^ is written as 



.4, 



(m x d^m) ■ a = A a a a 



H^Q.1 



(156) 



where summation over a = x, y, z is suppressed. The gauge field is explicitly 
obtained in a vector notation with respect to the spin index (a) as 



— dfj,9 sin — sin 9 cos fyd^ 



d/j.9 cos <p — sin 9 sin (pd^ 



\ 



A^ee + A% - A%n. (157) 



/ 



(1 — cos 9)8^4) 
Here the components of gauge fields are defined as 



A° = eo-A. 



1 



E^ 4 



± 



i^e.-A^-^Te^A^ 



-sm9d^ 



where 



(158) 



Aj = A«±iA» 



-e^ (±id M - sin0<9 M 0) = e ± ^(Aj ± L4j). 



(159) 
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Some useful relations between gauge field and spin vector are 



d ll n = 2(-A°et + A*e e ) = 2A U , x n 



(d lk nxn)=2(Ale 9 + A*e 4 ) 

d^n x d v n = An{A e At - A 9 ^) 



A fl = -nx d^n - $ BM n = -(d^e^ - d^n + d 4> e z 
where Q-q^ = |(1 — cos6)d^(j) = —A^ ■ n. 



(160) 



After gauge transformation, the electron Lagrangian is written in terms of the 
a-electron as 



(Tx 



h 2 



ti 2 

iha)a 1 Va| 2 + epa^a — Mcva z a 

2m 



+i— J2( aiA iVia - (Vi(f)A 
2m ; 



h 2 

a) A 2 a)a — hxvAoa 

2m 



(161) 



Defining a Fourier transform as a(x) = -7^Y^k a kG lk ' x and A"(q) = y J d 3 xe iq ' x A"(x) 
the expression becomes 



L° c = Y^ iha ka ( d t ~ efca) a ka - H A , 
k 

where 

tka = — e F - crM, 

2m 

is the electron polarized uniformly along the z-axis and 



(162) 



(163) 



kqa 



Y, h \J„ Ik + | j ■ A™(-q) ) a ] k+q a a a k 



+^Y, A n-q-p)At(p)a\ +q a k 

pi 



(164) 



is the interaction with the gauge field (the Greek suffix \x runs over x, y, z, t 
and i runs over x, y, z). Here 



Uh) 



h 



kA 



m 



(165) 



for [i — x,y, z, t. The interaction with the gauge field is shown diagramatically 
in Fie;. 15. 
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Fig. 15. Diagrammatic representation of interaction between domain wall and con- 
duction electron. Solid lines are electron Green's functions and wavy lines denote 
interaction with the magnetization configuration (e.g. domain wall). 

The Lagrangian (162) is very useful for our purpose, since the perfect adia- 
baticity (represented by uniform polarization in Eq. (163)) and non-adiabaticity 
(gauge field) are clearly separated (Fig. 16). 



* t \ 

4tf>H 



Fig. 16. By gauge transformation, electrons (denoted by their spin s) interacting 
with a non-uniform localized spin structure (S) are transformed to electrons uni- 
formly polarized and interacting with the SU(2) gauge field, A, localized around the 
inhomegeneity. Thus the non-adiabatic component causing scattering is separated 
from the adiabatic component. 

The electric current density, Eq. (117), is modified by the gauge transformation 
and the electromagnetic gauge field as 



e 2 h 



/ d 3 x — -(c f Vi c) A cm4 y2c\.c k 

J m 2 m , 



eh 



m 



E 



ha k a k + Y^ AtiqWk+q^a-k - eA em)i a k a k 



qa 



(166) 



Thus the interaction with external electric field, Eq. (115), is given by (Fig. 
17) 



-"cm / j 



iehE; ■ 



m 2 Q 



1 JClot 



E 



k I 



hala k + J2 Af(q)ai +q a a a k 



<\q 



+ 0(E 2 jl67) 



Qo^O 
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k + q 




Fig. 17. Diagrammatic representation of interaction with the applied electric field, 
E, denoted by dashed line. 

8.2 Electron spin density 



Calculation of electron spin density is carried out in the gauge transformed 
frame by the standard diagrammatic method. What we obtain there is the 
spin density in the gauge transformed frame, s, defined by 

s(x,t) = (a ] aa). (168) 

It is related to the physical spin density, s, as 

s(x, t) = (c ] ctc) = (a ] U ] (rUa] = 2m(m • s) - s. (169) 

Let us define 

s = s e e 8 + s z n + s^e^. (170) 

Each component is written in terms of components of s as 



8, = -\j:e**3± 

^ = -^E(T)^^, (171) 

where s^ = s x ±is y . Eq. (171) can easily be checked by use of m = — sin |ee + 
cos |n and cos 9e e + sin 9n = (cos 0, sin 0, 0) = (e^ x e z ) . 

Let us see how the force and torque on the spin structure (Eqs. (136) (137)) 
are represented in terms of sg and s^. Using Vn = e d V9 + e^ sin #V0 and 
n x eg = e<b, n x e^ = — eg, we obtain 



F=-M f d 3 x (sgV6 + S<i) sm9V(f)) = 2M f d 3 x (-s e A* + s^A e ) (172) 
t = M f d 3 x ( S(j) ee - sge^) . (173) 

We see here that s z , the perfectly adiabatic component, does not contribute 
to the force or torque. For the rigid wall we are considering, V0 = 0, only 
components F z and r 2 affects its dynamics, and we obtain 
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F z = -M f d 3 xs e V z 6 
t z = -M f d 3 x sin 6 s^ (174) 



Thus the force and spin-transfer torque on a rigid wall arises from sg and s^, 
respectively. 



8. 3 Electron Green's functions 

The spin density s(x,t) is defined in terms of an a-electron as 

s a (x, t) = (a\x, t)a a a(x, t) ' 



= — ihtv[cr a G < (x,t,x,t)], (175) 

where 

G< a ,(x,t,x',t') = l -{aK'{x\t')a a {x,t)), (176) 

is a lesser component of the Keldysh Green's function [105,106] and the trace 
(tr) is over the spin index. This Green's function is an extension of the standard 
equilibrium Green's functions (retarded and advanced Green's functions) to 
a contour-ordered Green's function defined on the complex time plane. (For 
more details see §B.) The contour-ordered Green's function is defined as 

G ay (x,t,x',t') = -UT c a ff (x,t)a)Ax',t')) , (177) 

where t, t' are defined on a contour on the complex plane and Tq denotes 
ordering on the Keldysh contour C (Fig. 18). The average here denotes both 
quantum expectation value and averaging over random normal impurities. 
This Green's function contains information on retarded, advanced, and lesser 
Green's functions. 



C 



[£ 



3-^ 



Fig. 18. Keldysh contour in the complex time plane used in non-equilibrium Green's 
functions. 

In the free case, the lesser component is defined as 

g< a ,(x,t,x',t') = ^(aK'(x',t')a a (x,t)) Q = g< al (x-x',t-t'), (178) 

where ( ) denotes the expectation value by use of the free Hamiltonian, H , 
but in the presence of impurities (H imp ). By standard ladder summation over 
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multiple impurity scattering, the free lesser component is obtained in the 
Fourier space as 



*£(«)= /(«)(£(«) -ifcCa;)), 



where the advanced and retarded Green's functions are 



(179) 



9W) 



u - e k - i; ' 



2r 



;iso) 



and g T = (<? a )*. 

Calculation of s is done by standard perturbation theory with respect to the 
gauge field A^, up to the lowest, i.e., linear, order. We define the Fourier 
transform of a Green's function as 



Gka,k+q,(i> (t, t ) = — % \Tca ka (t)a k+q,a'{t ) 

= - fd 3 x fd 3 x'e- i(k - x ^ k ' x ' ) G (T)al (x,t,x',t'). (181) 



The full Keldysh Green's function defined on a complex contour satisfies the 
Dyson equation similar to the conventional time-ordered and retarded Green's 
functions: 



ihd t G kkl (t,t')=8(t - t') ({(*,,(*), atfc/CO}) + ^ (T c [H,a k (t)}ah,(t')Yl82) 



Here H = H + H cm + Ha (H imp is already included in r of free Green's 
functions, e.g., in Eq. (179)). Evaluation of the anti-commutation relation can 
be done using [AB, C] = A{B, C} — {A, C}B for any three operators A, B, 
and C as 



[h, a k (t)} = - y: 



q<\ 



h 2 



J2hJJk+l) A«{q,t)a a a k+q {t) 



+ ^E A t( ( l-P^) A t(P,t)ak +q (t) 

p,i 



;i83) 



We then obtain 
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ihd t G kk '(t, t') = 5(t - t')5 k ,k> + e k G kk ,(t, if) + ti£ E J n ( k + f ) A "(<1> t)a a G k+q , k ,(t, if) 



q /MX 



ieEj 



~ E E 4 Q (<? - p, t) A q (p, *)gw(*, + o(e 2 ). 

Am qp ai 



By multiplying [ihd t — e k ] l = g k , we obtain the Dyson equation on contour 
C, 



184) 



G kk '(t, t') = g k (t - lf)S h ,h> + ^E / dhg k (t - ti 



E J, (k + |) Aftq, h)a a + E ^^° 41 (^,,0 + E 4 ? fo *iK 

2 >i2\ 



G k + q , k i(ti,t J 



+ 0(£ 2 ,1 



0b(* - 0**,* + ^E J m (^y^j J c dt i9k(t - ti)A%(kf - MiK#fe'(*i - £') 



ieh tt-^ _ . ( k + k r 






./c 



c ./c 



^'^(fc' - fe, ^(i - tx)a a gw{tx - t 2 )g k ,(t 2 - if) 
+he in ° tl A°(k' - k, t 2 )g k (t - h)g k (ti - t 2 )a a g k ,(t 2 - if) 

+0(E\A 2 ). 

Here we neglected a trivial term which does not contain A. Green's functions 
are written by use of matrix representation in spin space, 
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G kk ' 



G k -,k'+ G k -^ k r 
and free Green's function is diagonal, 



(186) 



9k 



9k+ 
g k - 



187) 



The Equation (185) is for Green's function defined on complex time contour, 
and is not physical. To evaluate spin density, we need to take a lesser compo- 
nent, G < , of the contour ordered Green's functions. As described in §B, lesser 
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component of product is mixed with a lesser component and retarded and 
advanced components as 



c 



dtxAit.t^Bit^t') 



< poo 



/oo 
dhiA'it, t 1 )B<(t 1 ,t') + A<(t, tjB^h,?)), 
-oo 

(188) 
while retarded and advanced components are the products: 



/ dt 1 A(t,t x )B(ti,l/) = / dM r (Mi)-B r (ti,t')- 

JC J J-oo 



189) 



We define Fourier transform of gauge fields as 



A$(q,n) = / dfe-"^,*). 



(190) 



and electron Green's function as (r denotes <, r, a) 



Cfc,fc'(M ) 



do; r dui' 



iujt „iu. ' t , /-it 



2vr J 2tt 



-e e 



CLfc'^V) 



'1911 



or GX ]fe ,(u;,u/) = J dt J dt'e iuJt e- iuJ,t 'G T kk ,(t,t'). In general cj ^ a;', electron 
absorbs energy (frequency) from the electric field (and from spin structures if 
dynamical). Taking the lesser component, Fourier representation of Eq. (185) 
is written as 



Gfc, fc+ q(w, ^ + ft) = 0faA,ofo,o + ^H ^ ( fc + 9 ) ^(<?, ft)^^^fc+q^+n 



tt/t 



«'e?l 



— — ^ EiAf(q, ^ - fi o) [5 , fcu;0-a5 , fc+ 9 , tl ;+c] < 
m "0 ™ 



mftn 



£^jJfc+fW(q,ft-ft 



X [(fc + q)igku)&a9k+q,LO+tt-n 9k+q,u)+n + kig k ^9k,LO+tt &a9k+q,u 



2 /|2> 



+0^,^ 



(192) 



We consider a slow spin dynamics and assume that the spin gauge field has 
only a zero frequency component, i.e., A"(q, ft) = 5^ !0 A"(q). This is justified 
when ftr <C 1. We then obtain 
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u ai 

+^E^(* + f)^) 



Q/H 



X [{k + q)igkuj^a9k+q,u}9k+q,LO+n + fcj#fc,u;<7fc,u'+no C,r aS , fe+g,W+fio] 

+0(£ 2 ,A 2 ). (193) 

5.^ iSpm densities 



The spin density in the rotated frame is calculated by multiplying the lesser 
Green's function by a Pauli matrix and taking a trace. Only ± components 
are necessary, whose Fourier transforms are given as 



duj 



5j = -2i Urn / ^£tr[0r±G< fe+q (u,, w + fi )] 



? ±(0) , ~±(1) 
s q i °q l 



(194) 



where a x ±ia y = 2a±, s^ ' and s^ 1 ^ = s^ 1 ^ + s^ lb - ) are the equilibrium and 
current-driven parts, respectively. Each contribution is given as 



s±(o) 



4fc/|E&,W^h 






a °^°t^T fi o 






Bn9 



0fc-§,= F ,w-^0fc+§,±,a,-i£yfc+f,±,w+i£ 



5^ 



1< 



0fc_*.X.a,-2tt0» 



Sj< lb > = lim V ^ 
where 



"o-o^ fi 



< 



?fe- f,T,w- ^#fc+f,±,w+^ 



(195) 

(196) 



Bt(k,q)=2^J2J,(k)A±(q) 
These contributions are shown in Figs. 19 and 20. 



(197) 
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k-i 



Sf = at 




Fig. 19. Electron spin density without current at the lowest order in gauge field 
(represented by wavy lines). 





^-f-f^V 


~(1) 

s± = 






t+foH^ eX 




Fig. 20. Diagrams representing the current-driven part of the electron spin density at 
the linear order in both gauge field (represented by wavy lines) and applied electric 
field (dotted lines). The first two processes are self-energy corrections , s ^ \ and 
the last process is from the correction of the current vertex, s ( >. 

In this section, we will derive the leading contribution in the clean limit, r —* 
oo. (When spin relaxation is concerned, we need to estimate the next-order 
contribution of order of (e^r) -1 . See §14.) Using Eq. (188), the equilibrium 
part is obtained as 



s±(o) 



n C ci i j 



(19* 



By use of 



9ko(u) ^ fkJ(uj-e ka ), 



(199) 



we obtain the equilibrium contribution as 



*t (0, = ^E^(9) 

fc 



/fc+f ,± ~ /fc-§,T 



(200) 



The first contribution to current-driven part is obtained by use of Eqs. (180) (188) 



as 
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*t (la) = EE^ 



'(« + t)«m 

+f(u 



k + 



w-%9 



k+l±,u,-^ 



'fc+f.i.uM-^ 



'fe+f,±,' 



2 



S'fe-|, :F ,a)-^5'fe+f,±, W -^ 



^fe-f^-HaSfc+f.i^-Hp 



f,±,w- 






2 
2 



2 



'*-*,*,<"- 



.a^+fk^fe+S.i.w+Sl 



?fe_|, :F ,a)+^ 9 'fe+f,±, W +^ 



.fMfi 1 ! 



_• 



(201) 



The dominant term arises from terms containing both g r and g a , since the 
contribution from g r only is higher order of (e^r) -1 after summation over the 
wave vector (k), as is well-known in electron transport phenomena. This is 
easily understood from following example. In the case of (e^r) -1 <C 1, the 
summation over the wave vector can be replaced by an energy integral as 



\ roo 

TtE^L = / <M 
V u J-e F 



[e - cMY + (j- 



Va 



de- 



e'-h ^r 



2nu a T, 



(202) 



where v(e) is the energy- dependent density of states and v a = v[aM\ Within 
this approximation, the product of g T results in J2k(.9ka) 2 — v o I^°oo ^ e ( e ~~ 
2 = 0. We therefore see that the dominant term is simply proportional 



2t 



to / f uj + -yj - f f to - -f J = Q f(u) (Q — »• 0), as is familiar in low energy 
transport properties. Therefore, we obtain 



*T la) -EE^/V) 



ojk i 



* + I 



9k 



•f, T ,a;fl , fe+f,±,a;fl'fe+f,±, 1 



2 i k 2' 



T,a;fl , fc-f,=F,a;5 , fc+f,±,( 



The contribution s^( lb ) is of higher order in (e^r) -1 than s^ la \ since it has 
a smaller number of Green's functions [50]. Thus the current-induced part is 
given by s^ 1 ^ ~ s^ la \ Since we are interested in low temperatures, /c^T/e^ <C 
1, we can replace f'(uj) ~ —6(u) in Eq. (203), and hence 



s±(i) 



2n 



EE^f 



k i 



X 



k + 2 ) . 9k-^ T 9k+^,±9k+^± + [k ~ 2) . 9k-2, T 9k-z, T 9k+z,± 



(203) 



(204) 



63 



where g k 2 = #£,_ g uJ=Q . Using the identities 



9k,a9k,a — ^ T \9k,a 9 k, a) 



1 



5 , fc-§,= F 5 l fc+§,± 



9k-z^9k+z± 



e fc+ f 



e fc _a =F2M + 27 



r-(fl'fc-§, T - 9h+i,±) 



e k+ i -e k _i T%M 



(fl'fe-f ,t _ S'fe+f ,±> 



(205) 
(206) 

(207) 



where 7 = -, we obtain, in the limit of r — ► 00, 



s±(i) 






IX 



X 



niJ H I ( rS 



k i ""' fJ- 

q 



9k,± 9k,± 1 7 ( 1 , Q\ 9k,T 9k,^f 



2JJ e k+q -e k ±2M 



2/ e k+q -e k T2M 



-™q* 9l, T J v \ k + 7;) 5 ^ +q ~e k T 2M) - gl ± J„ [-[k + j-)) 5(e k+q -e k ± 2M) 



(208) 



The first line in square brackets arises from the real part of {e k+q — e k ± 2M + 
27) ~ x etc., while the second line arises from the imaginary part. 

The adiabatic limit can be obtained easily from the results Eqs. (200) (208) by 
setting e k+q ~ e k . The result is (s^ = (Sj<°) + s±M)|( ad )) 



~±(ad) 
b q 



MV 



E 



At(q)(f k+ -f k -)+i 



eEjAt(q )T 
dirm 2 



k (g T k + - 9l- 



± (sAt(q) + f ■ AHq) 



(209) 



Here d is the dimensionality, 



s = n_i_ — n_ 



j s = — (n+ - n_)e 2 rE = j + - j. 



m 



(210) 



are the electron spin density and spin current density (both defined without 

spin len; 

density. 



1 k^ 

spin length of ~), respectively, where n a = -^ is the spin-resolved electron 



Going back to the general case, spin polarization in real space, s ± (x), is ob- 
tained by Fourier transform as s ± (x) = J2 q e ~ iqx Sq- After some calculation, 
the equilibrium contributions, s e (x) and s> (x), are obtained as 
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s^ix) 



8®(X) 



M- 



T?£ e ~*" ((«* x e z ) x ■ A (q)) X r(q) + Mx) ■ A (q)) X T(q) 



,(°). 



M- 



YjT*** Mx) ■ A (q)) x? } (q) - ((e, X e z ) x • A (q)) X f{q 



(o), 



(211) 



where we used \ £ ± e^^A^(q) = (e (x)xe z )-A M (q) and \ E±{Ti)e Ti(p{x) A±{q) 
e^(cc) ■ A^(q). The correlation functions are given by 



,(°) ( 



fi 



k± 



2_hM 
V tl ±2M I 2k -i+i 2 



xi 0) (<?) 



k± 

2HM ^7T 



2d) 



£-(/ fe+ -/ fe -)<5(±2M 



k± 



2k ■ q + q 7 
2m 



f212l 



Similarly, the current- induced contribution is calculated as 



sP(x) 



s { l ] (x) 



M 

EjEi 
M 



Q 



-tqx 



-iqx 



((c, x e 2 ), • A 4 (g)) x { i\q) + M*) ■ A,(q)) x^iq) 



(213) 



where 



x^fa) 



X? ] (9) 



erM 



EM* 



9k± 9k± 



3irm 2 Vj^V~ V" ' 2;;^ fc+q -e fc ±2M 
erM 



£ (±?H U ■ (k + f)) 5( efc+q - e fe ± 2M)ife - <^ ± ) 



(214) 



37rm 2 1/-^ V 2 7 V" V" 



Integration over k in Eq. (212) is carried out to obtain 



X?\q) = sx?\q) 

~(o)/ n _ 3 1 ly 

Xi w -4(3 + C 2 )5 2 Tl 



■^(g 2 -l)(g 2 -C 2 )ln 

2g 



+ i)(?±0 



(?-i)(gTC) 



±(i±C)(C±<f) , 

(215) 



where g = \q\/(2kp), s = n + — n_ = f^2"C(3 + C 2 ) is the electron spin density, 
k F = \{k F+ + k F _), ( 



1 + 0(g 2 ). Similarly, we obtain 



,-,(o) 



and xi (?) is normalized to be xi (q) 



r,(o). 
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xf\q) 



k F 3 (l 



-Y 

16tt q^ 



i+C ~~ ~_ 
dkk9{kq 

i-C 



? 2 ±C)I) 



k F *$\ 



Q), 



(216) 



where the integral is to be taken only in the regime k > \ ^(q 2 ± () |. Correlation 

functions, Xi an d X2 \ are plotted in momentum space and real space in Fig. 
21. 
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Fig. 21. Plot of correlation functions, \\ (q) an d X2 w)> an d their Fourier trans- 



.--,(0), 



,~,(°) 



forms, x\ ( x ) an d X2 ( x )i describing the nonlocal component of spin density (and 
torque) in the absence of current. x\ has a finite adiabatic component (q = 0, i.e., 
local component), while x 2 does not. 



The summation over k in Eq. (214) is carried out using ~ J2k = I dev(e) f_ a dc ° s k . 
where 8k is the angle of k measured from the q direction. We carry out the en- 
ergy integration first, assuming that the dominant pole arises from the Green's 

function and the residue contribution from l ,„,. is neglected. After some 

calculation, we obtain 



Xi («) = — (n+-n-)xi(q) 
m 



*f> 



1 



3 + e? (1±c) 



1 H ; r^ In 

2(l±C)g 



q)(q±0 



q)(q tC) 



,(217) 



r,W 



where n± = -g^f- is the spin-resolved electron density, xi is normalized to be 



?fl)t 



M) 



Xi (q) — 1 + 0(q ), and thus xi (q)E = -xi (q)j B - We similarly obtain 



±<-M)(z\ 



7&%) 



m 2 M 2 a 9 st (q) __ a n 



Qimkp e q 



xPb), 



(218) 
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where 



and 



x? } (g) 



2 1 + 3C 2 g 



M?) = 



1 (&f+ — fci?_ < \q\ < kF+ 
otherwise 



(219) 



(220) 



represents the regime of Stoner excitation, o"o = e 2 nr/m is the Boltzmann 
conductivity, and n = n + + n_ is the total electron density. As is obvious, for 
small q, X2 = 0. Correlation functions, Xi an d X2 > are plotted in momentum 
space and real space in Fig. 22. Note that the Xi~ an d X2-terms are proportional 
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Fig. 22. Plot of correlation functions xi (<?)> X2 (?) an d their Fourier transforms, 
%i (x) and x\ ( x )> describing nonlocal components of spin density and torque in 
the presence of current. x\ has a finite adiabatic component (q = i.e., local 
component) while x% does not. 

to spin current and charge current, respectively, only in the adiabatic limit 
(q = 0), but are not necessarily so when nonadiabaticity sets in, since x~i(q) 
and x.2{q) can depend on polarization ( in a complicated manner. 

In the real space representation, the spin densities are given as 



sf\x) = 


1 

~M 


s ( :\ X )= 


1 



and 



d 3 x 
d 3 x 



sx ( i\x - x') ((e^ x e^ • A (x')) + fc/xf } (* - x') (e^x) ■ A {x')) 
sx { i\x - x') (e^x) ■ A (x')) - k F 3 xP(x - x') ((e x e z ) x ■ A (x')) 



:221] 
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sP(x) 



sPix) 



J2 / d 3 x' f s xi(x - x') ((e^ x e z ) x ■ A^x')) 



eM 
+M\x - x')(e^(x) ■ Aitf)) 



■^Hfd 3 *' [f s X { i\x - x')(e+(x) ■ Mx')) 



-j'x-Pi* ~ x ')(( e <t> x e z) x ■ Ai( x ')) 



(222) 



The electron spin density induced around domain wall is summarized in Fij 
23 




Fig. 23. Spin polarization of conduction electrons (denoted by small arrows) around 
a domain wall (large arrows), and its effect on the domain- wall dynamics. Top: 
The adiabatic component, s z , along the localized spin does not affect the dynamics. 
Middle: Under current, the spin-transfer torque is induced by the polarization s^ 
out of the wall plane, which induces rotation of localized spins within the plane, 
resulting in translational motion, X. Bottom: cfto is induced by the component sq 
in the ^-direction in the wall plane. Such polarization arises from spin relaxation 
processes [47]. Non-adiabaticity also induces 4>o [26] but the spin density profile is 
nonlocal and is different from the one due to spin relaxation. 



68 



9 Torque and force from electron 



9. 1 Torque 



The torque acting on a localized spin at x is given by 



r(x) = M (s<f,e0 - sge^) 

= - I d 3 x [sxi\x - x) [e g (x) (e^(aj) ■ A (x')) - e^(x) ((e^ x e z ) x ■ A (x'))] 
-kF 3 X 2 ( x ~ x ') i e e(x) ((e x e z ) x ■ A (x')) + e^x) (e<f,(x) • A (x'))] 
+ Y,-Xi\x - x') [e e (x)(e 4) (x) ■ A^x')) - e^(a?)((e^ x e z ) x ■ A^x'))] 



Yl —X2 ] ( x - X 1 ) [e 9 (a;)((e ? i x e z ) x ■ Ai{x')) + e (j} (x)(e (t> (x) ■ Ai(x'))} 



This torque is the sum of an adiabatic part and a non-adiabatic part, 

T (x) = T ad (x) + T na (x), 



(224) 



where the adiabatic part r ad is a purely local part, given later in Eq. (230), 
and the nonlocal part is defined as 



T n Jx) = T(x)-T ad (x). 



(225) 



(223) 



It is in general non-local due to the non-adiabatic correction. The total torque 
on the whole structure is given by its integral, r = J c1 3 xt(x). 

In the case of uniform polarization, 4>{x) = <po, Eq. (223) becomes somewhat 
simpler using (e x e z ) x ■ A^x') = (e x e z ) x , ■ A^x') = Aj(aj'): 



T(X 



(x)=cf>o 

3s ~(i) 



drx 



,~,(°) 



sx\ (x — x') (n(x) x A (x')) 



J2^ X r(x-x')(n(x)xA i (x')) 



3^(°) 



+k F i X2 (x - x') [e e (x) {e e {x) ■ A {x')) + e^ (e^ • A {x'))] 



J2 —xf(x - X 1 ) [e g (x) (e g (x) ■ Ai(x')) + e^ (e • A* (a/))] 



3 ~(i) 
e 

where we used n = eg x e^. The ^-component is given as 



(226) 
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r z (x 



{x)=4>o 



sm6 (x) fd 3 x sx ( i\x - x')At(x') + ]T ~Xi\x - x')Af(x') 



n L 

3 ,~,(1), 



VxH* - *') (e 9 (x) ■ A (x')) + £ -A } (x - x') (e e (x) ■ A^x')) 

a e 



(227) 



9.2 Adiabatic limit 



The expression for the spin density takes the form of an integral such that the 
spin density at a given point is determined by all other points in space. The 
adiabatic limit is the only exception. In this limit, the momentum transfer q 
from the gauge field is negligibly small compared with the electron momentum 
k, and so we can estimate all correlation functions Xin (o) at q = 0. Thus, 
X ( i\q) = #(<?) = 1 in Eqs. (211)(213), and #(q) = #(q) = 0. Namely, 



,~,(°)( 



,-5W 



Xi'(x-x') = x\ > (x-x') = 5(x-x') 



,~,(°) 



X% '{x - x) = x\ >{x - x) = 0. 



(228) 



We therefore obtain, by noting ((e^ x e z ) x ■ A M (aj)) = (eo(x) ■ A^x)), 



4 ad) N = -T7 



sPUx) 



1 
M 

1 
M 



sA e (x) + ^A°(x) 



sAt(x) 



3 -±A+(x) 



(229) 



The local torque in the adiabatic limit, given using eq. (173) as T ad = M J d 3 x (s^ eg — sf e^ 
is then given as 



T ad (x) = nx (sAq + ^A; 



sd n 



1 



js ■ V TO 



(230) 



The second term is a spin transfer torque, and the first term represents the 
renormalization of the magnitude of the spin. 
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9. 3 Force 



Let us look into the force exerted by the electron. The total force is given as 
a sum of equilibrium and current- driven parts, Fi = F i + F^ , arising from 
s 1 ^ ) and s^ 1 -*, respectively. We first look into the force from the current, F^ 1 '. 
Using Eq. (208) and noting g r ka ~ — mb~(e ka ) on fc-integration, we obtain 



eE k rA 

± kq jk 



f = -EEE^(±«% ± : 



* (* + I), e^-T»±2A (4(«)^(" 9 ) - A 1{q) AH- q : 



-*ft (*+§). *(«** - c* ± 2A) (j4±(«)j4f (-«) + AJ(q)Af(-q 

= Ff + 5F?\ (231) 

where the first and second terms are defined as 



^eEkTA 
m 2 V 



± kq jk 

xnq k (k + |) 5(e fe+g - e fe ± 2A) (Af(q)A?(-q) + Aj(g)i4f(-« 
eE h rA 



^ 1} = EEE^I^(±)^ 



1 

2/ j Cfc+g _ e fc ± 2 A 

(232) 



"*> * + % , r- _ ... + , A (4(gW(-g) - Aj W ^(-,)) . 



The first term F ref vanishes in the adiabatic limit (g = 0) while the second 
term SF^ 1 ' remains finite. We now demonstrate that the term F rei represents 
the reflection of the electron due to the spin structure. In fact, we find, by 
using e k+q - e k = ±q ■ (k + § ), that F rcf is 



AttpF tA 2 

Ff = ~ E E E 2y S(e k± )5(e k+q -e h ± 2A)Af(q)A?(-q). 

± kq j 

(233) 

This force is proportional to the resistivity due to the spin structure, which 
will be discussed in §15.4 on the basis of the Mori formula. The resistivity is 
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calculated as (in the case of current along ^-direction) 

AttA 2 1 
P* = "^^TfE \A° z (q)\ 2 6(e k+q ^ - e^)^). (234) 

kqa 

We therefore see the relation between force and resistivity 

F rcl = 6 —p s n 2 V = eN cPs j, (235) 

m 

where N e = nV is the total electron number. This is the result in Ref. [26] 
extended to a general spin structure. 

The resistance due to spin structure, R$ = j^p s , is related to the reflection 
probability R, according to the four-terminal Landauer-Biittiker formula [107], 
as Rs = ^i^R' anc ^ nence Eq. (235) indeed relates the force to the reflection 
of the electron. Equation (235) can also be written, using the density of states 
N(0) = r $$%- (neglecting spin splitting), as 

F Tei = ^eV s N(0)2k F h^, (236) 

where Vs = R$I is a voltage drop due to spin structure and I = Aj is the 
current (A being the cross section). This equation clearly indicates that the 
force is due to the momentum transfer (2kph per electron, with frequency 
^y) multiplied by the number of electrons that contribute to the resistance 

CleV s N(0)). 

What is the origin of the other term SF^ 1 ' (Eq. (232)) ? Its meaning becomes 
clear in the adiabatic limit. In fact, in this limit, 5F} reduces to _p Ha11 where 



^ Hall = EEE^^(^ 2 (Af(q)A?(-q)+J$(q)At(-q] 

± kq j 



= ^E- [d 3 x(A+( X )A-(x)-A-(x)At(x)) . (237) 

Using 

(Af(x)Aj(x) - Ar(x)A}(xj) = -%\ sin e(d i ed j </> - djOd^), (238) 



and 

we see that 

uiHall 



d t S x djS = nS 2 sin Q{d i Qd j 4> - d^d^), (239) 



L.Y, 3 -f J d 3 xS ■ &S x djS) = -2ir-£ J -f®ij, (240) 



253 
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where 

$y = ^J d3xS ■ ( d i S x d 3 S ) ( 241 ) 

is a vortex number in a plane perpendicular to % and j directions. In the case 
of a thin system (with thickness d), this reduces to 

$xv = n v d, (242) 

where 

n v = ^Jd 2 xS-(d x S x 9,5), (243) 

is a topological number in two dimensions. This force is, in fact, a back reaction 
of the Hall effect due to spin chirality [108-110], and was derived by Thiele 
[111] and by Berger [35], then rigorously in Ref. [112], and also assuming a 
vortex structure in Ref. [113]. The Hall effect due to vorticity arises from the 
spin Berry phase [108] or spin chirality [109]. For a domain wall, the vorticity 
is zero and so _F Ha11 = 0. Defining the gyrovector of the vortex, G, as 

2-kS 
G, = -^e ijk $ jkl (244) 

we can write the Hall force as 

Pn 3 

fM = = "M (ixG) ' (245) 

where Pj represents direction of spin current flow. 

There is another force arising from spin relaxation, the /3 ST term. Comparing 
the definition of the (3 term, Eq. (14), with Eqs. (130)(129), we see that the (3 
term corresponds to spin polarization of 

s » = "2^7« ' v) "' < 246) 

This spin polarization exerts a force, according to Eq. (136), of [48,47,112,50] 

F p = ^^srPj, (247) 

where fi = J d 3 x (V z n) 2 . For a planar wall, 

ft = ^N w . (248) 

These three forces are schematically illustrated in Fig. 24. 
One may think this strange, but there is a force arising from the equilibrium 
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Hall 




Fig. 24. Schematic summary of three forces acting on spin structures. Electron 
reflection pushes the structure along the current flow (F rci ), and the Hall effect 
due to spin chirality results in a force in the perpendicular direction (F ). For a 
domain wall, ,F Ha11 = 0. Spin relaxation results in a force along the direction of the 
current (F@). 

component of spin polarization. This force, given by (see Eq. (172)) 

FP=2M [d*x(-sPA! + S fAl), (249) 



with s e and s^ given by Eq. (221), arises when the spin structure is moving, 
and represents the effect of relative velocity between the spin structure and 
the current- carrying electron. Let us consider again the adiabatic limit for 
simplicity. Using Eq. (229), the force in this limit is given by 



Fl CR = F/ 0) | adiabatic = 2s I d 3 xn ■ ( A x A t ) 
s 
~2 



f d 3 xn ■ (d n x d t n) = (* )«, (250) 



where the vector $ is defined as (3>o)i — &oi of Eq. (241). This equilibrium 
contribution to the adiabatic force thus has again a topological meaning but 
in two dimensions including time and space. The total force due to current is 
summarized as 

F = (eN e p B + P^Psr) 3 + P^j x G + 2tts$ + SF, (251) 

where 

5F = SF® - F Ha11 + F (0) - F ren , (252) 

is an additional non-adiabatic correction to the force (besides F ref ). In the 
adiabatic limit it reduces to 



P> fin 

F (ad) = p M + p • xG + 2irhs $ . (253) 

2eb 2eb 
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10 Effective Lagrangian (adiabatic limit) 

10.1 Spin transfer torque and Berry phase 



In the gauge transformed frame, the effective Lagrangian of the localized spin 
is obtained simply by taking the expectation value of Ha, Eq. (164) [114]. In 
the adiabatic limit, this value is given by 



a\o- a a k 



(H A ) = hY / (A a M + -k.A«(q) 

hqa \ 

= J d 3 x (hsA z (x) + —j ■ A z {x)\ , (254) 



(only z spin components are finite at the lowest order in A). Therefore the 
effective spin Lagrangian reads (neglecting dissipation) 



Ls° = L s — (Ha) 
d 3 x 



h - ,3 



S tot 0(cos£ - 1) + ^-[{j ■ V)0](cos£ - 1) 



H s , (255) 



where S to t = S + ^- is the total spin including electron spin polarization. 
This Lagrangian for the adiabatic case can be written as 



where 



£s eff = / ^ fiStot[(dt + ^s ■ V)0](cos^ - 1)] - H a , (256) 

Pa 3 , , 

- S ^ J - (257) 

is the drift velocity of the electron spin. The solution has the form f(x — v s t) 
(f is any function), which clearly indicates that the adiabatic spin transfer 
torque induces a stream motion of localized spin structure with velocity v c . 
This applies to any spin structure. We can also say from Eq. (256) that the 
spin transfer torque induces a spatial spin Berry phase. 



10.2 Nucleation of domain wall by spin current 



As seen from Eq. (256), the spin current favors a magnetic configuration with 
spatial gradient, or more precisely, with finite Berry-phase curvature. It is thus 
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expected that a large spin current destabilizes a uniform ferromagnetic state. 
This is indeed seen from the spin-wave energy around a uniform ferromagnetic 

state [114-116], 



fijf) = ^ v /(/ e 2 A 2 + ^(^2 + i + K ) + k . Vg _ (258) 

The first term is the well-known spin-wave dispersion with an anisotropy gap. 
The effect of spin current appears in the second term as the Doppler shift 
due to the drift velocity of electron spins [115]. For sufficiently large v s , Q k 
becomes negative for a range of k. This means that there exist states with 
negative excitation energy, indicating the instability of the assumed uniformly- 
magnetized state [114-116]. The wavelength of the unstable mode around the 
uniform magnetization is k = —(K(K + K±)/J 2 ) 1 ^ 4: [116]. The critical current, 
j nuc , is given by [115,116] 



2eS 



2 



Jnuc = -f-T^d + Vl + K). (259) 

The true ground state under a large spin current was found to be a multi- 
domain state [117]. Namely, the energy of a domain wall becomes lower than 
the uniform ferromagnetic state when a spin current exceeds a critical value 
Jnuc- Nucleated domains flow at a velocity v s if K± <C K. 

When K <C K±, as is usually the case of a film or wide strip, vortex nucleation 
occurs at a lower current than j nuc , as was shown in Ref. [118]. 



11 Torque and force on domain wall 

In this section, the torque and force on a rigid domain wall are estimated. 
The gauge fields for a domain wall is given by 



A e = -- sin 6 (f) 

^o = 2y sin ^o 

Ai = -8 itZ ^sme , (260) 

where sin# — [cosh x ~ x }~ 1 - Noting that sm6o(x) is odd in z — X and 
/ d 3 xj(2 (x — x') = X2 (q — 0) = 0, we see that the total torque on a wall is 
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obtained from Eq. (227) as 



' \\ 



— / d 3 x / d 3 x sin 6 (x) sin 9 Q (x') 



s- 



*n* - ^) - £ ^-*n* - *o 



A" 1 v ' VeA 



< 



(261) 



This is the full expression of the torque acting on the wall including non- 
adiabaticity (but without spin relaxation). Separating the adiabatic contribu- 
tion, it can be written as 

Tw = fir G? " sji ) + 5t ' (262) 

where 5t represents non-adiabatic contributions. 

For domain walls, the adiabatic Hall force term _p Ha11 vanishes since the domain 
wall has no vortex charge. The total force on the domain wall is given by 



F w = F rci + F ren + F? + 5F 

= eN cPw3 + ^ S O + F? + 8F, (263) 

where p w is the resistivity due to the domain wall. 

From the above argument, the equation of motion of a rigid wall is therefore 
obtained as (Eqs. (134) (135)) 



0o+a T = ^vb (Fw+Fpin) 

X-«A0 o = ^sin20 o + ^r w , (264) 

or 



X sa 3 ■ 

(j) + a— = / ref + f + /pi,, - —0 O + Sf 

X - aX^ = v c sin 2^ + f g (pi - sx) + ^fr. (265) 

where /( ref >P 1IX >0) = ^ fM.p™.^). We see that the contribution to the O term 
from F ren represents the dressing effect of localized spin due to conduction 
electron polarization. Namely, the spin constituting the domain wall is now 
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sa 3 



given by the total spin 

Stot = s - 

and the equations become 



.X A 



StoJo + aS- = — (F ref + FP + SF + F pin ) 



K A Po 3 A 

S tot X - aSA0 o = -t-S 2 sin 20 o + — j + tt^t. (267) 

2ft 2e n,iv w 

The spin renormalization factor appears also in the equation for X, and this 
contribution guarantees correctly that the spin-transfer effect arises propor- 
tional to the relative velocity (^j — sX) between the wall and conduction 
electron spin. These terms were considered phenomenologically in Refs. [35,24] 
and [54], but these effects naturally appear in our formulation. 



12 Domain wall dynamics 



In the following calculation, we simply write Stot as S. In non-adiabatic contri- 
butions, we include F ref , which is qualitatively important in dynamics, while 
we neglect SF and St, since these affect the dynamics only quantitatively. The 
equation of motion of the wall we consider is given by 



X _a? 

T~2^SA' 



^0 + a— = TT-TTV^wj + /pin 



a 3 



X - «A0 O = v c sin 20o + ^Pj, (268) 



where 



e 2 



(3 W = -nAXR w + pp BI , (269) 

n 

represents the total force acting on the wall (R w = p w L/A is the wall resis- 
tance), 

v c = -^-, (270) 

and 
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f pin = -—£X6(Z-\X\). (271) 



The coefficient /3 W is the effective j3 that acts on a rigid wall. In other words, 
a rigid wall feels electron spin relaxation and non-adiabaticity (reflection) 
exactly the same. The equation of motion in terms of dimensionless parameters 
is written as 



di (X — a0o) = sin 20 o + Pj 

d~ t U + aJt) = -V Xd(£/X -\X\) + M, (272) 



where i = tv c /X, X = X/X, Q = Q pin X/v c = 2 -f\^ P = P,J = ^fcj, 

I02 _ \/2J? - AJiL^^ or^ T/_ = A^02 t2 _ ^^iE/£ x2 



T/ — 1Q2 — T/ ^ A — 4 _Kq_ / A \ J if _ Mw Q2 t2 _ ft" ""p 



4 u c A ' 



The equation of motion can be written as 



d~ t X = —^— (sin20 o + {P + ap w )j - aV X9(£/X - \X\)) (273) 

1 + or K J 

d i^ = T^—2 (O 8 * - Pa ^ ~ a sin 2 ^ ~ VoXe^/X - \X\)) . (274) 

We will solve this equation of motion below for the case of steady current. A 
moving domain wall, free from the extrinsic pinning effect, is described by the 
above equation with Vq = 0. There is an analytical solution in this case as we 
see just below. 



12.1 Solution without extrinsic pinning 

Let us consider first the case without pinning potential of extrinsic origin. In 
this case the analytical solution is easily obtained. In fact Eq. (274) has a form 
of 

d^o = B - Asm2(j) , (275) 

where 



A- a 



2 



1 + a 
B J? Pa)~j 

1 + a z 
The solution to this differential equation is given by 

tan0 o = — (A - u cot(u(t - t ))) , (277) 

B 
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where cot 9 = ^ 



LO 



= VB 2 - A 2 



a 



l + a 2 \ 



4- -1, 



(278) 



and 



P- 



(279) 



In the case of imaginary uj (i.e., A 2 > B 2 ), Eq. (277) is still applicable by 
replacing uj — >• «|u;| and tan(u;(t — to)) — > *tanh(|u;|(t — to))- We require the 
initial condition 0o(O) = at t = 0. Then we find the solution as 



B tan ut 

tan 0o = , 

uj + A tan out 



(280) 



and thus sin20o, which appears in the wall velocity, is obtained as 



sin 20o : 



A + B sm(2ut - ■&) 
B + A sin(2u;t - ■&) 



_J_ 



-Pb' 



/?« 



^-^i 



l 



- \2 
J- 
ja 



^ - PJ j + sin(2cjt - •&) 
- 1 



i - sin(2c;t - 0) 



(281) 



where sin^ = -^ = ^ w * - - , cos?? = ^. We see that an anomaly appears when 



ui switches from real to imaginary, i.e., at \j\ = |j a |. Above and below |j a |, 
the wall dynamics is quite different. For \j\ > |j a |, the wall velocity (273) has 
an oscillating component. This oscillating dynamics is similar to the behavior 
in a high magnetic field, known as Walker's breakdown [97]. In contrast, for 
lil < |ja|, 0o reaches its stable angle (given by Eq. (285)), and simple sliding 
motion of wall with constant velocity (Eq. (286)) is realized. 

The average velocity for \j\ > \j a \ is obtained by use of Eqs. (273) (281) and 



dt 



C 



TJo C-sin(274) VC^^llC 



(282) 



(for \C\ > 1) to be (T = *) 



SO 



x )= ~rJ + i , „ 2 s g n 



(3 W , , 1 
a 1 + a< 

a 1 + a z j 



\ 



3_ 

yjlX 



2 



1 



-J + T—iHf-Jl (283) 



When j ^> |j a |, the velocity becomes 



P + af3 w ~. 
l + a 2 



x)->VrxJ S>U.D- (284) 



Below j a , a; is pure imaginary and sin20 o approaches at t — » oo a steady value 
of 

Sin20 o ^!=(^-Pb, (285) 

and so the terminal velocity is obtained as 

x) = —j. (286) 

1 a 



This speed can be larger than the spin-transfer limit, (X) = Pj, if ^ s - 
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large, but this high speed is realized only at a small current, j < |j a |, i.e., 

A special case arise when — = P, where simple sliding with (fi = and 
X = Pj is realized [54]. 

These behaviors of wall speed are seen in Fig. 25 (in the presence of extrinsic 
pinning potential). The most significant effect of /? w is to shift j a . When /3 W < 
aP, j a remains close to the intrinsic threshold, but when f3 w exceeds aP, j a 
becomes negative, meaning that no anomaly appears in the velocity curve. 



12.2 Pure spin-transfer case f/3 w = 0) 



In this subsection, we look into the case of /3 W = 0, i.e., the dynamics driven 
by purely spin torque. As we see, a significant feature of the domain wall, 
intrinsic pinning, appears in this regime. 
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Fig. 25. Wall velocity as function of current for £1 = 0.5 and a = 0.01. A jump 
in wall velocity is seen at threshold, j = j c . A crossover from I-a)(j c ~ Cl//3 W ) to 
I-b)(j c ~ Ct) regime is seen at (3 W ~ -g ~ 0.1. Intrinsic pinning occurs only for 
/3 W = in the present case of Q < 1 (see Table 2). Linear behavior of velocity seen 
for large /? w (> 0.1) is the perfect adiabatic limit realized for j 3> j a (Eq. (284)). 

12.2.1 Intrinsic pinning 

In the case of /3 W = and Vq = 0, the wall velocity becomes 



X 



l+a J 



f - Qlf 



Q<~3l) 

Q > f c ) 



^287) 



and f c = -i becomes a threshold current density [26] (Fig. 26). This effect, 
which could be called an intrinsic pinning, arises from the fact that a domain 
wall is a collective object that can deform (i.e., can develop 0o) and absorb 
spin torque. In dimensioned quantities, the intrinsic threshold is at 



P Q2 

& = *KP K * 



(288) 



The intrinsic pinning effect is robust against extrinsic pinning [26]. In fact, 
the dynamics near threshold is studied based on the equation of 0o derived by 
eliminating X from Eq. (272); 



(1 + a 2 )<9|0o + ad i( j) [2 cos 20 o + V J + V sin 20 o + jV P = 0. (289) 

This equation indicates that 0o is like a classical particle moving in a potential 

v n 



v A 



sin 20 o + jP(f> 



(290) 



1 + a 2 V2 
From Eq. (289), we see that the energy barrier for 0o vanishes when 

j c ~ P~\ (291) 
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Fig. 26. Domain wall velocity in the case of (3 W = with not too strong extrinsic 
pinning, Vq < a -1 . A finite threshold appears due to the intrinsic pinning due to 
perpendicular magnetic anisotropy energy. 

irrespective of extrinsic pinning strength. Once 0o escapes from a local mini- 
mum, its average velocity is given by Eq. (289) as 



dm - — • 

a 



(292) 



This velocity corresponds, via Eq. (274), to a maxmum displacement of the 
wall inside the pinning potential, 



X n 






aV 



Since a is small, \X- 

depinned as soon as 

even in the presence of extrinsic pinning if 



-i> 



aax i easily exceeds £ (unless Vq > a 
0o runs. Thus the threshold current is given by j 



(293) 
and the wall is 



Vo< 



o 



(294) 



Let us note that the domain wall moves even below the intrinsic threshold. 
Actually, the shift of the domain wall can be quite large. Here we consider no 
extrinsic pinning and /3 W = 0. For j < j a (j a coinsides with intrinsic threshold), 
0o reaches the value given by Eq. (285). Equation (272) then indicates that 
the shift of 0o is associated with the shift of the wall itself as 



AX 



X_ 

2a 



sm 



1 Pi 



(295) 



Therefore, even for a current 2% of the intrinsic threshold, the domain wall 
can move over a distance A (if a = 0.01). Such motion at very low current 
would be enough for applications. For this, a very clean sample is required, 
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Fig. 27. Threshold current j c plotted as function of pinning frequency fi for /3 W = 
with a = 0.01 and P = 1. We see that intrinsic pinning persists if £7 < \/^- 

since even weak extrinsic pinning can result in creep motion in the low current 
regime [88]. 



12.2.2 Extremely strong pinning case 

When Vq > -, the threshold current is governed by the extrinsic pinning. In 
fact, Eqs. (273) (274) in this strong-pinned case (and /3 W = 0) reduce to (within 
the pinning potential {\X\ < £) and at large current) 



d~ t X~ 



1 



0, 



l + a 2 
1 



Wo 



l + a s 



Pj - aV X 



:VnX. 



(296) 



We thus see that stream motion of the wall occurs if j c > j%, where 



(297) 



The behavior of the numerically determined j c in the case of /3 W = is plotted 
in Fig. 27. 



12.3 Effect of force fl3 w ) 



Now we look into the case of non- vanishing (3 W . We will see that there are three 
different behaviors depending on the extrinsic pinning strength. The analysis 
here is based on Ref. [51]. 
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12.3.1 Weak extrinsic pinning regime : I) 

Under weak current, j < 1, 0o remains small and the wall dynamics is well 
described by X only. This is the first regime I). Linearizing the sine term in 
Eq. (272) as sin20o — 20o 5 <Po can be eliminated to obtain a simple equation 
for X as [20,119] 

(1 + a 2 )d 2 X + -diX + Q 2 X = F, (298) 

where 1/f = 2cm 1 + |VoJ and 

F = 2/3J, (299) 

is a dimensionless force due to current. We consider a case of steady current 
and weak damping; 2f2f > 1 A general solution to Eq. (273) (inside the pinning 
potential) is given as 

X(t) = ^ + e~^ (Acos&i + B sm&i), (300) 

n 2 



where ft' = JVt 2 — -^ and A, B are constants. The initial condition required 

is X(0) = and dfX(0) = Pj. The second condition on the wall speed comes 
from the first equation in Eq. (272) (with 0o(O) = and a ^C 1), and is the 
most important consequence of the spin-transfer torque; namely, spin-transfer 
torque gives initial speed to the wall. With these initial conditions, we obtain 
a solution as 

X(t) = ^&f (l - e-* (cosQ't + ^sinn'i^ + ^e"^ smQ'i. (301) 

The first part is governed by a force from /3 W and the second is driven by a 
spin-transfer torque term. In the case of small /3 W , the spin torque contribution 
leads to a maximum displacement 

|A max | ~ |j, (302) 

while 

|X max | ~ *£z~j, (303) 

if (3 W is large. (We assumed here that damping is weak (Q!t ^> 1).) The 
first regime corresponds to regime I-a) and the second to I-b). The threshold 
current in each case is given as 

Jl a) - |f, (304) 
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and 



& ~ M- (305) 



The crossover between regimes la) and lb) occurs at 



/V ^ jCt. (306) 



In terms of dimensioned quantities, 



and 



Jc ~ Pa" h A ~ S V^, Jc ' [ ° 



. Ib) = _^X_^ = SeVoJ_)^_ = P^V^X, 
Jc " Av c a" \(3X ~ ha" \(3 W \ £ ' " S |/3 W | K ± ^ { ] 



Note that simple comparison of the pinning force and F in Eq. (298) gives a 

1 SeVg 1 A 2 

2 hcfi |Av| 5 



result correct up to a numerical factor, j c lh > ~ I^Mi^_ *_ 



The pinning strength Vq is experimentally accessible by driving the wall by 
a magnetic field. A magnetic field B along the easy axis adds a term in Eq. 
(265) 

fs = 9J ^B, (309) 

where g = 2. By a simple comparison of pinning force and magnetic field, Vq 
is written in terms of the depinning magnetic field B c as 

V = |^bS c |, (310) 

and so j c Ib ^ is simplified to be 



12.3.2 Intermediate regime (intrinsic pinning) : II 

When the extrinsic pinning is not weak, the threshold current becomes j c > 
0(1). For such a case, 0o no longer remains small, and the wall dynamics is 
governed by this angle variable. In fact, the effective mass of a </>o-" particle" is 
given by 1/Vo [31] (see Eq. (312)), and it becomes lighter than the correspond- 
ing mass of an X-" particle" given by 1/K± when extrinsic pinning is strong, 
Q > 1. By eliminating X from Eqs. (272), we obtain 

(1 + a 2 )df 0o + adrfo (2 cos 20 o + V ) + V sin 20 o + jV P = 0. (312) 
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Thus /3 W does not affect the dynamics of 0o- (Correctly, this feature is specific 
to a harmonic pinning potential, and unharmonicity results in the appearance 
of /3 W . In fact, the j3 term is eliminated from the equation of motion if one 
replaces X in Eq. (272) by X' = X - %^-j (i.e., shift of stable point of X). 
Even in the unharmonic case, nevertheless, we have numerically checked that 
the /3 W does not lead to an important contribution in this regime.) 

From Eqs. (312) (293), the threshold is roughly given by j c ~ P~ l , and is 
actually found numerically as 



rsj 



0.7 xp- 1 . (313) 



This story is the same as the case of /3 W = 0, Eq. (291), but with a different 
numerical factor (0.7 here instead of 1 for f3 w = 0). This difference comes from 
a different definition of threshold current with and without /3 W . In the analysis 
of the j3 w = case, even if X escapes from the pinning center at current 
j > 0.7P -1 , the terminal velocity vanishes if j < P~ l , since the motion 
stops due to the intrinsic pinning effect (i.e., 0o reaches a steady value and X 
becomes zero). On the other hand, if /3 W ^ 0, steady motion of X is possible 
as soon as X escapes from the pinning. This is the reason the threshold value 
in the intermediate regime is different for /3 W = and /5 W ^ (Fig. 28). If 
/9 W < 0, or more precisely, if the relative sign between /3 W and P in Eq.(272) is 
negative, the /3 w -term will drive the depinned wall back to the pinning center, 
and the threshold in this regime is given by the intrinsic value j c = 1. 



12.3.3 Strong pinning regime : III 

Equation (293) indicates that for extremely strong pinning, Vq > a" 1 , the wall 
is not always depinned even after O escapes from the potential minimum. 
Depinning occurs at 

Jc ~ -j- j, (314) 

as in the case of (3 W = 0. 

Figure 29 shows a numerically determined threshold. It is clearly seen in Fig. 
28(a) that behaviors for j3 w = and j3 w ^ are quite different except for the 
extremely strong pinning regime (Vq > 100 ~ I /a). 

Results of threshold current is summarize in Table 2. It is interesting that 
such a simple set of equation of motion results in so rich behaviors. 





(a) 



Fig. 28. (a): Threshold current j c plotted as function of pinning frequency fl = \J1Vq 
for several values of (3 W with a = 0.01 and P = 1. (b): Threshold current in weak 
pinning regime. Fitted curves for [3 W > 0.1 are j c oc fl 2 . For small /3 W (< 0.02), j c is 
linear in Q. 
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Fig. 29. (a): Threshold current j c in weak pinning regime plotted as function of 
pinning strength Vq for several values of /3 W with a = 0.01 and P = 1. For j3 w > 0.1, 
j c ~ 2/3 i s linear in Vq. (b): Weaker pinning regime. For small /? w (< 0.02), j c oc y Vq. 



Table 2 

Summary of threshold current for strength of extrinsic pinning. The third, fourth, 
and fifth columns indicate the origin of threshold (extrinsic (E) or intrinsic (I)), 
critical current, and mechanism of depinning (either spin transfer (ST) or force 
(/3 W )), respectively. The last column is a "good" variable to describe depinning. 
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12.4 Wall speed 



The wall speed after depinning is obtained as function of current based on 
Eq. (283) as in Fig. 25. The anomaly at j a are seen in Fig. 25 in the case of 



(3 W = 0.001 (at j — gg ~ 1.1) and (3 W = 0.02 (at j = 1), and the behavior is 
in agreement with numerical result of Ref. [48], where the anomaly of velocity 
under current was first reported. This anomaly is essentially the same as the 
anomaly under magnetic field, known as Walker breakdown [120,97]. 



12.5 Effect of external magnetic field 



Let us here briefly consider the effect of external magnetic field along easy axis. 
From Eqs. (265)(309), the field replaces (3 w ~j by (3 w j + b, where b = - 9 ^ A = 

— 2 '^g z (g = 2). The average wall velocity with V = is then given if j > j a 



P w3 + b sgn[[(P - %)j] - |] 
+ 1 + ft 2 \ 



a 




1, (315) 



X) - M±^, (316) 



(317) 



by 

and by 

N ' ft 

if J < ja., where the anomaly now occurs at 

~ b ± ft 

Ja ^ Pft - & 

the ± denotes sgn[(Pft — /3 W ) — 5/j a ]. 
12.6 Comparison with experiments 



Let us here try to explain the experimental results [70] assuming the wall 
there is planar and rigid. We first consider the case of regime I-a). Assuming 
£ ~ A, we estimate the pinning potential from the measured depinning field 
B c = 0.01 - 0.1 T as V = 0.34 x (10~ 2 ~ 10" 1 ) K = 4.7 x (10^ 26 ~ 10" 25 ) J, 
i.e., i = 1.4x (10~ 3 ~ 10~ 2 ), and so j c Ia) = (0.21 ~ 0.67) x j\. This value is 
still too big to explain the experimental value. Velocity jump is estimated as 
[51] Av 1 ^ = ^ x 839 m/s, so an extremely small (3 W (^ ~ 4 x 10" 3 ) is required 
to explain the experimental value of Av ~ 3 m/s [41]. If we assume regime 
I-b), the threshold is j c Ib ^ = t-j-t x 2.8 x (10^ 3 ~ 10^ 2 ) x f c . The experimental 
value could be reproduced if p w = 0.1 ~ 1. But such a large value of (5 W 
cannot be explained within the current understanding that j3 w arises from 
either non-adiabaticity [26,48] or spin relaxation [47]. Instead, Av cannot be 
explained by use of the above V assuming I-b) , as it predicts too large a value 
of Av Ib ) = 10 3 m/s. Thus, honestly, none of the above predictions based on 
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a rigid ID wall neglecting the temperature rise due to heating is successful in 
explaining the experimental result for metals quantitatively. 

There are some possibilities to resolve the disagreement. The most probable 
one would be the heating effect of the current, which has been reported to be 
important [70,77]. The estimate of V by use of experimental B c could be an 
over estimate if the effective barrier height Vq is greatly reduced by heating 
under current, while such heating does not occur under a static magnetic field. 
Let us estimate the pinning potential that gives the experimental value of j c . 
Assuming regime I-a), the experimental value of j c /Jl = 0.02 is reproduced if 
^ = |- = 1.3x 10~ 5 , which corresponds to V Q = 3 x 10~ 5 K= 4.5 x 10~ 5 T. 
This is two orders of magnitude smaller than the value extracted from B c . For 
I-b), we have \x = (3 W x 10~ 2 . From the experiment, Av/(—j c ) = (3 m/s)/(67 
m/s)= 0.05. This value is equal, for regime I, to — , so (3 W = 5 x 10~ 4 if 
a = 0.01. So in case I-b), ft = 5 x 10 -6 . Thus, assuming either regime I-a) 
or I-b), the experimental results could be explained by a pinning potential 
extremely weakened by heating, j&- = 10~ 6 ~ 10~ 5 . 

As far as the transverse wall is concerned, the assumption of rigidity seems 
not essential since the analytical results for a rigid wall [26,50] and numeri- 
cal simulation including deformations [48,76] predict similar behaviors. For a 
vortex or a vortex wall, the behavior is very different from that of rigid wall, 
as we see next. 



13 Motion of vortex 



In this section we briefly consider the motion of a vortex in two dimensions 
(Fig. 30). The dynamics of a vortex is very different from that of a planar 
(transverse) domain wall. Actually, the canonical momentum of vortex posi- 
tion X in the x-direction is Y, the position in the ^/-direction. This leads to 
vanishing of the intrinsic threshold, which is significantly different from the 
case of a domain wall, whose canonical momentum is 0o- This fact may explain 
why vortices or vortex walls move more easily than (transverse) domain walls 
[42,76]. 

The Lagrangian of a single vortex is given by Eq. (256) with a single vortex 
solution centered at X(t) = (X(t),Y(t)). It is given as [113] 

L v =-G-(XxX)-G-(v s xX)-U(X), (318) 

where G = Ge z {G = ^r&xy) is the gyrovector of Eq. (244), and the poten- 
tial energy U is from H§. The equation of motion in two dimensions is thus 
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Fig. 30. Examples of vortex structure. (Core size is very small.) The first two struc- 
tures are realized on a small circular disk [121,122], since the magnetic charge ap- 
pearing on the edge is zero. Structures denoted by blue and red have topological 
charges of n v = 1 and n v = — 1, respectively. 

given by [113] 

X-v s = ^e z x(F-aX), (319) 

where F = —^U(X) and a = aff dj d 2 x[{V9f + sin 2 0(V0) 2 ]- As seen, the 
vortex has a velocity perpendicular to the applied force, X _L F, which might 
sound strange but is an interesting feature of a topological object. The force 
due to current is given by Eq. (251), and so Eq. (319) reduces to 

a 3 ( f a 3 a ■ \ , 

X =2^- + e " X (>-4^; X j' (320) 

where / = 8 <%*, (2Se 2 N c p s + /iP/? sf ). The perpendicular force on the vortex 



■•• u 



of the adiabatic origin (j x G in Eq. (251)) is thus simply a spin torque, which 
induces the vortex to flow in the current direction. This is easily understood 
since the spin torque and the time-derivative term of the Landau-Lifshitz 
equation are combined into a Lagrange derivative along the current flow, d t —* 
(d t — 2^5 j s ■ VJ, which indicates that spin transfer torque drives any spin 
structure along the spin current. (This is seen from Eq. (256).) 

A very important consequence of vortex motion described by Eq. (320) is 
that the vortex (or vortex wall) has no intrinsic threshold [113]. In fact, as 
soon as any small current is applied, Eq. (320) indicates flow of the vortex, 
if extrinsic pinning is absent. This seems to be consistent with experimental 
and numerical observations [42,76]. 

Let us consider a single vortex in a film in more detail. The Hamiltonian of 
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localized spin is modeled as 

S 2 d 



H v = ^ J d 2 x[J[{V6) 2 + sin 2 #(V0) 2 ] + K ± cos 2 6]. (321) 

Since we cannot obtain the vortex solution analytically, we approximate it as 



7T 



V , n n ) 2 



r > A, 



tan-^dh^, 6={ ^_ 2 ,„,. ' .\ (322) 

2(1-6-1)' 



x " li ^=iy(l-e-'W)r<A, 



where r = v^ 2 + y 2 and A v is the size of the vortex core. Then the gauge field 
is given as A i ~ ±i Xl )^ $ e~ r ' Xv 6(X V — r), neglecting the small contribution 
from outside the core. Choosing the current direction as x, we obtain the 
Fourier transform as (A is the area of the film) 

\ 2 2 

A±{q) = ±ztt^(v^(1 - 2g 2 A v 2 ) + 2zgA v )e"^. (323) 

The resistivity due to the core is then calculated as 



^ / K 2 m 2 k F A 2 \ v 2 2c2kF 2 K 2 
Pv ~ 2 eVA 6 ' [6 ] 

where we used the fact that the resistivity is dominated by the contribution 
from q < k F+ — k F _ = 2k F (. The nonadiabatic force is then given as 

F ref ~ 3 -h (— V (k F \ v ) 2 de^ 2kF2x ^. (325) 

e \e F J 

On the other hand, the force due to the topological Hall effect is given by 

F Ha11 ~ 3 -h (326) 

e 

for a vortex with n v — 1. 

Let us estimate the magnitude of the Hall effect. The Hall resistivity is given 

as 

27iS 2 Pn v 

*■» " ^i^r- (327) 

The Hall conductivity is given by a xy = <J p , xv ' p0 , 2 ( cr = Pq 1 is the Boltz- 
mann conductivity) and hence the ratio of the Hall current to applied current 
is obtained as 

IL = EEL „ 2nS 2 Pn vT ^-. (328) 

Jo Po k F A 

Thus, the deviation of the electric current becomes significant if the ratio — 
is of the order of unity. The deviation of current due to the Hall effect would 
be large in clean samples. 
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Current-induced dynamics of a vortex in a disk was observed in Ref. [122]. 



14 Spin torques in LLG equation 



In this section, we describe the present theoretical status of the microscopic 
derivation of Landau-Lifshitz-Gilbert (LLG) equation in the presence of cur- 
rent. We focus on magnetization dynamics that is slowly varying in space and 
time, as described by the LLG equation with local torques. Here 'slow' means 
slow compared to the electronic scales, so it is satisfied quite well in most cases 
for metallic systems. One important aim here is to calculate microscopically 
the /3 sr term and damping due to the electron spin relaxation. Throughout 
this section, we consider general spin configurations. 



14-1 General 



The LLG equation is given by 

^ = 70 H - xM + ^Mx^ + T ol , (329) 

at M at 

in terms of the magnetization vector M. The first term on the right-hand 
side represents the precessional torque around the effective field H e g, with 
7o(= ^p = — jp- > 0) being the bare gyromagnetic ratio (without modifica- 
tion by conduction electrons). The effective field includes the ferromagnetic 
exchange (gradient energy) field, the magnetocrystalline anisotropy field, and 
the demagnetizing field. The second term represents the Gilbert damping, 
coming from processes that do not involve conduction electrons, and is thus 
present even in insulating ferromagnets. The effects of conduction electrons 
are contained in the third term, T c \, called the spin torque in particular. This 
term comes from the s-d exchange coupling H s d to conduction electrons, and 
is given by 

T cl = -Mn(r)x((r(r)} ne . (330) 

(We wrote here explicitly ( ) nc to indicate that the expectation value is taken 
in the non-equilibrium state with current flow or dynamical magnetization. 
See §14.2.2.) 

For notational convenience, we introduce a unit vector, n, whose direction is 
in the d-spin direction, hence is opposite to magnetization direction, 

hS 
M = - la -rn. (331) 



a 3 
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Note that the magnetization is, by definition, the magnetic moment per unit 
volume, hence \M\ = H^qS/o?. In terms of n, the LLG equation is written as 

n = 70-Hcff x n + a n x n + t' ch (332) 

where we have put T C \(M) = —(HS / a 3 )t' cl (n) . The dot represents the time 
derivative. 

For long-wavelength, low-frequency dynamics, it may be sufficient to consider 
spin torques that are first order in space/time derivatives. Let us call such 
torques as adiabatic torques. In the presence of rotational symmetry in spin 
space, they are expressed as 



'ad 



5S 



fo°- V) n- j3 sr nx (v s ° ■ V) n - a sr (n x h) - — n. (333) 



,5' 



The first term on the right-hand side is the celebrated spin-transfer torque 
[114,123,115,117,49], where 

< = £- s i, ^ (334) 

is the (unrenormalized) "spin-transfer velocity," , with j s being the spin-current 
density. The second term, sometimes called the '/3-term', comes from spin- 
relaxation processes of electrons [47,48,54,28,29,64,53]. Here (3 sr is a dimen- 
sionless constant. The third term is the Gilbert damping, which also results 
from spin relaxation of electrons. The fourth term contributes as a "renormal- 
ization" of spin, as seen below. 

If the magnetization varies rapidly, we have in addition a non-adiabatic torque, 
r° a , which is oscillatory and nonlocal (see §14.3.4). The total torque may thus 
be given by the sum of the two, 

.0/ i _o/ 



<i = < + r n :. (335) 

The LLG equation (332) is then written as 

Crn\ 
1 + ~o ) h = 7o-Heff x n ~ ( v s- v ) n ~ Pa n x ( v s ° -V) n - (a + aw) nxfi+T^- 

(336) 
Here we have transposed the 'spin renormalization' term to the left-hand side. 
We define the total ( "renormalized" ) spin as 



S tot = S + 5S, (337) 

with 5S being the contribution from conduction electrons (Eq. (266)), and 
divide both sides of Eq.(336) by S to t/S. Then we arrive at 

n = 7if c ff x n — (v s ■ V) n — j3 sr n x (y B -V) n — a (n x h) + Tj' a , (338) 
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where 7 = (S/S tot ) 70, a = (S/S to t)(a + a sr ), T„ a = (S/S to t) T° a , and 



5 



a 



,3 



u, 



5 V ' = 2^~ j " (339) 

is the "renormalized" spin-transfer velocity. (These torques are schematically 
shown in Fig. 10.) 

In the parameter space of the LLG equation, Eq. (338), the manifold of a = (3 sr 
(with T^ a = 0) provides a very special case for the dynamics, and there has 
been a controversy whether the relation a = (3 SV holds generally or not. If a = 
/3 sr , the following peculiar dynamics are expected, (i) Any static solution, n(r), 
in the absence of spin current is used to construct a solution, n(r — v s t), in 
the presence of spin current v s . (ii) The current-induced spin- wave instability 
does not occur [28]. 

The relation a = /3 sr was originally suggested in Ref. [54] based on the Galilean 
invariance of the system. Although one may argue that the Galilean invari- 
ance should be valid for the long-wavelength and low-frequency dynamics in 
which the underlying lattice structure is irrelevant, the a and /3 sr come from 
spin-relaxation processes [47], which are usually intimately related to the lat- 
tice, e.g., through the spin-orbit coupling. The problem is thus subtle, and one 
has to go beyond the phenomenological argument such as the one based on 
the Galilean invariance. Instead, a fully microscopic calculation, which starts 
from a definite microscopic model and does not introduce any phenomeno- 
logical assumptions, is desired. The present section is devoted to outline such 
attempts. 

At present, only a single model, where spin-relaxation processes are intro- 
duced by magnetic impurities, has been examined, with a result that a 7^ f3 ST 
(even for single-band itinerant ferromagnets) in general. This will be surveyed 
in §14.2 and §14.3. Studies on other models, hopefully with more realistic 
spin-relaxation mechanisms, are left to future studies. Readers who are not 
interested in theoretical details but the results can jump to §14.2.3. 



14-2 Small- amplitude method 



14.2.1 Microscopic model 

Let us first set up a microscopic model. We take a localized picture for fer- 
romagnetism, and consider the so-called s-d model. It consists of localized d 
spins, S, and conducting s electrons, which are coupled via the s-d exchange 



95 



interaction. The total Lagrangian is given by L tot = L$ + L c \ — H s d, where Ls 
is the Lagrangian for d spins (Eq. (21)), 



d fv 
dt 2m 



ih-k; + 77— V 2 + e F - V imp - V BI 



(340) 



L e \ = I d 3 x c' 

is the Lagrangian for s electrons, and 

H sd = -M f d 3 xn(r)-a(r) (341) 



is the s-d exchange coupling. Here, c* = (cl,c!) is the spinor of electron cre- 
ation operators, and cr{r) = c*{r)crc{r) represents (twice) the s-electron spin 
density, with er being a vector of Pauli spin matrices. We have put S = Sn 
with the magnitude of spin, S, and a unit vector n, and M = J s dS with J s d 
being the s-d exchange coupling constant. (We define n to be a unit vector in 
the direction of spin, which is opposite to the direction of magnetization. ) 



Fig. 31. Diagrammatic representation of spin relaxation due to spin flip scattering 
by impurity spin. 

The s electrons are treated as a free electron gas in three dimensions subject 
to the impurity potential 

V irap + V sr = uJ2$(r- Ri) + u s J2S imp>j -(rS(r - R'^. (342) 

« j 

The first term describes potential scattering. The second term represents 
quenched magnetic impurities, which is aimed at introducing spin relaxation 
processes (Fig. 31). The averaging over the impurity spin direction is taken as 

"imp,i"inipj = o^imp"*i" ■ (343) 

The damping rate of s electrons is then given by 

7 ff = -— = nriiU 2 ^ + -n s u 2 s Sl ''(2u a + v a ). (344) 

Zr a 3 

Here n\ (n s ) is the concentration of normal (magnetic) impurities, and v a = 
mkF a /2n 2 h 2 (with hkp a = ^2meF a ) is the density of states (per volume) 
at energy e-pa = £f + oM. (The subscript a =f, [ corresponds, respectively, 
to a = +1,-1 in the formula, and to a =|, t or — 1, +1-) We assume that 
7cr ^C e-pa and 7 CT <^C M, and calculate the torques in the lowest nontrivial 
order in ^a/e-pa and 7 CT /M. 
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In this section, the electron lifetime is treated as spin-dependent, since this 
feature becomes essentially important in calculating the effect of spin relax- 
ation consistently. (Determination of /3 sr and a terms requires much care, and 
so phenomenological argument can easily lose consistency.) 



14-2.2 General framework 

The spin torque from H s d is given by 

t A (r) = Mn(r) x (<r(r)) M . (345) 

This is related to T cl of Eq.(330) via t cl [n] = -T cl [M] , and to t' cl of Eq.(335) 
via 

t e i = -rt' cl (346) 



a 3 



with a being the volume per d spin. 



The calculation of spin torque is thus equivalent to that of s-electron spin po- 
larization, (cr(r)) nc , or precisely speaking, its orthogonal projection (<Tj_(r)) ne 
to to. (We define &± = a — n(n-er) and <x^_ = & — z(z-a). Note that (<r) nc in 
Eq.(345) can be replaced by (<r_i_) ne . ) The expectation value (• • -) ne is taken 
in the following nonequilibrium states depending on the type of the torque. 

(a) Nonequilibrium states under the influence of uniform but time- dependent 
magnetization. This leads to torques with time derivative of to, namely, Gilbert 
damping and spin renormalization. 

(b) Nonequilibrium states with current flow under static but spatially-varying 
magnetization. This leads to current- induced torques, namely, spin-transfer 
torque and the /3-term. 

In the presence of spin rotational symmetry for electrons, adiabatic spin 
torques, which are first order in space/time derivative, are expressed as 

r a ° d = a h + (a- V) to + b (n x n) + n x (b • V) n. (347) 

The corresponding s-electron spin polarization is given by 

(cr_L) ne = — [b n + (b-V) n-a (nxh) -nx (a -V) to] . (348) 

To calculate the coefficients a M and b^ microscopically, it is sufficient to con- 
sider small transverse fluctuations, u = (u x ,u y ,0), \u\ <C 1, around a uni- 
formly magnetized state, to = z, such that to = z + u + 0(u 2 ) [28,124]. Then, 
up to 0(u), Eq.(348) becomes 

(cr±)nc = -^ [b u + (b-V)u - a (z x u) - z x (o- V)«] . (349) 
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This equation can be regarded as a linear response of er± to u, and the co- 
efficients, a^ and b^, are obtained as linear- response coefficients. (Precisely 
speaking, (cr_|_) ne due to current is calculated as a linear response to the ap- 
plied electric field [52].) 



14. 
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Results 
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2Mr, 



(350) 

(351) 

S 
a = a sv + a 0) (352) 

0+00 

(353) 



where p s — n-f — nj_, and j' s = cr s -E = jj — j^ is the spin current, with cr s = 
(e 2 /m)(n^Ti — njTj_) being the "spin conductivity". (n a is the density of spin-cr 
electrons.) We have defined the spin-relaxation time r s by 

H 4:7V 

— =— n s u s Sf mp (u+ + z/_) (354) 

To o 



As expected, only the spin scattering (~ r -1 ) contributes to a and f3 BT , and the 
potential scattering (~ niu 2 ) does not. (For a, the second term on the right- 
hand side of Eq.(352) comes from processes that do not involve s-electrons.) 

For a single-band itinerant ferromagnet (as described by, e.g., the Stoner 
model), the results are obtained by simply putting S = in Eqs.(350)-(353), 
and by using the spin polarization of itinerant electrons for 5S. Even in this 
case, we see a sr 7^ j3 sr ; however, it was pointed out [28] that the ratio 

&E _ Ps a 1 + j (M\ 2 



(355) 
a sr M(z/+ + z/_) 12 \e F J K ' 

is very close to unity. Even in this case, if we make the impurity spins anisotropic 
by generalizing Eq.(343) to 



S\ (a,(3 = x,y) 



SfSf = «%^x <! — v " ( 356 ) 
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we obtain 

Psr __ 3Dj_ + S z (V\7\ 

a sr 2(5i+^f) 

which ranges from 1/2 (for Sj_ < Sf) to 3/2 (for Sj_ > Sf). From this 
example, we can learn that the ratio f3 sr /a sr is very sensitive to the details 
of the spin-relaxation mechanism. In reality, a non-electron contribution, a , 
exists, and so a = /3 sr predicted in Ref. [54] would never happen. 



The results obtained based on the phenomenological equations [47] can be 
written in the form 



5S h . . 

azL = sW^' (358) 

whereas /3zl = /3 ST is the same as Eq.(353). Thus, it predicts a = (3 sr for single- 
band itinerant ferromagnets, S = 0, which is, however, in disagreement with 
the present microscopic calculation, (350)-(353), predicting a ^ (3 sr in general. 



14-3 Gauge field method in the presence of spin relaxation 



In the previous section, we considered small-amplitude fluctuations of mag- 
netization, and calculated the torques in the first order with respect to these 
small fluctuations. In this sense, the spin torques calculated there are limited 
to small- amplitude dynamics. (Only for systems with rotational symmetry in 
spin space, where the form of the torque is known as Eq.(347), is this small- 
amplitude method sufficient to determine the coefficients, hence the torque.) 
In this section, we describe a theoretical formalism that is not restricted to 
small-amplitude dynamics, but can treat finite-amplitude (arbitrary) dynam- 
ics directly [52]. 



14-3.1 Adiabatic spin frame and gauge field 

To treat finite-amplitude dynamics of magnetization, we work with a lo- 
cal/instantaneous spin frame (called "adiabatic frame" in the following) for s 
electrons whose spin quantization axis is taken to be the local/ instantaneous 
rf-spin direction, n [125-127]. The electron spinor a(x) in the new frame is re- 
lated to the original spinor c(x) as c(x) = U(x)a(x), where U is a 2 x 2 unitary 
matrix satisfying c^(n-cr)c = a)o~ z a. It is convenient to choose U satisfying 
U 2 = l. 

Since d^c = £/(<9 M + Wd^^a = £/(<9 M + iA^a, the Lagrangian for the a- 
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ih (dt + iA ) + — (Vi + iA t f + Ma z - V imv - V m 
2m 



electrons becomes 

L[a] — / cPxa 1 

' (359) 

The original electrons moving in time-dependent /inhomogeneous magnetiza- 
tion are thus mapped to new electrons moving in a uniform and static mag- 
netization Ma z but there arises a coupling to an SU(2) gauge field 

Ap = -irffaU) = Ay a = A^cr. (360) 

Here A M is a measure of temporal (/i = 0) or spatial (/x = 1,2,3) variation 
of magnetization. (We use the vector (bold italic) notation for the spin com- 
ponent. The space-time components are indicated by subscripts such as //, v 
(=0,1,2,3) or i, j (=1,2,3).) 

Let us introduce a 3 x 3 orthogonal matrix H, representing the same rotation 
as U but in a three-dimensional vector space, and satisfying detTl = 1. Note 
that Tlz = n, Tin = z, rferc = Tl (a) era) , and that 1Z(axb) = (Ha) x (Tib) for 
arbitrary vectors a and b. Then the spin-torque density, Eq.(345), is written 

as 

t cl (x)=MTl(zx(cr(x)) nc ), (361) 

where cr(x) = (a)cra) x . 

Since the gauge field A yi contains a space/time derivative of magnetization, 
one may naturally formulate a gradient expansion in terms of A M to calculate, 
e.g., the torque (or spin polarization). In particular, the adiabatic torques are 
obtained as the first-order terms in A^. 

(*±) m = -| [a,Ai + \{z x A^)\ . (362) 

Here a± = er — z(z-cr) and Aj; = A^ — z(z-A^) are the respective trans- 
verse components, and the sums over /i = 0, 1, 2, 3 are understood. From the 
identities, 

TlAf l =-^nx(d,n), Tl(z x A$ = \d,n, (363) 

together with Eq. (361), we see that Eq. (362) leads to the adiabatic torque 
density r a ° d of Eq. (347). 

The processes contributing to the f3 term are shown in Fig. 32. [52] 
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Fig. 32. Diagrams contributing to the (3 term. The electron Green's functions here 
are denoted by thick lines, indicating that they include self-energy from spin flip 
interaction, H ST . The last three diagrams are vertex corrections due to spin flip 
scattering. Since we are interested in adiabatic limit, the gauge field does not change 
electron momentum. 

14.3.2 Results 

If we regard Eq.(362) as a linear response to the gauge field A^ appearing 
in Eq.(359), the coefficients a M and b^ are calculated as linear-response co- 
efficients, and 5S, v s , and f3 ST thus obtained coincide with those obtained by 
the small-amplitude method, Eqs.(350), (351), and (353). However, it predicts 
a sr = 0, namely, it fails to reproduce the Gilbert damping. 



14-3.3 Gilbert damping 

The above puzzle with the Gilbert damping can be resolved if we note that the 
impurity spins, which are static (quenched) in the original frame, become time- 
dependent in the adiabatic frame. Namely, the spin part of V BI is expressed 

as 



^imp,j ' C " C 



Sj(t)-a^cra 



where 



Sjit) = n(t)Si 



impj 



(364) 



(365) 



is the impurity spin in the adiabatic frame, which is time-dependent. (This 
fact is expressed by V ST in Eq.(359).) Actually, we can obtain the gauge field 
from this time dependence as 



m)n(t)} 



a/3 



1e a ^Al. 



(366) 
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Explicit calculation of (cr±) nc in the second order in Sj(t) (nonlinear response) 

gives 

2ttH 

leading to the Gilbert damping, with damping constant exactly the same as 
Eq. (352). The Gilbert damping term is shown diagramatically in Fig. 33. 



(<x±(*)> ne = -^jn s u 2 8 Sl mp ul{zxA^(t)), (367) 



The present calculation provides us a new picture of Gilbert damping. While 
the s-electron spin tends to follow the instantaneous ci-spin direction n(t), it is 
at the same time pinned by the quenched impurity spins. These two competing 
effects are expressed by the time dependence of Sj(t) in the adiabatic frame, 
and this effect causes Gilbert damping. Namely, the Gilbert damping arises 
because spins of s-electrons are 'dragged' by impurity spins. 

Generally, any terms in the Hamiltonian leading to spin relaxation break spin 
rotational symmetry of s electrons, and thus acquire time dependence in the 
adiabatic frame. Therefore, the same scenario as presented here is expected 
to apply to other type of spin-relaxation processes quite generally. 




Fig. 33. Diagrams contributing to damping a in the gauge transformed frame. The 
gauge field (wavy line) is induced by dynamical impurity spins (represented by 
double dotted line). 

14-3.4 Non-adiabatic torque 

The non-adiabatic torques come from large-momentum processes, and are ex- 
pected to be insensitive to impurities, so let us put V imp — V ST — for sim- 
plicity. If we retain the full q-dependence in the response function, we have, 
among others, a contribution 

(&±(q)U = XiM) E i (* x Kj) + ■ ■ ■ > ( 368 ) 

where the function Xij(o) * s nonzero only when |fcpT — &FJ.I < q < &ft + &FJ.- in 
particular, Xijio) = in the vicinity of q = 0, and hence the long- wavelength 
approximation cannot be applied. The resulting torque 



rL{x) - 



MEi J d 3 x' Xij (r - r')n(x)Af(x'), (369) 



is characterized by an oscillatory function Xij{ r — r ') i n rea l space, and is es- 
sentially nonlocal [62,26,49,50]. This is the non-adiabatic, momentum-transfer 
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torque due to electron reflection [26], generalized here to arbitrary magneti- 
zation texture. For details, see §3.2.1, §9.1 and Ref. [50]. 



14-3.5 Further references 

The spin-transfer torque in the form of Eq.(333) was first derived by Bazaliy 
et al. [114]. The /3-term was first derived by Zhang and Li [47] based on a spin 
diffusion equation with a spin-relaxation term included, and by Thiaville et 
al. [48] as a continuum limit of a special type of torque known in multilayer 
systems [128]. 

Duine et al. put the present microscopic (small-amplitude) calculation into the 
Keldysh formalism, and developed a functional description of spin torques, 
which will be used for finite-temperature and/or fluctuation dynamics [53]. 
An attempt with the Boltzmann equation is carried out in [64]. Some devel- 
opments in the treatment of Gilbert damping can be seen in [124,129]. The 
non-adiabatic torque is studied in [62,26,49,50]. The domain-wall resistivity is 
studied in [127,130]. 

The effects of spin-orbit coupling is studied on domain- wall resistance [131] 
and domain-wall mobility [66] in a model of ferromagnetic semiconductors. 

Microscopic understanding of each spin torque is also an important issue. In 
particular, the dissipative torques (a and (3 terms), coming from spin relax- 
ation of electrons, should be understood as material-dependent quantities; for 
each real system, we need to identify the dominant spin-relaxation mechanism, 
and clarify the dependence of a and (5 on material parameters. Development 
of first-principles calculational methods for spin torques is also desired for the 
purpose of material/device design. 



15 Domain wall and electron transport 



In this section, we briefly review electron transport theory in the presence of a 
domain wall. Electron transport coefficients such as resistivity and Hall coeffi- 
cients represent the reaction of current- induced forces on magnetic structures. 
We calculate these coefficients in linear response theory, and show that the 
current-induced forces we have discussed in §9.3 are indeed proportional to 
them. For experimental developments, see Ref. [17]. 
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15.1 Resistivity in magnetic metals 



Since more than a century ago a number of studies has been carried out 
on electric transport properties in ferromagnetic metals. They revealed many 
remarkable features that are not seen in non-magnetic metals. One of the most 
notable would be the hysteretic and anisotropic behavior of the resistance in 
the magnetic field (magnetoresistance) observed at small magnetic fields of 
< 1 T, which has been already noted more than a hundred years ago [1,132]. 
The magnetoresistance in the case of the field H parallel to the current / 
takes a minimum at a finite value of the field (~ 200 Oe for instance for 
the case of Ni and Fe) . If the field is applied perpendicular to the current, 
the curve of magnetoresistance is reversed; namely the resistivity shows a 
maximum at a certain field and decreases as the field deviates from that value. 
The hysteretic behavior of the magnetoresistance is due to the fact that the 
resistivity is mostly governed by the total magnetization vector M, not the 
applied magnetic field. (The total magnetization is given as the average of 
localized spin vectors, M = —g^y J d 3 xS.) The observed resistivity p as a 
function of the field, H, has been shown to be well fitted by a phenomenological 
relation p ~ p + Ap ani < cos 2 ^ >, where p is the field independent part 
and Ap ani measures the strength of the anisotropy in the resistivity [2,133]. 9m 
is the mutual angle between the local magnetization and the current, which 
depends on the magnetic field, and the brackets denote the average over the 
sample. In most ferromagnetic metals, the anisotropy Ap an i is positive. This 
anisotropic behavior of the resistivity is called anisotropic magnetoresistance 
(AMR). 




Fig. 34. Anisotropic magnetoresistance (AMR) effect. Resistivity depends on the 
angle between the magnetization M or localized spin S and applied current, i. 

The mechanism of AMR was argued to be spin-orbit interaction by Smit [133] 
and McGuire and Potter [2]. They pointed out that the conduction (s-) elec- 
tron is coupled to the magnetization due to the scattering into the magnetic 
(d-) band and the spin-orbit interaction there, and that this process gives rise 
to a spin asymmetric lifetime, which depends on the angle between the cur- 
rent and the field, resulting in an anisotropy. According to their arguments, 
the positiveness of Ap ani is explained if the resistivity is dominated by the 
electron having the minority spin. 
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15.2 Domain wall contribution to classical magnetoresistance 



The magnetoresistance observed in bulk materials is mostly understood well 
in terms of the AMR effect, which assumes that the magnetization changes 
very slowly and hence the electron feels only the average magnetization [2]. 
However this assumption may not be good in small magnets that contain many 
domains with different directions of local magnetization. In fact the boundary 
of these domains is a structure called a domain wall where the spins rotate 
spatially within a finite distance. Such domain walls lead to scattering of the 
electron, which is not taken into account in the standard AMR argument. In 
this paper we exclusively consider the effect of domain walls on the electronic 
resistivity. Of particular interest would be the case of a small sample less than 
a typical domain size ~ 1 ~ 10 /xm, since there the magnetic properties can be 
described in terms of the domain wall configuration. In fact in such samples the 
magnetization process as the field is swept will be described by the nucleation 
of one or a few domain walls, motion of the walls by depinning followed by the 
annihilation. The effect of domain walls on the resistivity would then appear in 
the magnetoresistance as discrete jumps in the measurement with the field is 
swept. (Such events of domain walls are faster (e.g., the speed of domain wall 
motion is estimated to be about 182 m/s in submicron wires of NiFe [134]) 
compared to the sweep speed.) These effects are similar to the Barkhausen 
noise [135]. The jump will not be seen very clearly in bulk samples, since 
there the contribution from a macroscopic number of domain walls, whose 
nucleation and pinning energy may differ, will be summed up in the observed 
magnetoresistance and thus the contribution from each domain wall will not 
be visible. Jumps in the magnetoresistance was first observed in 1994 in a Ni 
wire of about the diameter of 300 A[136,19]. Nowadays resistance of a single 
domain wall has been measured in many systems [17,85]. 

Let us here give a rough estimation of the effect of the wall on the classical elec- 
tron transport. "Classical" here means that the conductivity (or resistivity) 
is evaluated at the lowest order in h, which is determined by the probability 
of reflection by the wall. The most important parameter of the wall on the 
transport is the thickness of the wall, A. This quantity is determined by the 
competition between the exchange energy, J, which aligns the neighboring 
spins and the magnetic anisotropy energy in the easy axis, K, which tends 
to make the wall thinner to minimize the deviation of spins from the easy 
axis, as A = yJ/K (Here J has dimensions of J/m 2 .) Thus A depends on the 



material and also on the sample shape since K depends on the shape. In the 
case of 3d transition metals such as iron and nickel A ~ 500 ~ 1000 A[19], 
and A ~ 150 A in Co thin films [21]. This length scale is very large compared 
with the length scale of the electron, kp 1 ~ 0(1 A) (kp being the Fermi wave 
length of the electron). Thus in such materials the conduction electron can adi- 
abatically adjust itself to the local magnetization at every point as it passes 
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through the wall, resulting in a very small scattering probability by the wall. 
The classical resistivity in the Boltzmann sense, which is proportional to the 
reflection probability by the wall, is thus expected to be negligibly small. This 
was explicitly shown in 1974 by Cabrera and Falicov [137] by calculating in the 
clean limit the reflection coefficient based on the one-dimensional Schrodinger 
equation for the electron coupled via exchange coupling to the magnetization 
whose configuration is a domain wall. They obtained for the case of thick wall, 
kp\ ^> 1, the expression _R W oc exp(— 2nk F X) for the wall contribution, and 
showed that the wall can have a large effect only if the wall is extremely thin 
(kp\ < 1) and if the conduction electron is strongly polarized by the magne- 
tization. Such a thin domain wall was fictitious at that time, where 3d metals 
in the bulk were mostly considered, but can be realized at present [23,22]. 

In the past 10 years there has been a renewal of interest in the theory of 
the classical resistivity due to the domain wall [138,127,139-147], because ex- 
perimental control of the wall thickness [23] and precise measurements are 
becoming possible. (Resistance due to a single domain wall in GaMnAs was 
experimentally determined to be ~ 1 Q [85].) In Ref. [127], the domain wall 
resistivity was studied based on a linear response theory. There the effective 
wall-electron interaction was derived and treated perturbatively. The calcula- 
tion corresponds to the perturbation expansion with respect to 0(/cfA) _1 . The 
effect of the finite mean free path due to the impurity scattering was taken 
into account consistently. As far as the classical transport is concerned, the 
calculation confirmed the result of Cabrera and Falicov. It was pointed out by 
Levy and Zhang [139] that in the presence of a spin asymmetry in the electron 
lifetime, which is the case in most ferromagnets, the domain wall can have a 
substantial effect on the classical resistivity by mixing the two spin channels 
with different resistivities, in agreement with experiments on Co films at room 
temperature [21]. The effect of lifetime asymmetry has been further discussed 
in detail in Refs. [140,141]. The domain wall resistivity in the ballistic limit has 
been discussed by use of realistic band structures in Ref. [142] . It was shown 
there that the existence of nearly degenerate bands at the Fermi level in real 
magnets enhances the classical resistivity due to the wall. The calculation of 
domain wall resistance by use of Landauer's formula was compared with that 
by linear response in Ref. [148]. 

Of recent particular interest is an atomic scale contact of magnetic metals 
in the ballistic region, as experimentally realized in Ref. [23]. In such narrow 
contacts the profile of the wall is determined mostly by the shape rather than 
the anisotropy energy of the bulk magnets, and thus the wall is trapped in a 
contact region, which is typically on a nm scale [142,143,149]. The adiabaticity 
does not hold in such cases of small kpX comparable to 1, and the wall can 
have a large effect on resistivity [137,142,143]. Indeed a magnetoresistance of 
200% was observed in 1999 in ballistic Ni nanocontacts [23] where the number 
of the channels is less than N < 10, and this was interpreted in terms of a 
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strong reflection by a nanometer scale domain wall [143]. In dirty contacts, 
the effect is not so large, since the wall scattering is smeared out by impurity 
scattering [150]. More recent studies report magnetoresistance over 1000% 
[151]. Theoretical studies, however, suggest that such high magnetoresistance 
is not of electric origin and could be due to mechanical effects in the contact 
region [152]. Whatever the mechanism, such high magnetoresistance would be 
useful for applications. 



15.3 Quantum electron transport and domain wall 



Most experimental studies have been carried out on low resistivity materi- 
als, and theoretical studies also have mostly been based on classical transport 
theory (i.e., neglecting 0(h/epr)). However, besides the classical transport, 
there is another important aspect of electronic transport at low temperature. 
This is the effect of the quantum coherence among electrons, which modi- 
fies the low-energy electronic properties significantly. The effect becomes im- 
portant in disordered metals with high resistivity, where £ becomes shorter, 
£ = {hkpTJm) being the elastic mean free path. The electron wave scattered 
by such normal impurities can interfere with the incoming wave, leading to 
a standing wave. This is called weak localization and the resistivity in this 
case is enhanced due to the quantum interference [153,154]. This correction 
becomes large in small dimensions such as in wires since there the interference 
becomes stronger. The interesting point of this situation is that the electronic 
properties are very sensitive to a small disturbance because of the presence of 
coherence. For instance in non-magnetic metals of micron size, even the mo- 
tion of a single impurity atom has been shown to change the low-temperature 
conductance (G = o~A/L, A and L being the cross-sectional area and length of 
the system) of the entire system by disturbing the coherence [155]. The mag- 
nitude of the conductance change turns out in most cases to be a universal 
order of e 2 /h [156]. 

Because of this sensitivity these conductance fluctuations as a consequence of 
quantum interference have already been used as a probe in studies of various 
mesoscopic non-magnetic metallic or semi- conducting systems. For example 
the telegraph noise due to a two-level oscillation of a defect in Bi films has been 
investigated and it turned out that the oscillation at T < 1 K is governed by 
the quantum tunneling subject to the dissipation from the conduction electron 
[157]. Such measurement of the quantum transport properties has been proved 
to be a useful probe also for studies of mesoscopic spin-glass systems [158,159] 
and the magnetization flip of mesoscopic magnets [160,161]. 

It would be natural to expect that the domain wall also affects the quantum 
transport. The effect of domain wall scattering on the quantum correction to 
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the conductivity was first studied in Ref. [127]. It was shown that a domain 
wall destroys the interference among the electrons and the wall contributes 
to a negative quantum correction to the resistivity. In disordered 3d transi- 
tion metals, this quantum correction can overcome the classical contribution 
(i.e., reflection), and thus a wall may in total lead to a decrease of resistivity. 
Numerical simulation also supports the negative resistivity contribution from 
the wall in disordered thin wires [162]. It has been also pointed out that the 
geometric phase attached to the electron spin as it passes through the wall 
can also cause an important dephasing effect, which would become important 
in multiply connected geometries [163,164]. The effect of dephasing due to the 
magnetic origin has been considered also in a thin film of metal sandwiched 
by ferromagnetic layers [165]. There the internal magnetic field at the inter- 
face causes dephasing in the conduction layer in the presence of the spin-orbit 
scattering, which contributes to a positive magnetoresistance. 



These theories on the quantum transport in magnetic metals are based on the 
existence of electron coherence in magnetic metals. In non-magnetic metals 
such coherence has been observed as an appearance of weak localization, for 
instance, in thin films of Cu [166] and non-magnetic metals with magnetic 
impurities [167]. Magnetic systems have been considered as not suitable for 
coherent electron transport, since the magnetic field (or magnetization) and 
magnetic disorder result in decoherence. However, in ferromagnets, magnetic 
disorder should be frozen by a strong internal field (magnetization) and would 
not destroy coherence as theoretically discussed in Refs. [168,169]. The strong 
internal field in ferromagnets (M ~ 0(1 T)) would not affect the coherence 
in a mesoscopic case like a sufficiently small wire, since there the magnetic 
flux penetrating through the wire can be small enough. Nevertheless, there 
have not been so far many observations of electron quantum coherence in fer- 
romagnetic systems. A negative contribution to resistance in the presence of 
domain walls was observed in GaMnAs at 4.2 K [170], which was explained 
as due to electron decoherence due to the domain wall predicted in Ref. [127] 
in the weak localization regime. Semiconductor ferromagnets were found to 
show further evidence of quantum transport: universal conductance fluctu- 
ation [171]. Observation in metallic systems seems more difficult. Aharonov- 
Bohm oscillation was observed in a 500 nm FeigNigi ring [172,173]. Mesoscopic 
resistance fluctuation due to phase arising from a domain wall was observed 
in Co nanoparticles [174]. 



In the present paper, we will not discuss further the subject of quantum trans- 
port. 
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15. J,. Linear response theory of resistivity and Hall effect 



In this section, we derive the resistivity (both diagonal and off-diagonal (i.e., 
Hall) components) due to the general spin texture on the basis of linear re- 
sponse theory, using the Mori formula [175-177]. The Mori formula is valid in 
the clean limit, r — > oo, where the resistivity is dominated by spin structure. 
It relates the resistivity p si j (ij being spatial directions) to the correlation of 
random forces in the weak scattering case as 

fe 2 n\~ 2 1 

*« = {-^) £^Mw,(fcO - xjjM- ( 37 °) 

Here, Xjj^d = -(h/fiV) < Ji( iu >t)Jj(-iue) > with j, t = dJi/dt = ±[H, J;], 
where H is the total Hamiltonian and Jj = — J2k ^i^feCfc i s the total current. 
The operators here are defined in the imaginary time, r = it. The correlation 
function Xjj{^) in Eq. (370) denotes an analytical continuation of the corre- 
lation function calculated for imaginary frequency, i.e., Xjji^) — Xjji^e ~* 
fru + iO). 

The nonconservation of the current (i.e., finite J) arises from the scattering 
by the spin texture. In fact, Eq. (118) leads to 






-2M^2Ar a a\ +q a ff a k 



+|- (Af((2k + q)-q)- qi (A a ■ (2k + q)) - hq t A«) a\ +q a a a k 
2m 

(371) 

where we neglect higher order terms in A. We consider a static spin texture 
(i.e., A = 0), and then 



- 22 (-) M EE4 :ff «WA (372) 

k,q cr 



The Fourier transform in imaginary time r, defined as Ji(iue) = J e luJeT Ji(r), 

2ir£ 



where Ui = ^ is a bosonic thermal frequency, is given as 



Ji(wi) = -z2 ( — ) M J2 J2 A a ^k+ q ,n+e(T<rak,m (373) 

k,q,n o" 



where a kn — ~~i= ft ^ e luJnT ak(T)dr and uj n represents the fermionic frequency 
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n = ^~j — '-. The thermal correlation function is then obtained as (Fig. 35) 
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O i j i ^i)= - n {—) yE^(9) A r(-9)«E G ^, n+ i,-An ff , (374) 
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where the imaginary time Green's function is defined as 

1 



G 
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(375) 




k+f co„+(q 



Fig. 35. Diagram representing the resistivity in the Mori formula at the lowest order 
in the gauge field. 



The summation over u n is carried out to obtain 
V m J V 



xu<m > = -^"I'kAn^n-.) ^^ . (376) 

kqa £k+q,-a — ^k,a — VjJf, 



Thus the resistivity is obtained as 
AM 2 1 



Psij = lim 
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E Im 
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. (377) 



W e k+q,-a — £k,a ~ UJ — iO _ 

Then, choosing the current direction as z, we may rewrite the resistivity as 

4ttM 2 1 
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X; \A:(q)\ 2 5(e k+g ,_ a - e kjff )5(e k!(7 ). (378) 
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Thus, the reflection force of Eq. (233) is proportional to the resistivity due 
to spin, as we explained for Eq. (235). The fc-summation can be carried out 
easily as 

(379) 



v„ra 



V 



Y^S^k+q-a ~ £k,a)5(e k , a ) = ~f 9 st (q) 



2k Fa q 



where q = \q\ (9 s t(q) is defined in Eq. (220)). 



In contrast to the diagonal components of resistivity, which vanish in the 
adiabatic limit, the off-diagonal (i.e., Hall) components remain finite in this 
limit. In fact, we find 
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(We have retained the antisymmetric component, \{p S ij — Psji), in the second 
line.) This is the Hall effect caused by the spin chirality, or the spin Berry 
phase [108], demonstrated in the slowly varying case [109,110]. Comparing 
Eqs. (253) and (380), we see that the force in the adiabatic limit is exactly 
due to the Hall effect from spin chirality: 



2P 2 77 2 
F Hall = /e n y- ■ 

= e 2 nY,JjPsji- (381) 



15.5 Kubo formula 



When impurity scattering dominates the resistivity, we need to use the Kubo 
formula. The wall contribution to resistivity is calculated perturbatively The 
conductivity a for the current in the ^-direction is calculated from the imaginary- 
time current- current correlation function 

Q(t) = ^< T r J z (r) J z (0) >, (382) 

where T T denotes the time order in the imaginary time and the bracket rep- 
resents the thermal average; 

< T r J z (r)J z (0) >= Tr[pT T J z (r)J z (0)]. (383) 

Here p = e~^ n jZ , Z = Tre - ^, j3 = 1/(]cbT). The Fourier transform of Q(t) 
is written as 



Jo 

= — (J z (tu i )J z (-tu i )), (384) 

where uji = 2n£/P is the Bosonic Matsubara frequency and J z {iuj() = /q dre~ WlT J z {t) 
The Fourier transform of the electron is defined as a n = (l/y/]3) Jq dTe luJnT a(r), 
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u> n = (2n — l)im/ j3 being a Fermionic Matsubara frequency. Here, thermal fre- 
quencies are written as ug, u^,- ■ • or ou n , uj n >, • • • , where subscripts £, £',■■■ are 
used for Bosonic and n,n', ■ ■ ■ are for Fermionic frequencies. We consider in 
this section a static domain wall solution, given by 9q{z) and <fi = 0. In the 
imaginary time, the current is then written by use of Eq. (166) as 

J z (iu e ) = — Y, Y,( k zai, n +e a k,n -^J2 ^((iWk+q^+^^n • (385) 

171 k u>„ \ Z q ) 

The electron operator carrying a thermal frequency of u n + oog is denoted as 
a n+ £. The Kubo formula relates the correlation function (384) to the conduc- 
tivity as [106] 

a=lim-lm(Q(hu; + iO)-Q(iO)). (386) 

Here Q(ftuj + z'0) is the retarded correlation function obtained by analytical 
continuation, Q(hu> + iO) = Q(iu>e — > ftw + iQ), and Im denotes the imaginary 
part (iO denotes a infinitesimal imaginary part). We estimate the correction 
to the conductivity due to the wall to the second order of gauge field. In 
this section the classical (Boltzmann) contribution (denoted by er c ) is calcu- 
lated. The quantum corrections represented by maximally crossed diagrams 
are considered in the next section. 

As is well known, the zeroth-order term of Q is obtained as 



Qo(iue) 
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= ~\ YFq ^2 k z Gkn(rGk,n+e,a, (387) 

and this contribution leads by use of Eq. (386) to the classical conductivity 
due to the normal impurity, o"o = e 2 nT/m (n being the electron density) [106]. 
Here the imaginary time electron Green's function {G^na = — \ a kncr^kna)) 
includes the effect of the impurity and is given by [106] 

Gkna = —, — H T~V\ ' ( 388 ) 

i{uj n + ^sgnH) - e ka 

where sgn(n) = 1 and —1 for n > and n < 0, respectively. (The Green's 
function carrying a frequency of uj n + ujg is denoted by Gfc, n +v) 

Let us include the scattering by the domain wall. The first-order contribution 
of A q vanishes. The second-order contributions to the Boltzmann conductivity 
are shown in Fig. 36. The process Qi arises from the correction of both of the 
two current vertices by the wall, 5 J, and Q3 is due to the correction of one 
of the current vertices and a interaction with the wall. Q2 and Q4 are the 
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Fig. 36. The contributions to the Boltzmann conductivity in second order with 
respect to the domain wall. The interactions with the wall are denoted by wavy 
lines. The current vertex expressed by crosses represents Jo and the cross with 
wavy line represents 5 J. 



self-energy due to the wall and Q§ is the vertex correction. We here consider 



a domain wall with 



0, i.e., A. 



Ad„6)e y . The contributions in Fig. 36 



are written as (we write A y = A y z (q)) 
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(389) 



The contribution from the wall needs to vanish in the limit of the vanishing 
Zeeman splitting, M — > 0, since no scattering occurs there. This becomes 
obvious after summing the contribution Qi_5 = J2i=i 5 Qt- Using 
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(390) 



we obtain Q1-5 = Q c + Q' c , where [147] 



Q c (iue) 



and 
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The term Q c is dominant and the contribution from the term Q' c turns out to 
cancel with the effect of the shift of the electron density, which is calculated 
later. 

The summation over the Matsubara frequency, u> n , in Eqs. (391) and (392) 
can be carried out by use of contour integration (see §C) and the contribu- 
tion to the Boltzmann conductivity from the five classical processes, o"i_5 = 
lim w ^o im (Qi_5(^ + iO) — Q c (i0))/cu, is obtained as cri_ 5 = cr c + a' c , where 
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and a' c is the contribution from Q' c ; 
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Besides the processes in Fig. 36, there is another contribution to the classical 
conductivity, which is due to the change of the electron density in the presence 
of a domain wall [140,141]. (In this paper the chemical potential of the system 
is fixed as the domain wall is introduced, considering a constant voltage. In 
[140,141], on the other hand, the electron number is kept constant by shifting 
the chemical potential, which corresponds to a measurement under constant 
current. The difference is small if kpX ^> 1.) The correction to the electron 
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density due to the interaction with the gauge field is written as 
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This shift of the electron density leads to a correction of the zeroth-order 
conductivity, o"o — ► u + 5cr c , where 5a c = e 2 rSn/m is obtained from Eq. (396) 
as 



5a c 



h ( eh 
167T \ m 



E 






M-a 



(k z +q/2)k z 






+ UYM + a 



(k z +q/2)q 
2m 



(397) 



It turns out after the A;-summation that a' c + 5cr c vanishes in the case of 
kpX 3> 1 [147]. Thus the classical correction to the conductivity due to the 
wall is given simply by a c . 

To proceed further we neglect quantities of 0((q/kF) 2 ) and approximate €fc± (? /2,= FCT 
6^ ± [(h 2 k z q/2m) + aM]. This is because the momentum transfer, q, is lim- 
ited to a small value of q < A -1 due to the form factor of the wall, \A y n \ 2 oc 
[cosh(7rq , A/2)]~ 2 , and we are considering the case of a thick wall, k F \ 3> 1. 
The result of fc-summation is (see §C) 
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where x = 7iq\/2 and l a = hkFcrT/m is the mean free path of the electron 
with spin a. We have introduced a density of the wall, n w = 1/L (we have one 
wall in the sample of length L). The extension of the result to the many- wall 
case is straight-forward as long as the distance between walls is larger than A. 
Using this result the contribution of the wall to the Boltzmann resistivity is 
given as 



Pc = (o"o + c c )~ 



(Tn 



|gc| 



(399) 



The last expression is valid if the contribution from the impurity is much larger 
than that from the wall, namely if |cr c |/cr ^ 1. (But see also the discussion 
at the end of this section, before Eq. (404).) 

Here we consider the case of a ferromagnet with weak disorder, and assume 
two further conditions: 

Mr/h > 1, (400) 
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which indicates that the effect of the Zeeman splitting is not smeared by the 
width of energy level, and 

mMX/k F h 2 > 1. (401) 

The second condition is satisfied if the Zeeman splitting is not too small. Both 
inequalities would be satisfied in the case of d electron. In this case the classical 
correction due to the wall is obtained as (Eq. (C.19)) 
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(Mr/h > 1, m.MX/k F h 2 > 1). (402) 

Although the last equality of Eq. (399) is true only when |<x c |/o"o <C 1, the 
result of Eq. (399) with a c calculated from the Kubo formula (398)has a finite 
limiting value at r — ► oo of 

\a r \ m 2 A 2 ^ f°° dx 
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which is identical with the result of the Mori formula, (378), and thus is correct 
in this limit. Furthermore if we assume a finite lifetime for the electron in the 
Mori formula (374), which is without justification, we obtain 

, ,_ 2 h 3 fhA\ 2 l „, „ , 2 l , t jN 

P< = M ZT 77 E l^l 2 - 77TWT7 ~7^- ( 404 ) 

\ / kqa 

This is shown to be equivalent to Eqs. (398) and (399) by a similar calculation 
as in §C. These results may suggest that the expressions (398) and (399), which 
are justified only for |er c |/o"o <^ 1, are valid for any value of r. 



15.6 Spin chirality mechanism of Hall effect 



Here we will give an intuitive explanation for the Hall effect due to vorticity 
in terms of the Josephson effect of spin [178-180] . A persistent current in a 
metallic ring is an equilibrium current that can be induced when the time- 
reversal symmetry is broken. Such a current appears in the presence of a 
magnetic flux through a normal ring. The effect is due to a U(l) phase factor 
attached by the flux to the electron wave function. Here we show theoretically 
that a permanent current is induced in a conducting normal ring just by 
attaching three ferromagnets, without magnetic flux through the ring [178]. 
This surprising effect can be seen in nano-scales at low temperatures. The key 



116 



here is the non-commutativity of the SU(2) spin algebra, which breaks the 
time-reversal symmetry, and leads, in the presence of electron coherence, to a 
permanent electron current (spin Josephson effect [179]). We also show that 
this persistent current gives an intuitive explanation of the anomalous Hall 
effect due to the spin Berry phase in frustrated magnets [109,178]. 



15.6.1 Spin chirality and time reversal symmetry breaking 

The electron has spin 1/2 (i.e., has two components), and the spin obeys SU(2) 
algebra. The algebra is represented by three 2x2 Pauli matrices Oi ii = x, y, z) 
satisfying the commutation relation 

[a^aj] = 2ie ijk a k , (405) 

where e^ is the totally antisymmetric tensor with e xyz = 1. When a con- 
duction electron in a conductor is scattered by some magnetic object, the 
electron wave function is multiplied by an amplitude A{n) = ae lf3n ' cr = 
a(cos j3 + i(n ■ cr) sin/5), which is generally spin-dependent and is represented 
by a 2x2 matrix in spin space. Here a and (3 are complex numbers and n is 
a three-component unit vector representing the magnetization direction. We 
consider in this paper only classical, static scattering objects, and assume that 
the n's are constant vectors. 

™2 



(a) (b) 



Fig. 37. A closed path contributing to the amplitude of the electron propagation 
from x to x. At Aj, the electron experiences a scattering represented by an SU(2) 
amplitude, A(rij). The contributions from one path (a) and the reversed one (b) are 
different in general due to the non-commutativity of .A(raj)'s. 

Let us consider two successive scattering events represented by A(n{) and 
A(n 2 ) (Fig. 37). Due to the non-commutativity of <T;, the amplitude depends 
on the order of the scattering events: A{ni)A{n2) ^ A{n2)A(n{) in general. 
Various features in spin transport, which are under intensive pursuit recently, 
arise from this non-commutativity. Actually, 



A(ni)A(n 2 ) =a 2 (cos 2 (3 — (rti ■ n 2 ) sin 2 /? + isinfl cos(3(rii + n 2 ) ■ cr 

-zsin 2 /3(ni x n 2 ) ■ cr), (406) 

has an asymmetric part, 
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A(ni)A(n 2 ) - A(n 2 )A(n 1 ) = -2ia 2 sin 2 f3(rii x n 2 ) ■ cr. (407) 

This relation indicates that the spin current with polarization {n-y x n 2 ) flows 
between two spins, i.e., 



j: cx (m x n 2 )°e 12 , (408) 

where ei 2 represents a vector connecting two localized spins. This spin current 
is spontaneous (equilibrium), and represents an exchange interaction between 
rii and n 2 transmitted by the conduction electron [181]. In fact, when the spin 
current given by Eq. (408) flows, spin accumulation oc (ni x n 2 ) accumulates 
on rt\ and n 2 , which rotates n\ and n 2 within the rii-n 2 plane until the two 
spins become parallel or anti-parallel to each other. In the case of smooth spin 
structures, the equilibrium spin current, Eq. (408), reduces to [182] 



J s >(nxV,nf. (409) 

This non-commutativity of two spins does not affect the charge transport, 
since the charge is given as a sum of the two spin components (denoted by 
tr), and trL4(ni)A(n 2 )] — tr[A(n 2 )A(ni)] = 0. An anomaly in the charge 
transport arises at the third order. We have, by virtue of Eq.(405) and the 
relation tr [oiOj] = 25^, 

ti[A(n 1 )A(n 2 )A(n 3 )]-ti[A(n 3 )A(n 2 )A(n 1 )} = 4a 3 sin 3 pn r (n 2 xn 3 ) = iC 123 . 

(410) 
This relation indicates that in the presence of fixed rij's with rii ■ (n 2 x n 3 ) ^ 0, 
the symmetry under time-reversal (more appropriately, reversal of motion) is 
generally broken in the charge transport. In fact, the relation (410) indicates 
that the contribution from one path, x — > X\ — > X 2 — > X 3 — > x (Fig.37a), and 
its (time-) reversed one, x — > X 3 — > X 2 — > Xi — > x (Fig. 37b), are not equal, 
and this difference results in a spontaneous electron motion in a direction 
specified by the sign of C\ 23 , namely, a permanent current. What is essential 
here is the non-commutativity of the SU(2) algebra. In fact, C\ 23 vanishes if 
all rij's lie in a plane, in which case the algebra is reduced to a commutative 
U(l) algebra. The degree of the symmetry breaking, rii ■ (n 2 x n 3 ), is given 
by the non-coplanarity, often called spin chirality. 



15.6.2 Charge current 

The spontaneous current above would be realized on a small conducting ring 
with three ferromagnets or magnetic dots with different magnetization direc- 
tions, Si, S 2 , and S3 [178]. The electron in a ring feels an effective spin po- 
larization when it goes through the region (i^) affected by the ferromagnets, 
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and the effect will be modeled by the exchange (spin-dependent) potential, 
V(x) = —Mrii-er for x G F t . Here M represents the effective exchange field. 
The equilibrium charge current in the ring is calculated from 

He 
](x) = —Im(\7 x -V x/ )tiG(x,x',T = 0-)\ x , =x , (411) 

where G(x, x',r) = — < Tc(x,r)c f (x',0) > is the thermal Green's function, 
e, m, c being the charge, mass, and annihilation operator of electrons in the 
imaginary time, respectively. G(x,x',t) is calculated perturbatively from the 
Dyson equation, G = g + gVG, where g represents free Green function. As is 
seen from Eq. (410), a possible finite current arises at the third order in V. 
By summing the contribution of the two paths, x — > X\ —* X 2 — > A3 — ► x 
and the reversed one, we have j(x) = —-B(x) ReCi 23 . Here Ci 23 is defined 
by Eq. (410) with a = iM, (3 = vr/2, and™ 



P P (it ) 

B ( x ) = II L dXj —f(u) \/ Xo lm[9oi9i2923934}\ X4=Xo=x (412) 



describes the electron propagation through the ring, which is common to 
both paths. In Eq. (412), f(u) is the Fermi distribution function and gij = 
g r (Xi — Xj,u) is the retarded Green's function of free electrons. Approxi- 
mating the transport along the ring as one-dimensional and neglecting mul- 
tiple circulation, we have g r (x,tu) ~ —iir(D/L)e zkF > x >, where Jzf is the Fermi 
wavenumber, D the density of states (~ 1/ep ; e^ = h 2 k F /2m being the Fermi 
energy), and L the length of the ring perimeter. The final result is given by 



2e^ cos(k F L) (-pj S 1 -(S 2 xS 3 ), (413) 



at zero temperature. Here J = ttWM/L with W being the width of the 
ferromagnets, and vp = hkp/m is the Fermi velocity. 

The current is thus induced by the spin chirality Si ■ (S 2 x £3) of the fer- 
romagnets. This quantity reduces to the Pontryagin index (density) for the 
case of a smoothly varying field S(x), which is also interpreted as the Berry 
phase of the spin. The effect of the spin Berry phase on electron transport 
has so far been investigated in the limit of strong coupling to S(x) where the 
electron spin adiabatically follows S(x) [183]. In contrast, the present result 
Eq. (413) is obtained in the opposite limit; we have treated the coupling to 
S perturbatively (weak coupling) and made no assumption of smoothness on 
S(x). 

The appearance of the current is due to the symmetry breaking of the charge 
(U(l)) sector, as in the case of the current in Josephson junction. But note 
that here the U(l) symmetry breaking was due to the non-commutativity of 
spin (SU(2)) sector ("spin Josephson effect" [179]). 
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The system of a persistent current arising from the spin chirality works as a 
novel quantum operation gate, where spin Qbits and flux (current) Qbit are 
combined [179]. 



15.6.3 Anomalous Hall effect 

The phenomenon predicted here is not restricted to artificial nano-structures, 
but will be present rather generally in metallic frustrated spin systems such as 
pyrochlore ferromagnets [184,185] and spin glasses [186-190], where finite spin 
chirality is often realized. The spin chirality was recently pointed out in the 
adiabatic limit to induce a peculiar anomalous Hall effect [108,191,192]. The 
present chirality-driven persistent current [178] affords an intuitive interpreta- 
tion. The circulating current starts to drift when the electric field is applied, 
in the direction perpendicular to the electric field [25] (Fig. 38), just as in the 
normal Hall effect. With the frequency of the circulating motion, read from 
Eq. (413) as 

n~^(i.)V(S a x5 s ), (414) 

we may estimate the Hall conductivity by a xy = a^Vtr. Here o"o is the classical 
(Boltzmann) conductivity, r is the elastic lifetime, and the dirty case fir <C 1 
is assumed. If the spin chirality is located uniformly on every triangle of size 
of inter-atomic distance (i.e., S\ • (S 2 x S3) = Xo an d L ~ l/k F ), we have 
o~xy/o~o — Xo^ 3r / e F- This result agrees with the result based on the linear 
response theory [109,178]. 




Fig. 38. Schematic picture of the chirality-induced Hall effect. Circulating permanent 
currents are induced local around spin chirality, and this circular current drift in 
the perpendicular direction when an electric field is applied. 



When the spin structure is slowly varying, spin chirality can be expanded as 
S(x) • (S(aJi) x S(x 2 )) ~ X)(zi -x) M (x 2 -x) v S{x) ■ (d„S{x) x d v S(x)), (415) 

and the spin chirality reduces to the spin Berry phase [110]. 

So far, the chirality mechanism are confirmed in spin-glasses [186-190], but 
the case of Nd2Mo2C>7 seems to be under debate [192,193]. 
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16 Current from magnetization dynamics 



In most parts of this paper, we have studied magnetization dynamics induced 
by an electric current. The opposite effect, conversion of magnetic informa- 
tion into electric information, also exists, and has quite recent developments. 
The electromotive force (voltage) induced as a counter action of spin-transfer 
torque was in fact pointed out by Berger in 1986 [35]. Later, Stern showed 
in a general context that the fictitious field of the Berry phase ($s) accu- 
mulated by electron spin induces a spin-dependent electromotive force (V) 
via Faraday's law as V = — $5 [101]. This argument was applied to the case 
of domain wall dynamics by Barnes and Maekawa [194]. Duine [195] argued 
that spin relaxation (/3 sr ) induces a new current contribution, proportional 
to /3 ST (S ■ ViS), where S is a localized spin and i is the current direction. 
The effect of a Rashba-type spin-orbit interaction was studied by Ohe et al. 
[196,182] in the weak s-d coupling limit. The current in this case was found to 
be ji oc Xso^ijz \S x S), where A so is the strength of the spin-orbit interac- 
tion, z is direction of Rashba field, and (• • •) denotes averaging over electron 
motion. The quantity (SxS) represents a spin damping and is related phe- 
nomenologically to a spin current across the interface in the case of junctions 
[197,198]. 

Such a spin-dynamics-induced current was experimentally observed in systems 
where the spin-orbit interaction is strong. Saitoh et al. [199] demonstrated 
that when the spin-orbit interaction is strong, an inverse process of the spin 
Hall effect would occur. Namely, the spin current pumped by spin dynamnics 
[197,198] would be converted into charge current by spin-orbit interaction. 
They confirmed this idea by observing charge current in a Pt layer attached 
to a ferromagnet with dynamical magnetization. The induced current was per- 
pendicular to the spin current flow and (S x S), consistent with the inverse 
spin Hall mechanism, and with the theory of Refs. [196,182] at least at the 
qualitative level. Experimental confirmation of electromotive force due to the 
pure spin Berry phase without spin-orbit interaction is challenging and inter- 
esting. Closer studies on these phenomena will lead to a unified picture of the 
interplay between electric/spin current and magnetization. 



17 Summary 



We have reviewed the theoretical aspects of current-driven domain wall mo- 
tion, including microscopic derivation of the equation of motion, wall dynam- 
ics, and brief discussion of experimental results. The effect of current arises 
from the s-d exchange coupling between the localized spin and conduction 
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electrons. Treating the non-adiabaticity perturbatively, we have derived fully 
quantum mechanical expressions for torques and forces acting on the wall in 
terms of Green's functions. The obtained torques are summarized in Table 3. 
With these results, we derived the equation of motion of the wall. The wall is 
assumed to be rigid and planar (one-dimensional), described by two collective 
coordinates, position X and angle 0o of magnetization out of the easy plane. 
Spin-transfer torque reflecting angular-momentum conservation was shown 
to contribute to wall velocity, and spin relaxation and non-adiabaticity were 
shown to produce forces on the wall, which induce <J)q. Solving the equation of 
motion, we found that there is a threshold current for driving the wall arising 
from hard-axis magnetic anisotropy energy K± and/or extrinsic pinning po- 
tential Vq. The threshold current is determined by K± in the intrinsic pinning 
regime. In the extrinsic pinning regime, it is determined by Vo, Kj_, and the 
force from the current (represented by /3 W ). 

Table 3 

Summary of torques induced by electrons. The solid lines represent the conduction 
electron Green's functions, the wavy lines are the gauge fields representing the 
domain wall, the dotted lines are applied electric fields, and the double-dotted lines 
represent spin relaxation due to spin- flip scattering by magnetic impurities. The 
Green's functions (thick solid lines) in the contributions to the (3 terms contain 
lifetimes due to spin-flip impurities. 



local contributions 



torque 



sn 



spin renormalization 




js ■ Vn 



spin-transfer torque 



a sr n x n 



Gilbert damping 



<Xf GXC^Oc 



/3 sr (™ x (j s ■ V)n) 



(3 term 



non-local contributions 




nref 



reflection force 



Current-driven magnetization switching has an advantage in downsizing of the 
devices compared with the Ampere's mechanism. As we have seen above, the 
domain wall motion is governed by the applied current density and material 
parameters. In contrast, the field created by the Ampere's law is proportional 
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to the current, which is proportional to the area of the system. Therefore, the 
conventional Ampere's mechanism requires high current density and becomes 
inefficient in small devices, while the efficienty of the current-induced mech- 
anism is not reduced. The same applies to the current pumped by the spin 
dynamics due to the inverse spin Hall effect, discussed in §16. These two new 
magnetoelectric effects, based on quantum mechanical material features, have 
therefore great advantage in high density devices. 

Our predictions are briefly summarized as follows. 

The wall dynamics depends much on the behavior of </> . For j <C j a (j a defined 
in Eq. (279)), 0o remains small and there is no intrinsic pinning effect, while 
f° r j ^ Ja) 0o develops and an intrinsic pinning effect arises due to K±. In 
most cases, the intrinsic pinning threshold seems rather high, and the extrinsic 
pinning regime would be suitable for low-current operation. This is realized 
when (3 W is finite if the extrinsic pinning potential is much weaker than K± 
(since, then, <fro ~ and the wall does not suffer from the intrinsic pinning 
effect). In other words, wall motion at low current would be realized in systems 
with 

• very clean samples (with very few or weak pinning sites) and 

• large /3 W realized by large spin relaxation (/? sr ) or thin wall (/3 na ) 

Reduction of thresholds for thin walls seems consistent with experimental 
results [80,22]. 

The wall speed is in most cases governed by the spin-transfer rate (determined 
by angular momentum conservation), and is equal to v w ~ v s = t~§PJ- This 
speed is estimated to be v s = 100 m/s for j = 10 12 A/m , which is high 
enough (corresponding to switching at 10 GHz if the device is 10 nm size), 
but present experiments on metals are far below this limit. Faster walls beyond 
this spin-transfer limit could be realized if /? w /a is large (Eq. (286)). 

Although many studies, both theoretical and experimental, have been done 
so far, current- induced magnetization reversal has still exciting problems to 
be attacked. The final aim is to achieve fast switching at low current density. 
Roughly speaking, domain wall speeds need to be enhanced by a factor of 
10 (so that v w > 10 m/s) and the threshold current density needs to be 
reduced by factor of ^ (i.e., below 10 11 A/m ). To do this, material choice 
and structural design will be important. For instance, spin-orbit interaction 
is expected to lower the threshold current greatly by inducing effective force 
(/3 sr ), and this possibility should now be investigated quantitatively taking 
account of material details. For this, a microscopic formalism as described in 
the present paper needs to be combined with first-principle calculations. In 
addition, the effects of deformation and fluctuation (spin waves), the case of 
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vortex walls, and details of the extrinsic pinning potential need to be studied. 
Systematic studies for realizing efficient domain wall motion are now under 

way. 
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A Spin Lagrangian 



In this section, we derive the Lagrangian of a classical spin by employing the 
path integral for quantum spin. 



A.l Spin coherent state 



Consider a spin 1/2, which takes two states | j ) and | I ) pointing, respectively, 
parallel and antiparallel to the z-axis. The 2-axis here is not a special one, but 
could be chosen arbitrary. So one should be able to construct a spin eigen 
state | n ) pointing in an arbitrary direction, 

n = (sin#cos0, sin 9 sin <fr, cos6). (A.l) 

This state is by definition an eigen state, with eigenvalue +1, of the spin 
operator projected in the direction of n, and satisfies 



(n-cr)\n) = \n). (A.2) 

Here 



cos 9 e ^ sin 9 
na = | | (A.3) 

e J * sin 9 — cos 9 
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is the Pauli spin matrix projected to n. The solution of the eigenvalue equation 
(A. 2) is given by 



| n ) = |^0)=e-^ 2 cos^|T) + e^ 2 sin^||) (A.4) 

up to an overall phase factor. 1 The relative magnitude of the two coefficients 
is related to 6, and the relative phase to 0. The state (A.4) is called spin 
coherent state [92,200] 

Note that it does not satisfy the (three-component) vector identity, 



<r | n) ^ n I n), (A. 6) 

but only the (single-component) scalar identity, eq.(A.2). The expectation 
value, however, satisfies the vector identity, 



(n | a | n) = n. (A. 7) 



The following completeness relation (actually, they form an over-complete set) 
is easy to verify, 



/din 1 r w r^ir 

— \n)(n\ = — snxddd d<j)\6, <f>)(0, (f>\ 
2tt 2ir Jo Jo 

HTXTI + IIXII 

= 1. (A.8) 



The coherent state (A.4) for S = 1/2 is obtained from the state | f ) by two 
successive rotations as 

\0, <j))=e- i i a *e- i % a »\ T), (A.9) 



1 This state is mutivalued as a function of n (though single-valued as a function 
of 9 and 4>), so the notation | n) may not be appropriate. However, the quantity, 
| n)(n\, which will be used below, is single- valued. The same applies to the state, 

|n)' = |0,0)' = cos||T> + e^sui^U>, (A.5) 

which differs from (A.4) by an overall phase factor and is used more often. This is 
due to the spinor nature of the state, and does not mean any mathematical difficulty. 
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namely, first by 9 around the ?/-axis, and second by around the z-axis. For 
general S, the spin coherent state is given by 



\9, 0)s = e-^ 2 e-^|,S). (A.10) 

Here, S is a spin-5 operator, and the state \m ) is defined by S z \m ) = m\m). 
One can varify that ( n | S \ n ) = Sn. 



A . 2 Path integral for spin 



In the formulation of quantum mechanics in terms of path integral, a central 
role is played by the Lagrangian of the system. 

According to Feynman [103], the probability amplitude that a particle located 
at Xq at time t = will be found at later time t = T at position Xf is given 
by the sum of amplitudes over all possible paths x(t) satisfying x{t = 0) = x 
and x(t = T) = Xf as 



( Xi \e- iHT/h \x )= [ Dx(t)e iS[x( - t)]/h . (A.li; 



Here H is the Hamiltonian of the particle, and the «S[x(i)] is the action (time 
integral of Lagrangian). The same expression holds for spin. The probability 
amplitude that a spin in a state n at time t = will be found in a state rif 
at some later time t = T is given by the sum of amplitudes over all possible 
paths n(t) satisfying n(t = 0) = n and n(t = T) = rif as 

(m\ e- iHT/h \n ) = f Dn(t) e iS[n{t)]/h . (A.12) 

Here H is the Hamiltonian for the spin, and <S[n(i)] is the action of spin. Let 
us utilize this fact to calculate the Lagrangian for spin [92,200]. 

Consider the quantity on the left-hand side of eq. (A.12). By dividing the time 
interval T into many (N ^> 1) tiny intervals e = T/N, we first write it as 



(rif | e~ iHT \ n ) = (n f | e- iHe e~ im ■ ■ ■ e~ iHe \n ). (A. 13) 

(For simplicity, we drop h for the time being.) Next, we insert the completeness 
relation (A. 8) between each neighboring factors e~ lHe as 



-iHe„-iHe „-iHe | „ \ _. f dn N-l f dn 2 f dn X 



(n f |e |n ) - j -^— . . . J — J -^{m\e \n N -r) 



x (n 2 | e~ lH£ | ni ) (m | e- Jiie | n ) (A.14) 
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We have used the variable n^ for the k-th position counted from the right. We 
regard t\. = ke as a discretized time, and \rik) as a snapshot \n(tk)) of a path 
\n(t)) defined on a continuous time t. The k-th. element is written, neglecting 
terms of 0(e 2 ) and writing tk = t, as 



(n(i) | e~ im | n(i - e) ) = (n(t) |(1 - iHe) l\n(t) ) - £j\ n(t) > j + C(e 2 ) 

= 1 - eUn{t) \j t \ n(t) ) + i{n(t) \H\ n(t) )| + 0(e 2 ) 

= l+ieL[n(t)} + 0(e 2 ) 

= exp{ieL[n(t)]} + 0(e 2 ). (A.15) 



Here we have put 



d 



L[n(t)]=ih(n(t) \ — \ n(t) ) - (n(t) \H\ n(t) ). (A.16) 

(JjL 

(We have recovered h.) Therefore, we have 

(A.17) 
Taking the limit, A^ — > oo, e — > (with Ae = T = constant), the right-hand 
side is written as 

f Dn(t)expl^ f dt'L[n(t'))\, (A.18) 

expressing the sum over all possible trajectories n(t). 

The Lagrangian is thus given by L of eq.(A.16). The first term on the right- 
hand side is known as the Berry phase [201,202]. In terms of 9 and 0, we have 



L[n] = -<p cos 8-H(n). (A.19) 

Here H(n) = (n \H(&)\ n) is the Hamiltonian. Thanks to eq.(10), spin oper- 
ators in the (quantum) Hamiltonian can be replaced by the c-number coun- 
terparts. Lagrangian (A.19) is for S = 1/2. For general spin S, the first term 
of the Lagrangian becomes HS<p cos 8, as verified with eq.(A.lO). 



2 If we use the spin coherent state |n)' of (A. 5), the first term becomes | <fi (cos 6—1). 
The difference ^ (f> is the total time derivative, and does not affect the classical 
equation of motion. 
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A. 3 Lagrangian for classical spin 



In this subsection, we consider a classical spin of magnitude S. The Lagrangian 
is given by 



L s = hSJ> cos0 - H(n). (A.20) 

The first term on the right-hand side is the kinetic term, 3 which is essential 
to lead to the equation of motion for angular momentum. It has the same form 
as the spin Berry phase in quantum mechanics, but the concept of phase does 
not appear in classical mechanics considered here. (The appearance of h is due 
to the fact that we have written the angular momentum as hS with S being 
a dimensionless number.) Comparison with the Lagrangian L = pq — H for 
ordinary particles implies that HS cos 9 = hS z and are canonically conjugate 
each other. More precisely, it specifies (in quantum language) the commutation 
relation of spin. Let us see this within classical mechanics. 

Let us define the canonical momentum conjugate to <fi by 

dT 

P^ = — f = hS z , (A.21) 

d(p 

and the Poisson bracket [94] by 

dA dB dB dA 



{A,B} 



PB 



1 (dA dB dB dA 



(A.22) 



h\d<j>dS z d<pdS z ) 
It satisfies {0, HS z }pb = 1- Other components of spin, defined by S f± = S x ± 



iS y = J S 2 — S% e ±l ^, satisfies {hS z , HS^pb = Ti^S'^, and thus 

{hSi, hSj} PB = Y^ e ijkhS k . (A.23) 

k 

Namely, the SU(2) algebra of spin has been obtained within classical mechan- 
ics. It is notable that this information is contained in the kinetic term of the 
Lagrangian. 

§Al-4. Equation of Motion for Classical Spin 



This has been known as the "kinetic potential" [97]. 
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The equation of motion is obtained by taking the variation of the action, 



S = / L(n, h) dt. 



(A.24) 



Under an infinitesimal variation, n(t) — > n(t) + Sn(t), we have 

dL d dL 



5S 



dn dt dn 



■5ndt = 0. 



(A.25) 



Since \n\ = 1, all three components of Sn are not independent, and we cannot 
equate their coefficients to zero. Since 5n is orthogonal to n (in the first 
order in Sn), we put Sn = n x 5w, and regard all the components of 5w as 
independent. Then, from 



5S 



'dL d dL\ 
dn dt dn 



■ 5w dt = 0, 



we obtain the equation of motion as 

(d dL dL\ 
\ dt dn dn I 



(A.26) 



(A.27) 



Using 



5S 5S 



1 5S 



5n Ge 56 ^ sin 6 6(f) 
the differentiation of the kinetic term L = hS<p cos 6 is calculated as 



(A.28) 



d 8Ln dLr 



HS (e^d — e e <j)sm8 



dt dn dn 

= hS(h x n). 

The equation of motion (A.27) then becomes 



(A.29) 



n x hSh x n + - — = 0. 
V on I 



Using nn = 0, we obtain the well-known expression 



(A.30) 



h = 7-Heff x n - 
Here we have defined the effective field by 



(A.31; 
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1 8H 1 8H ,. . 

with 7 being the gyromagnetic ratio. 

The equation of motion, (A. 31), describes precessional motion around the ef- 
fective field, if c ff, as shown in Fig. 13 (left), but the effect of damping (energy- 
dissipation) is not included. To treat the damping in the Lagrangian formal- 
ism, we introduce the dissipation function [94], 



W = ~^ = ^hSn\ (A.33) 

where a is a dimensionless damping constant, and generalize eq.(A.27) to 

nx f-— -— — ^|=0 (A 341 

y dt dh dn dh I 

The equation of motion then becomes 



h = 7-fi c ff x n — an x h. (A. 35) 

The second term on the right-hand side represents the dapming torque, called 
Gilbert damping. 



A. 4 Landau-Lifshitz- Gilbert equation 

So far, we have essentially considered a single spin. In this subsection, we 
consider the dynamics of magnetization of ferromagnets consisting of many 
spins. 

For simplicity, we consider a localized model for ferromagnetism, and assume 
that the magnetization is carried by localized spins, Si, located at each site 
i. Writing as Si = Srti(t) where S is the magnitude of Si and rii is a unit 
vector representing its direction, we assume that S is constant, and only n,i(t) 
can vary. Moreover, we assume that its space/time variation is sufficiently 
smooth, and adopt a continuum approximation, rij(t) — > n(r,t), and write as 
9 = 9(r,t), <f) = <p(r,t). (Namely, it is a function of continuum space r and 
time t.) Then the magnetization (magnetic moment per unit volume) is given 
by 

M(r,t) = -jH £ Si(t) = -^fn(r,t). (A.36) 

unit volume 

(For simplicity, we consider a simple cubic lattice, and put the volume per 
localized spin to be a 3 .) 
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As a microscopic model, we have the following Hamiltonian in mind, 

H 8 = -J y £3 i -8 j + l y £ ( K ^ S ?v ~ KS l) ■ ( A ' 37 ) 

Here, J is the exchange coupling constant leading to the ferromagnetic state, 
K (> 0) and K± (> 0) are easy-axis and hard-axis anisotropy constants, re- 
spectively. (In soft magnets such as permalloy, shape magnetic anisotropy due 
to long-range dipole-dipole interaction is important. Here, we assume that 
such effects can also be modeled by this Hamiltonian. It is of course straight- 
forward to take such long-range interaction into the LLG equation, which can 
be treated numerically, but its analytical treatment is almost impossible.) 

In the continuum approximation, the Hamiltonian is expressed as 

S 2 rdh 



Hs ~YJ a 3 



J(\7n) 2 + K ± n 2 - Kn\ 



(A.38) 



Here J = Ja 2 with a being a lattice constant . Similarly, the Lagrangian L$ 
and the dissipation function Ws becomes 

L s = HS f — f cos 9 - H s (A.39) 

J a 6 

W s = ^hsf^n 2 . (A.40) 

2 J a 6 

From these, the equation of motion is obtained as 

h = 7-Heff x n — an x h (A. 41) 

(Usually, this equation is written in terms of magnetization vector, M , rather 
than the spin direction, n.) Here, 

„3 stj a 

7 /feff = ^r^f = -r [-JV 2 n - Kn z z + K ± n y y] + 1 H cxt (A.42) 

is the effective field coming from Hs, 4 and we have added an external field, 
Hext- The equation (A. 41) is formally the same as eq.(A.35), but here n = 
n(r, t) is a field variable. 

A phenomenological equation describing the space/time variation of magneti- 
zation vector of ferromagnets was first introduced by Landau and Lifshitz in 



4 Here, the derivative means (spatial) functional derivative, and an additional factor 
a 3 appears compared to eq.(A.32). Mathematically, it comes from the factor 1/a 3 
ineqs.(A.39) and (A.40). 
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1935. They used an equation of the form, 

n = 7'Heff x n + a LL ^'H cS x n) x n, (A.43) 

to study the structure and dynamics of a domain wall, etc. The second tern on 
the right-hand side represents damping torque, called Landau-Lifshitz damp- 
ing. The Landau-Lifshitz (LL) equation (A.43) and the Landau-Lifshitz-Gilbert 
(LLG) equation (A. 41) are equivalent, and their coefficients are related by 
«ll = ol and 7' = 7/(1 + a 2 ). 

§Al-6. Effects of Conduction Electrons 

To study magnetization dynamics driven by electric/spin current, let us in- 
troduce conduction electrons. We treat conduction electrons (s electrons) and 
magetization (d spin) as distinct degrees of freedom (the s-d model), which 
are coupled by the s-d exchange interaction, 

H sd = -M fd 3 xn(r,t)-(J(Tc). (A.44) 

Here c^ and c are electron creation/annihilation operators (spinors), and (c^crc) 
represents spin density. (We have put M = J s dS with J sd being the s-d ex- 
change coupling constant.) Through H sd , the electrons exerts an effective field 

„, a 3 / 5H sd \ a 3 „,. + . ,.„,-% 

on magnetization, and thus the torque (times a 3 /HS) 

t' cl = jH' eS xn = ^Mn x (cW). (A.46) 

Here (• ■ ■ ) represents quantum-mechanical and statistical average. The LLG 
equation is then becomes 

h = ^H cS x n — an x h + t' cl . (A. 47) 



B Brief introduction to non-equilibrium Green function 

We consider the electron system evoluving under the Hamiltonian 

H(t) = H + H c (t), (B.l) 

where H e (t) represents the interaction Hamitonian driven by an external field 
which is switchied on at a certain time, to, much earlier than the time we are 
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to t\ p t'l t 

Kjtot "i 1 ! V 

4 ! d 4 ^l 



Cd 

Fig. B.l. The contour Ct Q t in the case of t\ < t^. Times t\ and t'i are defined on the 
contour Ct t- The contour runs on the complex-time plane from to on the slightly 
upper branch (C u ) to t on the real axis, and back again from t to to on the slightly 
downward branch (Cd). 

interested. This interaction is treated perturbatively. The Hamitonian H can 
contain interaction term which are also treated perturbatively. It is written 
therefore as H = Hq + H 1 , where Hq is perturbed part and interaction term to 
be treated perturbatively is H 1 . (The differece between H e (t) and H 1 is that 
while H 1 is there at t < to, H e (t) vanishs.) For t < t , the system is in the 
thermal equilibrium state determined by H and the expectation value of an 
arbitrary operator O is calculated by the density matrix 



e -/3(H-pN) 

Tr(e-^ H -' iN ^) 



P( H ) = ^,- m -,N^ (B-2) 



as (O)h = Tr (p(H)O). Here \x and N are respectively the chemical potential 
and the total number operator of the system. 

We introduce the following Green's function [100,203,204] 5 denned on closed- 
time pass contour C ioi as depicted in Fig. B.l. 

G(x,x') = -i(Tc t0t [a n (x)4 i (x')]) H , (B.3) 

where a n = (a n , , a n _) is the two-component electron operator in the Heisen- 
berg picture and T Ct t is the contour-ordering operator. We here put x = (x,ti) 
and x' = (x 1 , t^) for simplicity. Since there is four different combinations for the 
times t\ and t' x located on either of the two branches of C u and Cd, the contour- 
ordered Green's function G(x,x') contains four different functions (Fig. B.2): 



G(x,x') 



G c (x,x'), tx, t[ eC u 
G>(x,x'), hed, t[eC u 
G<(x,x'), tiGC u , tieC d 



G- C {x,x'), t u t[ GC d 
Here G c is the causal, or time-ordered Green's function and G 5 is the antitime- 



In this subsection, we put h = 1. 
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to tl tl t to tl tl t 

g c = ; I ; ^ U G> =t | i ; ^ 

to tl tl t to tl tl t 

g = . j « ; y~^ g^c — * F . . . y 



Fig. B.2. Four different Green's functions 
ordered Green's function, which are given by 

G c (x, x') = -i(T[a n (x)al i (x')]) H 

= -iO(ti - t' 1 )(a n (x)al i (x')) H + id(t'i - t 1 )(al i (x)a n (x')) H ,(B.5) 

Gc(x, x') = -i{f[a n (x)a\ i (x)]) H 

\ f-iw i -/i/j. j./\/_.t c_\_ /'_/\\ 

(B.6) 



-iQty'i - ti){a n (x)a} H (x)) H + i9(h - t' l )(a\ i (x)a n {x)) 1 



where T(T) is the usual (anti)time-ordering operator and 6(t) is the step 
function. The second G > and the third G K in Eq. (B.4) are respectively the 
greater and the lesser Green's functions, which are given by 

G > (x,x') = -i{a H {x)a\ i {x')) H , (B.7) 

G < (x,x')=i(al i (x')a n (x)) H . (B.8) 

These four functions are not independent of each other and related as 

G c (x, x) + G- C (x, x) = G < (x, x) + G>{x, x'). (B.9) 

Furthermore, the retareded and the advanced Green's functions are written 
by using G > and G K as 

G R (x,x') = -i9(t 1 -t' 1 ){{a n (x),al i (x , )}) H 

= 9{t x - t'JiG^x, x') - G<(x, x% (B.10) 

G A (x, x') = iQ(t[ - ti)({a w (x), a ] H (x')}) H 

= 6(t' l - t 1 )[G < (a;,x , ) - G>(x,x% (B.ll) 

We define the Fourier transform of the contour-ordered Green's function as 

G^(x,t u x>,t[) = e r t- r ^^-^'^^'^G^A^m 

, , / J -co Z7T J — oo 2i7T 
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Note that we here explicitly write down the spin indices a, a'. In the free 
Hamiltonian H = H = I]fc,o- £fco- a Lr a fc<x and the zero temperature, G > and 
G < are given by 

G ka,k>a>{u,u') = 2tv5(lo - uj')5 k ^5 a y g^ a (uj) , (B.13) 

g> a {u) = -2vrz(l - f(u))6(u - &,), (B.14) 

G L,fc'a'(^> ^') = 27r5 (^ - ^/)<*fe,fe' ( Wfl , £r(^), (B.15) 

^(w) = 27ri/(a;)«J(w-e*«r), (B.16) 

where g k(T {oj) and g ka {oj) are respectively the free greater and lesser Green's 
functions in the momentum space, f{uj) = (1 + e^)" 1 is the Fermi distribution 
function and £fco- = £k — ^Fo- with e Fcr = e F + aM. Using these and taking ac- 
count of (B.10) and (B.ll), we obtain the free retarded and advanced Green's 
functions in the momentum space as 



<?L» = ^—-tt, (B.17) 

<£» = (£(*, <"))*, (B.18) 

where is the positive infinitesimal coming from the step function 9 in Eqs. 
(B.10) and (B.ll). 



B.l Perturbation expansion for closed-time ordered Green's function 



In this subsection, we derive the general formula of the perturbation expansion 
for the contour-ordered Green's function. 

We first rewrite an,k(t) in the Heisenberg picture as 



Ow )fc (t) = Ulj (t,t )a h (t)U Ho (t,to), (B.19) 

where dk(t) is the electron operator in the interaction picture and UH (t,to) 
is the unitary operator defined by 



U Ho {tM = Te- l ^ dt ' {H ^ t ' )+H ^ t ' ) \ (B.20) 

with 



H [ Ho (t) = e «*o(t-*o)jyi( t ) e -*tf (t-to) 

H e H (t) = e iHo{t - to) H c (t)e- iHo{t - to) (B.21) 
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ci- 

do —i/3 
Fig. B.3. The imaginary time contour C\ . 

in the interaction picture with respect to the free Hamiltonian H . Note that 
we have chosen to as reference time so that two pictures coincide. Using the 
contour-ordering operator Tc t t , Equation (B.19) is rewritten as [203] 



anAt) = T c to t e 



-k ot ^v^v')\ k{t) 



(B.22) 



We thus obtain the contour-ordered Green' function 



Gk,k'(tl,tl) — —i(Tc t()t [0'H,k(tl-) a 'H,k'(ti)])H 



-i(T, 



Ctn 



-^^^(O+fiVO} 



"afc(<i)oL(<i) 



. (B.23) 



J / H 



The next step is to rewrite the density matrix p(H) as 

p(H) = e- /3 (^ 0+ ^ i ^^ Ar )/Tr(e" /3( ^ 0+ ^ i ~ AiiV) ) 

= e-K H °-^U Ho (to,t - iP)/Tr(e-K H °^U Ho (t ,t - i/3)), (B.24) 

where 



U Ho (to,t - i(3) = T c; e - i £r 0dt ' H "° {t ' ) 



(B.25) 



is defined on the imaginary time contour C\ as depicted in Fig. B.3. Combining 
with the following relation 



Trv , e 



■■'()' 



-i f dt'H'„ ft') -t f dt'J?« (t')\ 



(B.26) 



we thus obtain Gfe,fe'(^i)^i) i n the form 
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Cut 



u 

i— 



tl 



ti 



t 



• to—i/3 

Fig. B.4. The contour C\ t in the case of t\ < t^. 
-CO tl t'l t 

ct i i ' 



> 



Fig. B.5. The contour Cj. 



Gk,k'(ti,txj 



-iTr e-^-^Tci 



-i/i d?H\j(?) -if dt 'H% (?) , . , .' 



Tr e-«^-" w )T r , 



C *o* I e 



4 r = di'j?' (f) _if di /^o (f) 



-i T c i 



-*/ci d *'^L(*') -«L dt'H% (?) , s , ,, ' 
e Jc v H ° l a fc (*i)aU*i) 



;b.27) 



J/ H 



where the contour Ci t is depicted in Fig. B.4. In this paper, we do not consider 
initial correlations and therefore put to = — °o to ignore the contribution from 
the imaginary part of contour C\ Qt , namely, C\ ot —> C to t [204]. In this prescrip- 
tion, we can expand the exponents of Eq. (B.27) with respect to H l Ho (t)and 
Hfj (t) on equal footing and can use the Langreth method, which is provided 
in the next subsection. Finally we rewrite Eq. (B.27) as 



Gk,k'(ti,ti] 



i(T, 



c, 



-i f d?H\, (?) -i f d?H% (?) ,, , + , ,, , 

Jc t H y > e j Ct H y > a k (ti)a} k (t\) 



0-28) 

H 



where the contour C t is depicted in Fig. B.5. 



If we furthermore put t = oo in the contour C t , this is often called the Keldysh 
contour C K as depicted in Fig. B.6 and the corresponding Green's function is 
called the Keldysh Green's function [100,205]. The Keldysh Green's function 
consists of G r , G a and G K = G < + G < , which is fundamentally equivalent to 
the contour-ordered Green's function. In this paper, since our main evaluation 
is to calculate the lesser Green's function G < , we do not use the Keldysh 
representation. 
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Fig. B.6. The Keldysh contour C . 
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Fig. B.7. The contour Ct changed into C^ + C t > . The both contour are equivalent. 



B.2 Langreth method 



In this subsection, we briefly review the Langreth method [203,204,206], which 
is useful to evaluate perterbation expansion of the contour-ordered Green's 
function in Eq. (B.28). 

In a concrete calculation for perturbative expansion of the contour-ordered 
Green's function in Eq. (B.28), we need to take the lesser component of the 
following integration 



C(t 1 ,t' 1 )= / dTA(i 1 ,r)B(r,i' 1 ) 

JCt 



(B.29) 



where C, A, and B is the contour-ordered Green's functions, which are defined 
on the contour C* in Fig. B.5. Taking the lesser component of C(ti,ti) and 
changing the contour C t into C tl + C t > as depicted in Fig. B.7, the lesser 
component C < is given in the form 



C < (t 1 ,t[)= f dr A(t 1 ,r)B < (r,f 1 )+ f dr A < (t 1 , t)B(t, t[). (B.30) 



Futhermore the first term in the right hand side of the above expression is 
written as 
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r /•*! /■— oo 

/ dr A( y t 1 ,r)B < ( y r,t' 1 ) = dr A > (t u r)B < ( y T,t' 1 ) +/ rfr A<(t 1; r) J B < (r, t;; 



r*i 

dr9(t 1 -r){A>(t l ,T)B<(r,t' 1 )-A < (t 1 ,r)B < (T,t' 1 )} 



OO 



OO 



/CO 
rfr A r (t 1 ,r)5 < (r,t / 1 ). (B.31) 

-co 

In a similar manner for the second term in Eq. (B.30), we thus obtain 

/oo 
dr (A T (t 1 ,r)B < (r,t[)+A < (h,r)B^r,t' 1 )). (B.32) 

-co 

For the retarded and the advanced Green's functions, we obtain 

/oo 
dTA r (ti,r)S r (T,ii) 
-oo 
/oo 
dr A a (ti,r)5 a (r,ti). (B.33) 

-oo 

These formula are abbreviated as 



[AB] K = A T B < + A < B a , [AB] r = A T B r [AB] a = A*B*. (B.34) 

By using these, equations (B.32) and (B.33), and the Fourier components 
C v (uj,lj') are written as 



/oo 
d<rLA(t 1 ,T)£(T,t / 1 )]' ? , (B.35) 

-oo 

/°° d(, i" 
^-[A{u,u")B(u",u')] T >, (B.36) 

-oo 2tt 

where t] is denoted by 77 =<,r,a. The result (B.34) can be extended for a 
product of three or more functions as 



[ABC]< = [AB] r C< + [AB^C* 

= A'B'C* + A V B < C & + A<B a C a , (B.37) 

[ABCf = A T B r C T 
[ABCf = A & B & C*. (B.38) 

B. 3 Impurity averaged Green's function 



In this subsection, using the paturbation expansion of Eq. (B.28) and the 
Langreth method, we calculate the contour-ordered Green's function averag- 



139 



ing over the rondomly distributed impurities without spin-flip scattering. The 
Hamiltonian is given by 



Hq — 2^, e kcrO'ka a ka, 



k.a 



H l = if imp = Y^ V q a k+qa a ka, V q = U\ ^ C ~ 



iq-Ri 



ka 



j=l 



H e (t) = 0, 



(B.39) 

(B.40) 
(B.41) 



where we here consider the short ranged potential V(r) = U[5(r — Rj) with 
Ri being position of impurity. Since we can regard the averaged V q in terms 
of position of impurity as zero, the first order Green's function ^fejfev^i,^) 
with respect to H 1 is zero. The next is the second order Green's function 

#i2,fc'<r'(*i>*i)> which is S iven b y 

£,w(ti,()=VEWWj[ dt 2 f dt'J^-t^ih-Qgtlih-A) 

p JCt JC't 

= WW / dt 2 f dt'J^ih - t 2 )U a (t 2 - t' 2 )g£l(t 2 - t[), (B.42) 
where g k l is the free contour-ordered Green's function and 



n,(t) = ^Eft 



(B.43) 



is the self energy in the first Born approximation. We here have used the 
relation Vk- p V p - q ' = rtiufS^y, where n\ is the concentration of impurity. The 
diagram representing g k l{uj) is shown in Fig. B.8(a). 

Using the Langreth method (B.35), we can write 

dh / dt> 2 [g^ih - t 2 )U a (t 2 - t' 2 )g k %(t 2 - t[ 

-oo J — oo L 



(B.44) 



and this Fourier component is 

\9wAu, u')] T = 2tt5(cj - uj^Sk^yg^iuj), 



(B.45) 
(B.46) 



The diagram representing g k J.(u)) is shown in Fig. B.8(b). In a simmilar man- 
ner, we can calculate the higher order Green's functions and sum the diagrams 
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(a)0<°> = > (b) 5 (2) 



(c) 9 = 



Fig. B.8. Diagrams of the Green's functions. (a)The free Green's function. (b)The 
second order Green's function. A dot line represents the impurity potential and 
two dot lines with one vertex represents averaging over impurity positions. (c)The 
Green's function averaging over the randomly distributed impurities. 



as dipicted in Fig. B.8(c), flw(o>) satisfies the Dyson equation 

9U«>) = *£(") + 02(w)lMwW«). (B.47) 

Using the Langreth method, we obtain 



(C M) _1 - n r » ^ - £fe + £ FCT H- z 7ct 

(B.48) 



where 27 CT = l/r CT is the damping rate, which is coming from 

n^.(o;) ~ —vKn-^Vo = —i^a (B.49) 

with v a being the density of state at e-p a . 

Similarly, taking account of (B.48), the lesser component oi g^iuo) is given by 



(0)< , (OW < 
9ka ^ 9ka LL a9ka 

/(u;)(rt»-y£»)- (B.50) 



„(0)< , „(0)r n r n < , „(0)r n < a , _(0)< n a _a 
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C Details of linear response calculation 

C. 1 Summation over Matsubara frequency 



The summation over Matsubara frequency, u n , in the expression of classical 
contribution to the conductivity, Eq. (391), is carried out in a standard manner 
by use of contour integration as follows. The quantity we consider is 

Q e {iu t ) = -B^£ |^(g)| 2 /^(*u;«), (CI) 

kqa 

where B = {{eh) 2 /2m 2 ) and 

Ikqa{iut) = — } j G k _g na G h _^ n+t<T G h+ ^ n _ (T G h+ ^ n+ i_ (T . (C.2) 

P n 

This function is written by use of contour integration with respect to z = iuj n 

as 

Ik q *{iui) = -^—f dzf{z)G k _i jff {z)G k _i j(T {z+iuJe)G k+ i_ (T {z)G k+ i_ tT {z+iuj i ) 

(C.3) 
where f{z) = 1/(1 + e 13 *) and Green's function G kcm is here denoted by 
G ka {iu n ), and the contour Co goes around the imaginary axis in complex 
z-plane (Fig. C.l). Noting that G ka (z) — [z+ (l/2r)sgn(Im[z]) — e ka \~ l has a 
cut along Im[z] = 0, the contour can be changed to four paths C\ parallel to 
the real axis, namely z = ±(u/ + zO) and z = — iuJi ± {uj 1 + zO), where uj 1 runs 
from — oo to oo (Fig. C.l). We then obtain 



1 f°° 
Ikqa{iue) = / du> 

IX J-oo 



xlm 



f , t) 1 1 

(u/ + iu t + £- e k ^ a ) {uj' + iui + §- e fe+ a _ CT ) 
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1 
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' + £ - €*_!,„) (a/ + f - e fe+f ,_,) 
+/(u/ - ^) — — m — - 



xlm 



1 



w 



/ _|_ * 



ft 

2r 



e fc-f ,<t) (^ + 27 - e fc+f ,-er; 



(C.4) 



The classical conductivity is expressed by taking the imaginary part of the 
correlation function analytically continued to uj + iO = iug as 

. mIm Q £ ^ ±i 0)-Q £M 
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tx 



Fig. C.l. The contour Co (denoted by two dotted lines) of integration in complex 
z-plane which appears in rewriting the summation over the Matsubara frequencies 
(Eq. (C.3)), which surrounds the imaginary axis. Because of the function f(z) in 
the integrand, there are poles on the imaginary axis at z = (2n — l)iri//3 where 
n = 0, ±1, ±2 • • • . C\ denoted by four solid lines is a deformed contour parallel to 
the real axis. The imaginary part of two lines are z = =ti0 and the other two are 
z = —ioJi ± iO, where uji (denoted by a thick cross with circle) is a pole and iO 
represents small imaginary part. 

The imaginary part of If, qa (iuji = uo + iO) is obtained from (C.4) as 
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(C.6) 



By use of 



/(a/) - /(a/ + w) 



UCo> 



(C.7) 
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which holds at small uj and low temperature, we obtain 

1 1 



lml kqa (uj + z'O) 
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Im 



ih ih 

2 r fe -f ,C 2t fc fc+|,-CT 



(C. 



The classical contribution to the conductivity, (C.5), is then obtained as 



A 2 h 3 fehY 1 ^. , l2 



(e fc _| iCT + e k+ 



57YT* 



ml V 
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(C.9) 



C.2 Summation over wave vector k 



The fc-summation in (393) is carried out as follows. We neglect quantities of 
0((q/k F ) 2 ) and approximate e k ± q /2^a — £k±[(k z q/2m)+aA}. This is because 
the momentum transfer, q, is limited to a small value of q < A" 1 due to the 
form factor of the wall, \A q \ 2 oc [cosh(7rgA/2)]~ 2 , and we are considering the 
case of a thick wall, kp\ ^> 1. Then <r w is written as 



A 2 h 3 (eh 



Ur 



SlXT- 



g) hv A ^ 
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This function is written as 
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(C.12) 



where iV(e) = (Vm 3 ^ 2 /7r 2 ^/2h 3 )y/e + e^ is the density of states, and fc 



2m(e + €f)/H, 6 being the angle between k and the z-axis, and Fg(e) is 
defined by the last line. In the first line we have included the factor of two 
due to the summation over the spin. The integration over e can be carried out 
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by use of contour integration. In doing this we need to be a little bit careful 
since the density of states in three-dimensions N(e) oc ^/e + ep has a cut on 
the real axis running from e = — €f to e = oo. Choosing a closed path C% in 
the e-plane as in Fig. C.2, the e-integral in (C.12) is written as 

f" deN(e)F 9 (e) = \ <f deN(e)F 9 (e). (C.13) 

J — ep <L JC2 

The contour C2 contains four poles at e = e* >± = a f A -\ 2rn cos ^J ^ z if • 
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Fig. C.2. The contour of integration C2 in complex e-plane. There is a cut on the 
real axis from e = — cf to +00, because of the behavior of the density of states, 
iV(e) oc (e + e^) 1 ' 2 in three-dimensions. 



In ferromagnets with strong polarization we are interested in, A ~ 0(ep), and 
then neglecting 0{q/kp) contributions we may approximate e* ± ~ a (A 



^f 



i-^p, where hkp a = y2m(eF + crA) is the Fermi momentum of the polarised 
electron. The density of states estimated at the pole in the lower half plane 
(i.e., at e*_ ~ e 2,n e* + ) has a opposite sign as the upper half plane, i.e., 
N{e%_) ~ e ni N(el + ) = -N a (N a = N(e* a+ ) = N(a£)). We therefore obtain 



/OO / n- \ 3 

rfeiV(e)F e (e) ~ 2vr ^-J £ JV ff 



1 



(e^Lcose + vA) + (£) 



2 • 



(C.14) 



where N a = (VmkF a /27r 2 h 2 ). The integration over cos6 in eq. (C.12) is then 
easily carried out to obtain 



2m 2 Vt a _ 1 
Jq ~ rfP \ q 



tan^(^ + A 
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_ t 2r fh 2 k Fa q 

tan — 

n \ 2m 
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(C.15) 
The conductivity (CIO) is then obtained by use of the expression of A v z (q) 



(EJ^(g)P 



(tt/L) / rfg[cosh 2 (7rgA/2)]- 1 • • • ) as 
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(C.16) 



Changing the variable x = irXq/2 we obtain (398) 



To look into the asymptotic behaviors at large A it is useful to rewrite the 
result by use of tan" 1 x = n/2 — tan" 1 (l/x) for x > and tan" 1 x = —n/2 — 
tan _1 (l/x) for x < Oas 



Or 



Airh 



AVE 



7T 



00 dx 1 



00 dx 



tan 1 
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/a ct x cosh x Jo x cosh x (2A) 2 - (Z CT x) 2 + 1 

where A = At/H, l a = (2l a /ii\), and A a = 2A/l a = (iimA\/k Fa h 2 ). We 
consider the case A>1 and A CT ^> 1, which would be satisfied for d electrons 
in 3d transition metals. Then the conductivity correction is approximated as 
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(C.18) 



In the case of a thick wall (A ^> k F l ) with a finite r, the first term in (C.18) 
is exponentially small and thus neglected. The conductivity in this case is 



<y c 



e ^ ~ r°° dx 



8nh 



cosh x 



8nH 



-n. 



J2 Z~ (Ar/ft, mA\/k Fa h 2 > 1). 

(C.19) 



On the other hand if we take the limit of r —* 00 first, the first term in 
(C.18) becomes dominant and the result of the Mori formula is obtained [147]. 
The result of Cabrera and Falicov[137], obtained by calculating the reflection 
coefficient in the absence of impurities, corresponds to this limit. 
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